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Abstract 


Lnndform  morphometry  is  related  to  causative 
factors  of  climate,  mantle  characteristics,  vegeta¬ 
tion  density,  and  lithology.  Drainage  basins  ana¬ 
lyzed  are  in  nature  fluvial  development,  free  from 
obvious  structural  influence,  and  in  climates  rang¬ 
ing  from  arid  to  humid.  Cher  80  basins  in  Arizona, 
Colorado,  New  Mexico,  and  Utah  were  inspected  in 
the  field;  22  were  subjected  to  detailed  field  inves¬ 
tigations. 

Morphometric  properties  of  all  basins  were 
measured  from  recent  1:  24,000  topographic  maps, 
or  special  drainage  maps,  and  show  great  ranges: 
valley-side  slope  has  a  four-fold  range,  drainage 
density  and  basin  relief  about  100-foid,  and  basin 
area  5,000-fold. 

Climate  was  measured  by  the  Thornthwaite 
precipitation-effectiveness  (P-E)  index,  and  the  in¬ 
tensity  of  average  five-year,  one-hour  rains.  Lith¬ 
ologies  distinguished  are  shale,  clastic  rocks, 
schist,  limestone,  granitic  rocks,  and  acid  volcanic 
rocks,  assumed  to  indicate  roughly  the  proportion 
of  clay  in  residual  soils.  Mantle  and  cover  charac¬ 
teristics  are  summarized  by  relative  infiltration 
capacity,  wet  and  dry  soil  strength,  percentage  ot 
bare  area,  and  surface  roughness. 

Validity  of  statistical  inferences  is  affected  by 
sampling  procedures,  which  were  mainly  non- 
random.  Analysis-of-variance  tests  show  that 
shale  and  schist  produce  high  drainage  densities, 
whereas  granitic,  clastic,  and  acid  volcanic  rocks 
produce  low  densities,  with  no  significant  differ¬ 
ences  among  these  three.  Shale  and  schist  have 
low  infiltration  capacity,  whereas  clastic  and  gran¬ 
itic  rocks  produce  high  values.  Valley-side  slopes 
are  greatest  in  shale  and  limestone,  and  lowest  in 
schist,  on  the  average,  subject  to  differences  in 
erosional  history. 

Drainage  density  averaged  highest  in  basins 
having  gullied  slopes,  intermediate  in  unmodified 
basins,  and  lowest  in  basins  with  trenched  main 
channels.  Differences  were  not  judged  significant, 
even  with  adjustment  for  infiltration,  though  low 
power  renders  the  test  inconclusive. 


Using  punched  cards,  variables  were  analyzed 
m  pairs  for  significant  correlations  and  regres¬ 
sions.  Texture  of  topography  decreases  exponenti¬ 
ally  with  valley-side  slope.  Relative  relief  and 
ruggedness  number  corrected  for  order  increase 
exponentially  with  slopes.  Texture  questionably 
varies  inversely  with  relative  relief. 

High  correlations  exist  among  channel  fre¬ 
quency,  drainage  density,  and  frequency  and  den¬ 
sity  of  first-order  channels.  It  is  inferred  that 
average  channel-segment  length  of  any  order  var¬ 
ies  inversely  with  the  total  channel  length,  for 
areas  of  the  same  size  but  differing  drainage  den¬ 
sities.  Average  first-order  channel  length  is  a  con¬ 
stant  proportion  of  average  channel  length  of  all 
orders,  and  is  inversely  proportional  to  total  chan¬ 
nel  length. 

Multiple  correlation  analyses  show  that  valley- 
side  slopes  and  drainage  density  ace  related  to 
climate,  and  to  properties  of  mantle  and  vegetal 
cover,  as  measured  in  the  field.  Slopes  increase 
with  infiltration  capacity  and  P-E  index,  but  vary 
inversely  with  wet  soil  strength  and  runoff  inten¬ 
sity-frequency.  Greater  runoff  produces  lower 
valley-side  slopes.  Climate  enters  further  by  the 
increase  of  infiltration  capacity  with  P-E  index. 
Greater  infiltration  permits  interflow,  reducing 
the  ratio  of  slope  runoff  to  channel  flow.  Channels 
are  deepened  with  less  opportunity  for  erosion  of 
slopes,  producing  steep  slopes. 

Drainage  density  varies  with  percent  bare 
area  and  runoff  frequency-intensity,  but  inversely 
with  P-E  index  and  infiltration  capacity,  fitting 
theoretical  observations  of  Gilbert,  Davis,  and 
Johnson,  and  confirming  Horton’s  infiltration  the¬ 
ory  of  erosion.  Primary  controls  of  topographic 
texture  are  lithology  and  climate,  which  act 
through  the  agenc'  of  secondary  surface  proper¬ 
ties.  Accidental  changes  in  these  secondary  con¬ 
trols,  unrelated  to  climatic  changes,  will  quickly 
alter  the  drainage  density  to  another  value  than 
that  expected  from  consideration  of  climate  and 
lithology  alone. 
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Introduction  and  Acknowledgements 


Climatic  Influence  on  topographic  form  proper¬ 
ties  of  f  1  u v  tally  dissected  landscapes  has  long  been 
recognized  by  geologists  and  geographers.  Al¬ 
though  considerable  literature  exists  on  this  sub¬ 
ject.  three  vital  procedures  are  lacking  in  early 
studies:  (1)  Isolation  of  climatic  elements  and  ef¬ 
fects.  Climate  was  usually  referred  to  as  “arid" 
or  “humid."  in  regard  to  moisture:  “hot."  “tem¬ 
perate,"  or  “cold,"  In  regard  to  temperature.  (2) 
Investigation  of  relative  importance  of  processes. 
Knowledge  of  process  was  largely  taken  for  granted 
emphasis  being  placed  on  historical  interpretation 
and  classification.  (3)  Expression  of  laws  in  quan¬ 
titative  terms.  Research  was  largely  of  a  qualita¬ 
tive  nature  because  it  was  thought  that  features  of 
the  landscape  are  so  variable  that  mathematical 
treatment  would  be  impracticable.  Attempts  at 
quantitative  application  of  principles  were  resticted 
to  simple  cases. 

Recently,  however,  a  number  of  geomorpholo¬ 
gists  have  recognized  the  availability  of  mathema¬ 
tical  tools  capable  of  coping  with  the  variability, 
or  chance  element,  in  landforms.  Their  approach 
has  been  chiefly  to  examine  in  detail  the  variations 
in  properties  of  areas  similar  in  major  character¬ 
istics  and  history.  If  it  can  be  shown  that  these 
variations  are  too  great  to  be  of  random  nature, 
then  they  must  be  the  result  of  some  basic  differ¬ 
ences  in  rate  or  intensity  of  operation  of  process. 

In  these  recent  studies,  climate  has  been  for  the 
most  part  of  secondary  interest. 

The  main  purpose  of  this  study  is  to  provide  a 
foundation  for  further  research  into  the  effects  of 
availability  of  moisture  to  vegetal  cover,  and  of 
rainstorm  intensity  and  frequency  on  processes  of 
landform  evolution.  Consideration  is  given  to  pro¬ 
perties  of  the  surface  mantle  as  they  are  affected 
by  climate,  and  in  turn  as  they  affect  processes. 
Because  action  of  the  processes  is  too  slow  to  be 
observed  during  a  short  time,  relative  rates  are 
inferred  from  measurements  of  the  landforms 
themselves. 

The  factors  influencing  topographic  texture  and 
average  slope  can  be  considered  as  fundamental  in 
control  of  the  fluvial  landscape.  The  controls  most 
often  mentioned  in  descriptive  treatments  of  the 
subject  are  “climate,”  “rock  permeability,”  and 
“rock  strength.”  Structural  properties  of  the  rock 
such  as  jointing,  initial  or  available  relief,  stage 
in  the  erosional  cycle,  and  condition  of  the  surface 
soil  are  frequently  mentioned  as  well  (Thornbury, 


1954.  p.  128).  Factors  influent  mg  the  steepness 
tn  ’  mode  of  retreat  of  slopes  in  fluv tally- controlled 
landscapes  have  received  much  attention  in  testing 
Walther  Penck's  ideas  of  slope  development. 

Gilbert  (1877,  p.  109-118)  analyzed  the  processes 
of  slope  and  channel  equilibrium  as  effected  by  run¬ 
ning  water.  He  reasoned  that  flu vially  eroded  areas 
in  arid  regions  are  “degraded"  at  a  rate  limited  by 
the  slowness  of  the  production  of  fine  waste  mate¬ 
rial;  removal  by  transport  is  rapid  once  waste  has 
formed.  Thus,  differences  in  rock  texture  are  accen¬ 
tuated  in  arid  regions.  Greater  availability  of  water 
hastens  the  rate  of  production  of  waste  material, 
and  likewise  increases  the  rapidity  of  transport. 

But  the  greater  vegetal  cover  in  humid  regions 
slows  transportation  and  is  generally  effective  in 
causing  an  accumulation  of  waste.  Thus,  differences 
in  rock  texture  are  obscured  by  accumulated  waste 
in  humid  regions.  A  climate  intermediate  between 
arid  and  humid.  Gilbert  concluded,  must  produce 
the  most  rapid  and  effective  erosion. 

In  badlands  areas,  declivities  decrease  as  the 
amount  of  water  flowing  over  them  increases,  and 
slopes  are  greater  “in  proportion  as  they  are  near 
divides,  unless  they  are  very  near  divides"  (Gil¬ 
bert.  1877,  p.  117).  The  convexity  of  hilltops  evi¬ 
dently  puzzled  Gilbert  at  the  time  he  wrote  the  re¬ 
port  on  the  Henry  Mountains,  and  his  solution  to  the 
problem,  as  indicated  later,  was  the  agency  of  creep 
(Gilbert,  1909,  p.  344-350).  Little  has  been  added  to 
the  qualitative  knowledge  of  slope  equilibrium  since 
Gilbert’s  publications;  quantitative  study  of  these 
processes  described  has  only  just  begun. 

Walther  Penck  (1953.  p.  119-120.  and  130-133) 
evidently  believed  that  the  same  landforms  can  ap¬ 
pear  under  all  climates,  though  developing  at  dif¬ 
ferent  rates.  Climatic  differences  may  change  the 
relationship  of  mass-wasting  to  transportation,  but 
they  do  not  produce  any  fundamental  change  in  the 
denudational  forms.  Slope  gradients  depend  only  on 
the  intensity  of  erosion;  so  if  all  slopes  in  a  region 
are  of  the  same  steepness,  it  is  an  expression  of  the 
fact  that  erosion  has  acted  upon  them  with  the  same 
intensity.  It  follows  from  this,  according  to  Penck, 
that  gradients  are  determined  by  endogenetic  fac¬ 
tors,  and  less  so  by  rock  hardness.  The  form  of 
slope  profiles  reflects  changes  in  endogenetic  rates. 
Old  highlands,  separated  from  a  former  base-level 
of  erosion  by  a  topographic  break,  invariably  dis¬ 
play  concave  slopes.  Convex  and  straight  slopes  are 
associated  with  zones  of  intense  recent  erosion. 
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When  tin  niiMii  slope  gradient  corresponds  to  the 
maximum  gradient  for  the  rock  concerned,  coin  ex 
slopes  are  replaced  by  straight  profiles;  these  find 
Uietr  “purest  expression  in  the  steep  relief  of 
mountains  which  escaped  Pleistocene  glaciation, 
whether  these  are  covered  with  tropical  virgin  for¬ 
est  or  belong  to  arid  regions’  (Penck,  1953,  p.  132). 
Penck's  views  of  relation  of  slope  form  to  endogene- 
tu  factors  have  largely  been  discredited  in  univer¬ 
sal  application.  It  is  possible  that  qualifications 
may  have  originally  been  present  in  his  book  that 
have  been  lost  in  translation,  so  great  is  the  mass 
of  verbiage. 

The  first  truly  quantitative  approach  to  the  sub¬ 
ject  of  the  control  of  topographic  texture  was  by 
hydrologists.  Robert  E.  Horton  (1945,  p.  283,  285) 
defined  the  terms  drainage  density,  the  total  length 
of  channels  divided  by  the  area  drained,  and  stream 
frequency,  the  number  of  channel  segments  divided 
by  the  area  drained,  showing  that  together  they  de¬ 
termine  the  fineness  of  topography.  He  concluded, 
perhaps  erroneously,  that  areas  with  low  drainage 
density  are  poorly  drained,  and  areas  with  high 
drainage  density  are  well  drained.  Horton  em¬ 
phasized  the  importance  of  surface  characteristics 
in  determining  drainage  density,  particularly  in¬ 
filtration  capacity,  and  resistance  of  the  surface 
to  erosion  (1933,  p.  260,  and  1945,  p.  284). 

Forest  and  range  conservationists,  working  un¬ 
der  the  U.S.  Department  of  Agriculture,  have  ob¬ 
tained  much-njeded  quantitative  information  on 
infiltration  and  erosic.i  rates  related  to  land  use, 
type  of  soil  and  cover,  and  distribution  of  rain  dur¬ 
ing  storms  (Dortignac,  1951;  Love  and  Johnson, 

1952;  Love  and  Dunford,  1952;  Love,  1953;  Anderson, 
1954).  Of  special  importance  to  geomorphologists 
are  studies  of  rainfall- runoff  proportions,  in  cases 
of  known  accelerated  erosion,  in  both  forested  and 
grazed  watersheds.  Schumm  (1955)  has  demon¬ 
strated  that  sediment  loss  from  basins  can  be  re¬ 
lated  to  form  properties  within  the  basin,  and  to 
infiltration  capacity  of  the  soil. 

A.  N.  Strahler  has  emphasized  the  use  of  statis¬ 
tical  methods  in  quantitative  investigations,  rather 
than  deductive  analysis  (1950,  1954a).  Recently  he 
has  used  the  Pi  Theorem  and  dimensional  analysis 
in  combining  factors  that  rationally  can  be  said  to 
contribute  to  the  main  causes  determining  drainage 
density  (Strahler,  1956b).  Dimensionless  numbers 
wore  derived  that  have  the  property  of  independence 
of  size  of  corresponding  units  of  the  drainage  basins 
considered.  Thus,  if  drainage  density  can  rationally 
be  said  to  depend  upon  (1)  the  average  rate  of  flow 


stated  in  volume  per  unit  of  area  of  basin,  (2)  an 
erosion  proportionality  factor,  (3)  the  height  of  the 
basin  above  its  mouth,  (4)  density  of  surface  ma¬ 
terial.  (5)  viscosity  of  water,  and  (6)  the  accelera¬ 
tion  of  gravity,  then  these  quantities  can  be  arranged 
into  four  dimensionless  numbers.  These  are  (1)  rug¬ 
gedness  number  (product  of  drainage  density  and 
relief  of  basin).  (2)  Froude  number,  (3)  Reynolds 
number,  and  (4)  Horton  number  (product  of  erosion 
proportionality  factor  and  the  average  volume  rate 
of  flow  per  unit  area).  Full  evaluation  of  the  rela¬ 
tions  among  all  these  quantities  would  take  many 
years  of  research.  The  present  study  treats  only 
portions  of  these  relations  that  are  inherent  In  ex¬ 
pression  of  relief,  strength  of  surface  material, 
average  rate  of  runoff,  measures  of  roughness,  and 
extent  of  < 

Generalizing  from  these  discussions,  there 
seems  to  be  agreement  on  factors  influencing  slope 
development  and  equilibrium,  and  drainage  density. 
Considerable  disagreement  exists  as  to  which  fac¬ 
tors  are  dominant  in  effect  and  which  are  secondary. 
Where  topographic  features  and  related  quantities 
are  measurable,  it  is  generally  not  known  whether 
these  quantities  are  related  by  simple  proportion, 
by  power  functions,  by  exponential  functions,  or  by 
other  transcendental  functions.  The  writer  has  en¬ 
deavored  to  advance  the  understanding  of  these  fun¬ 
damental  questions  by  using  (1)  form  measurements 
from  topographic  maps,  (2)  simple  field  measure¬ 
ments  of  soil  properties  and  characteristics  of 
vegetal  cover,  (3)  analysis  of  climate  elements,  and 
(4)  combination  of  numerical  data  by  statistical 
analysis.  Correlation  analysis,  in  particular,  allows 
separation  of  major,  minor,  and  definitely  unre¬ 
lated  effects. 

This  investigation  formed  part  of  a  quantitative 
study  of  erosional  landforms  sponsored  through  Co¬ 
lumbia  University  by  the  Geography  Branch  of  the 
Office  of  Naval  Research  as  Project  Number  NR 
389-042,  under  Contract  N60NR  271,  Task  Order  30. 

The  writer  is  indebted  to  Professor  A.  N. 

Strahler  for  guidance  in  planning  and  carrying  out 
field  and  map  regard).  Professor  Howard  Levene 
of  the  Department  of  Mathematical  Statistics  of 
Columbia  University  advised  on  a  method  of  analy¬ 
sis  of  variance  and  critically  read  the  text.  Mr. 

Gary  Gordon  of  Mt.  V  ;rnon,  Iowa,  assisted  in  the 
field  during  two  months  of  the  summer  of  1955. 

Sara  Melton  assisted  during  two  field  seasons,  and 
helped  in  preparing  the  manuscript.  Styron  Douthit 
and  Steven  Haines  uf  the  University  of  Oklahoma 
aided  in  preparation  of  the  illustrations. 


Methods  of  Map  Analysis  and  Definitions  of  Terms 


Geologists  familiar  with  quantitative  landform 
analysts  will  understand  most  of  the  terms  used 
here,  although  some  are  original.  To  avoid  misun¬ 
derstanding,  all  terms  and  symbols  are  defined  in 
Table  1,  which  gives  units  and  dimensions  of  all 
quantities.  Map  dr  ta  is  included  from  3  fifth-order 
basins,  22  fourth- order  basins,  and  <3  third-order 
basins.  Field  me  surements  are  given  for  1  fifth- 
order  basin.  19  tourth-order  basins,  and  9  third- 
order  basins.  These  are  divided  into  58  third-order 
basins  for  analysis.  The  remaining  basins  were  in¬ 
spected  briefly  to  determine  the  general  condition 
of  the  vegetal  cover,  rock  type,  and  condition  of 
channels. 


MAP  ACCURACY 

The  amount  of  usable,  reliable  information  con¬ 
tained  on  topographic  maps  increases  with  the  ac¬ 
curacy  of  the  map.  Because  maps  differ  in  accuracy 
depending  on  scale  and  contour  interval,  and  oi.  .he 
method  by  which  they  were  made,  it  is  important  in 
comparative  studies  that  data  be  obtained  from 
maps  of  comparable  relative  accuracy.  Otherwise, 
comparisons  between  areas  differing  in  some  re¬ 
spect  will  be  in  error. 

Only  those  maps  published  by  the  U.  S.  Geologi¬ 
cal  Survey  at  the  scale  1:24,000,  and  compiled 
from  aerial  photographs,  have  been  used  in  this 
study,  with  three  exceptions:  (1)  maps  of  the  Neder¬ 
land,  Black  Hawk,  and  Morrison  quadrangles  in 
Colorado,  were  published  at  the  scale  1:24,000,  but 
were  made  by  plane-table  sux  veys  in  1942  and  1938, 
before  photogrammetric  methods  were  in  use  by 
the  U.  S.  Geological  Survey.  The  drainage  nets  of 
the  basins  studied  in  these  areas  were  drawn  from 
aerial  photographs  obtained  from  the  U.S.  Forest 
Service,  by  using  a  KEK  Stereoscopic  Plotter.* 
Considerable  revision  in  the  drainage  nets  of  basins 
in  those  areas  was  necessary.  (2)  A  small  badlands 
area  in  the  Chinle  shale,  north  of  Cameron,  Ari¬ 
zona.  is  not  included  in  any  quadrangle  map,  and 
could  not  be  adequately  shown  at  the  scale  1:24,000. 
Two  small  drainage  basins  in  the  badlands  were 
mapped  at  1  inch  to  50  feet  by  E.  D.  Koons,  S.  A. 
Eftitmiiii,  and  A.  N.  Strahkf .  In  the  summer  J 
1953.  When  examined  by  the  writer,  two  years  later, 
two  small  channels  had  become  obliterated  and 

♦Manufactured  by  Kail  and  Associates,  Denver,  Colorado. 


were  removed  from  the  drainage  net;  otherwise, 
the  map  was  extremely  accurate.  (3)  An  area  of 
badlands  in  the  Saguaro  National  Monument,  south¬ 
east  of  Tucson,  Arizona,  is  covered  only  by  the 
1904  Tucson  Quadrangle,  at  the  scale  1:125,000.  In 
order  to  use  the  area  in  the  study,  a  fifth-order 
basin  was  mapped  at  1  inch  to  100  feet  with  plane 
table  and  alidade  by  the  writer  and  assistant  in 
December,  1955.  '  he  drainage  nets  mapped  in 
these  three  areas  are  believed  comparable  in  ac¬ 
curacy  to  the  drainage  nets  shown  on  photogram- 
n.etric  maps  of  other  areas. 

In  all  basins  in  which  field  measurements  were 
taken,  the  drainage  net  was  drawn  on  the  map  to 
correspond  with  the  existing  channels.*  It  was 
found  that  in  most  maps  a  high  percentage  of  the 
first-order  channels  is  represented  by  the  small¬ 
est  cuspo  in  contour  lines.  In  a  few  areas  (e.g. 

Spring  Creek  basin,  Springville  Quadrangle,  Utah), 
the  smallest  cusps  represent  only  irregularities  in 
the  surface,  not  persistent  drainage  lines.  In  other 
areas  'n.g.  Quartz  Hill  basin,  Paxton  Springs 
Quadrangle,  N.  M.),  the  smallest  cusps  in  the  con¬ 
tour  lines  represent  second  or  third-order  channels. 
In  many  basins  shallow  gullies  established  on  the 
valley  sides  are  not  shown  by  the  contour  lines.  In 
accordance  with  these  considerations,  about  95% 
of  the  basins  studied  needed  at  least  minor  correc¬ 
tions  of  the  drainage  net  as  shown  on  the  contour 
map.  (See  Appendix  IV  for  photo-copies  of  portions 
of  maps,  showing  drainage  nets  of  all  basins  studied.) 

METHODS  OF  OBTAINING  QUANTITATIVE  DATA 

FROM  MAPS 

Many  of  the  basins  studied  are  rather  small,  as 
shown  on  the  original  maps.  In  measuring  features 

♦Horton  (1945)  used  stream  with  the  same  meaning  that 
channel  has  here  and  in  previous  reports  by  A.  N. 

Strahler  (1950,  1952a,  1953),  V.  C.  Miller  (1953),  K.  G. 
Smith  (1953),  S.  A.  Schumm  (1954)  and  D.  R.  Coates 
(1955).  The  term  channel  means  a  permanent,  clearly 
defined  trench  or  trough  clearly  showing  evidence  of 
scour  by  channel  flow  end  hounded  hy  valley  stiles  slop¬ 
ing  toward  the  channel  axis.  Excluded  by  this  definition 
are  ephemeral  or  seasonal  rills  cut  in  o'.,  erwiso  smooth 
u.hliJh  il  inf  nmlriblili'  iriroti  Ut  L1  >'  rill  8 
which  have  no  clear'v  defined  water  .-.hi..1,  and  drainage 
divide.  Channel  is  preferred  since  emphasis  is  placed 
on  the  depression  in  which  water  flows,  rather  than  a 
stream  of  water  itself.  Most  channels  studied  here  would 
be  classed  as  ephemeral  or  intermittent  streams. 
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1m  thrtff  Ixuiinn,  the  per  (‘enlace  of  error  can  tx1  re¬ 
duced  by  enlarging  tli»*  inapn  photographically; 
nuMxuroim’nl.H  on  the  prints  can  then  lx-  made  with 
greater  relative  preci.snm  The  maps  of  all  lias  ms 
studied  in  letatl  In  the  field  were  diluted  two  to 
eight  limes  depending  upon  original  size.  Prints 
were  made  by  a  Photostat  machine,  or  by  making 
full-st/.ed  negatives  on  la  rite  sheets  of  film  and 
printing  directly  on  blue-line  paper.  The  prints 
were  air  dried  Distortion  of  the  linage  was  presen’ 
In  almost  all  cases  but  by  marking  several  one- 
inch  Intervals  with  various  orientations  on  the  or¬ 
iginal  and  measuring  the  corresponding  interval  on 
the  enlarged  print,  it  was  established  that  differ¬ 
ential  shrinkage  caused  slight  variation  in  scale. 

In  two  perpendicular  directions,  parallel  to  the 
edges  of  the  print,  the  greatest  difference  In  scale 
was  1.6%. 

Measurements  of  length  weie  made  from  the 
enlarged  maps  by  means  of  a  Dietzgen  Map  Meas¬ 
ure,  No.  1719B,  scaled  to  1/32  inch,  and  on  which 
measurements  can  be  estimated  to  the  nearest  1/64 
Inch.  Hence,  on  an  enlarged  map  with  scale  of  about 
10  inches  per  mile,  measured  lengths  can  justifi¬ 
ably  be  recorded  to  the  nearest  0.01  mile. 

Area  measurements  were  made  with  a  Keuffel 
and  Esser  Compensating  Polar  Planimeter,  Model 
4242.  Highly  consistent  results  were  obtained  in 
repeated  measurements  of  the  same  outline.  Re¬ 
versing  the  position  of  the  pole  with  respect  to  the 
tracer  head  eliminated  a  small  systematic  error. 
Deviations  greater  than  0.02  square  inches  are  un¬ 
usual  in  individual  area  measurements,  hence,  on 
a  map  with  scale  100  square  inches  per  square 
mile,  it  is  justifiable  to  record  measurements  to 
the  nearest  0.001  square  mile. 


PROPERTIES  MEASURED  FROM  MAPS 

Of  all  the  information  that  can  be  obtained  from 
maps  by  measurement,  two  general  categories  may 
be  distinguished:  (a)  quantities  measured  (or 
counted)  directly,  with  no  more  computation  than 
multiplication  by  a  scale  factor,  and  (b)  quantities 
that  must  be  computed  from  direct  measurements 
because  they  cannot  in  themselves  be  measured. 
The  first  category  includes  length  and  area  meas¬ 
urements.  numbers  of  items  obtained  by  enumera¬ 
tion,  and  certain  angular  measurements.  Lower¬ 
case  Roman  letters  are  used  as  notation  for 
properties  of  length,  area,  or  number;  lower-case 
Greek  letters  for  angular  measurements.  To  de¬ 
velop  a  satisfactory  notation  for  the  numerous 
morphometric  properties,  full  use  of  upper  and 
lower-case  letters,  and  of  subscripts,  is  necessary. 
Subscript  notation  must  be  systematic  and  consist¬ 
ent.  Subscripts  are  appended  to  letter  symbols  to 
distinguish  measurements  of  the  same  kind  that 


are  iii.ulr  or.  features  differing  In  some  way,  such 
an  dj.  mmtlx-r  of  first-order  channels,  n-  mini- 
tx-r  of  second-order  channels,  and  so  forth  The 
notation  inherited  by  geomorphologisls  and  now 
employed  by  hydrologists  is  remarkably  poor  In 
the  use  of  subscripts.  The  principal  objection  Is 
Hurt  the  same  symbol,  (or  Instance  F.  Is  used  for 
such  different  quantities  as  force  and  frequency, 
with  the  difference  Indicated  by  a  subscript.  The 
writer  has  followed  a  conservative  practice  In  us¬ 
ing  subscripts  only  to  denote  the  order  of  the 
feature  upon  which  the  measurement  is  taken, 
such  as  nj,  number  of  first-order  channels,  and 
I)  .  total  length  of  first-order  channels. 

Channel  order  states  the  relative  magnitude  of 
a  segment  of  stream  channel  in  a  hierarchy  de¬ 
termined  by  arrangement  of  tributaries  with  re¬ 
spect  to  the  main  trunk.  The  system  of  ordering 
channel  segments  proposed  by  Horton  (1945)  and 
modified  by  Strahler  (1952b,  footnote,  p.  1120)  is 
followed  here.  The  smallest  permanent  channels, 
those  having  no  tributaries,  are  designated  first- 
order  channels.  At  the  junction  of  two  first-order 
channels,  a  second-order  channel  begins;  at  the 
junction  of  two  second-order  channels,  a  third- 
order  channel  begins,  and  so  forth.  A  channel  of 
any  given  order  can  directly  receive  tributaries 
of  all  lesser  orders.  The  trunk  stream,  or  master 
stream,  through  which  all  discharge  passes  has 
the  highest  order.  The  area  drained  by  a  c’  ?1 

of  order  u,  down  to  its  junction  with  anoth  in- 

nel  of  order  equal  to  or  greater  than  u,  is  g- 
nated  a  basin  of  order  u. 

Because  of  the  way  in  which  channel  ordering 
is  defined,  it  is  possible  to  state  the  minimum 
possible  number  of  channels  of  each  order  in  a 
basin  of  a  given  order.  Thus,  for  a  second-order 
channel,  at  least  two  first-order  channels  must 
be  present;  for  a  third-order  channel,  at  least 
two  second-order  channels  and  four  first-order 
channels  are  needed,  for  a  fourth-order  channel, 
two  third-order  channels,  four  second-order  chan¬ 
nels,  and  eight  first-order  channels  are  needed. 

In  general,  in  a  basin  of  order  uD  the  minimum 
number  of  channel  segments  of  order  u,  is 

n„  =  2  u“'u  (1) 

The  minimum  total  number  of  segments  of  all 
orders  is 

“u  Uu 

(2) 

ii  =  l  u  *  1 

Basin  area,  or  planimetric  area,  a,  is  the  area 
in  square  miles  of  the  outline  of  the  watershed  of 
a  channel  as  projected  onto  the  horizontal  map 
plane.  This  area  is  commonly  used  to  calculate 
the  total  amount  of  water  entering  the  system 
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Figure  1.  Inflection  angle  as  measured  on  various  forms  of  first-order  channels. 


during  a  rain.  Basin  area  is  a  fundamental  dimen¬ 
sional  property  in  the  majority  of  the  quantitative 
expressions  of  drainage-basin  characteristics. 

Basin  perimeter  is  the  length  along  the  major 
divide  as  projected  onto  the  horizontal  plane  of  the 
map.  The  symbol  used  here  for  the  perimeter  is  p; 
the  unit  of  measure  is  miles.  The  perimeter  is  a 
linear  measure  of  the  size  of  the  basin,  and  is 
largely  independent  of  the  texture  of  the  topography. 
It  is  used  in  conjunction  with  basin  area  to  give  a 
measure  of  the  departure  of  the  basin  from  a  true 
circle  (see  basin  circularity),  and  in  conjunction 
with  the  relief  to  give  a  measure  of  general  steep¬ 
ness  of  the  basin  (see  relative  relief). 

Basin  is  the  vertical  distance  between  the 

mouth  of  i-.it;  basm  and  tht  highest  point  in  the  ba¬ 
sin,  which  in  the  basins  m  tills  study  is  always 
located  on  the  major  divide.  The  symbol  for  relief 
is  r;  the  unit  of  meaourement  is  feet.  In  order  to 
find  the  relief,  it  is  usually  necessary  to  inter¬ 
polate  between  contour  lines  at  the  mouth  and  high¬ 
est  point.  This  is  done  to  the  nearest  one-fourth 
contour  interval,  since  it  is  desirable  to  avoid  the 
rigorous  but  clumsy  expression  of  r  given  with  a 
possible  range  of  two  contour  intervals. 

Total  length  of  channel  segments  of  each  order 
within  the  basin  were  measured  and  recorded  un¬ 
corrected  for  gradient,  in  units  of  miles.  The  sym¬ 
bol  used  for  the  total  length  of  channel  segments  or 


order  u  is  lu.*  The  mean  length  of  channel  seg¬ 
ments  of  order  u  is  ill,  where  n„  is  the  number  of 

nu 

segments  of  order  u,  and  is  represented  by  the 
symbol  1„.  The  total  length  of  channels  of  all  or¬ 
ders  in  a  basin  is  I  lu. 

Inflection  angle  of  contour  lines  is  the  angle 
which  a  contour  line  makes  with  itself  where  it  de¬ 
picts  a  channel  and  is  given  the  symbol  The  in¬ 
flection  angle  has  apparently  not  previously  been 
studied.  In  this  study,  the  minimum  inflection  angle 
of  each  first-order  stream  was  measured  from  a 
map  by  use  of  protractor  and  straightedge.  Figure 
1  jhuws  several  common  forms  of  first-order  chan¬ 
nels,  as  delineated  by  contour  lines,  and  the  inflec¬ 
tion  "ngle  for  each  is  illustrated. 

It  is  obvious  that  geometrically  the  gradient  of 
a  channel  can  vary  from  zero  up  to  the  angle  of 
slope  of  the  adjacent  valley  sides,  but  no  steeper. 
Increasing  difference  between  the  channel  gradient 
and  the  slope  of  the  adjacent  valley  sides  produces 
a  more  arute  inflection  angle.  Thus,  the  inflection 

♦Confusion  may  arise  because  of  the  identity  of  form  of 
the  ietter  1  and  numeral  1.  In  this  notation  all  subscripts 
refer  to  order  number,  hence  a  subscript  1  always  means 
the  numeral  one.  The  letter  1  is  always  followed  by  a 
subscript  denoting  channel  order  Therefore,  1  appearing 
without  a  subscript  always  means  one. 


mi  more  Ilian  a  general  indication  of  the  Inflection 
anglcri  can  lx-  obtained  from  map  measurement  uf 
ftriil-order  rlunneln. 
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i titles  of  first-order  ihannels  are  usually  greater 
than  of  higher  order  channels.  In  Figure  1  a.  the 
dunnel  gradient  is  nearly  the  name  an  the  slope  of 
the  \alley  Hide.  and  elves  a  broadly  obtuse  inflec- 
"  »  angle.  In  Figure  l  b.  the  gradient  in  some¬ 
what  less  than  the  slope  of  the  adjacent  valley 
sides  and  the  inflection  angle  is  acute,  in  Figure 
1  i  the  channel  gradient  in  much  loss  than  the 
adjacent  valley-side  slopes,  and  the  Inflection  an¬ 
gle  is  even  more  acute.  As  valleys  .ire  deepened 
the  angle  of  inflection,  as  shown  on  a  map  of  the 
channel,  may  be  expected  to  become  more  acute. 
Thus,  the  evolution  of  a  first-order  tin  sin.  Initiated 
by  a  small  gullv  ha\  Ing  as  Its  gradient  the  slope  of 
the  surface  upon  which  it  formed,  to  a  mature 
basin  form  In  which  channel  gradient  Is  consider¬ 
ably  lower  than  the  slope  of  the  adjacent  valley 
sides,  might  be  traced  by  changes  in  the  Inflection 
angle. 

New  tributaries  to  existing  channels  tend  to 
develop  perpendicular  to  the  contour  line,  so  as  to 
have  the  steepest  possible  slope.  The  projection  of 
the  new  tributary  on  a  map  makes  an  angle  with 
the  main  channel  equal  to  the  complement  of  the 
angle  which  the  contour  line  makes  with  the  main 
channel.  Main  channels  which  have  highly  acute 
inflection  angles  will  tend  to  develop  tributaries 
entering  at  nearly  right  angles.  The  relation  be¬ 
tween  angle  of  entry  (axil  angle*),  valley-side 
slope,  and  channel  gradient  was  found  by  Horton  to 
be 


cos£  = 


tan  y 
tan  0 


where  £  is  the  axil  angle,  y  the  channel  gradient, 
and  G  the  valley-side  slope  (Horton,  1945,  p.  349). 
This  equation  can  be  derived  from  equation  19  in 
Appendix  I,  by  setting 
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90  -  ^ 


The  relations  among  channel  gradient,  valley- 
side  slope,  dihedral  angle  between  valley-sides, 
and  the  inflection  angle  can  be  expressed  in  geo¬ 
metric  laws,  and  the  derivation  of  these  laws  is 
given  in  Appendix  I. 

The  measurements  of  inflection  angles  of  first- 
order  channels  taken  in  this  study  are  biased  by 
the  habit  or  convention  of  indicating  small  streams 
adopted  by  the  individual  cartographer.  It  is  there¬ 
fore  advisable  in  future  studies  to  take  inflection- 
angle  measurements  on  small  channels  in  the  Held. 
Reliable  results  from  map  measurements  mignt  be 
obtained  for  second  or  third-order  channels,  but 

♦Axil  angle  has  been  used  by  botanists  for  150  years  to 
denote  the  upper  angle  between  the  main  stem  anci  a 
branch.  Its  use  in  this  sense  is  analogous. 


DFH1VF1)  QUANTITIFS 

Derived  quantities  are  those  that  must  lx*  ob¬ 
tained  by  calculation  from  the  properties  measured 
or  cotinted  directly  from  maps.  The  accuracy  of  the 
derived  oment'es  depends  directly  on  the  accuracy 
of  measured  quantities,  and  their  frequency  dis¬ 
tributions  depend  on  the  distributions  of  all  the 
measured  Items  involved  In  their  calculation.  A 
greater  amount  of  freedom  Is  Involved  tn  their  es¬ 
timation.  For  these  reasons  It  is  desirable  to  dis¬ 
tinguish  between  measured  and  derived  quantities. 
Upper -case  Roman  letters  are  used  here  as  sym¬ 
bols  for  derived  quantities,  except  that  derived 
angular  quantities  are  represented  by  lower-case 
Greek  letters. 

Drainage  density,  D.  Is  the  ratio  of  the  total 
length  of  channels  of  all  orders  in  a  basin  to  the 
v  i 

area  of  the  basin,  — ,  and  is  therefore  the  num- 

ber  of  miles  of  channel  maintained  by  a  square 
mile  of  drainage  area  (Horton,  1945,  p.  283).  Drain¬ 
age  density  is  a  numerical  expression  of  the  fine¬ 
ness  of  texture  of  the  basin,  but  because  a  particu¬ 
lar  value  is  usually  nearly  constant  over  an  area  of 
homogeneous  lithology,  climate,  and  stage  of  de- 
\elopment,  it  is  independent  of  the  basin  size. 
Drainage  density  has  the  dimension  of  inverse  of 
length,  and  is  a  useful  indicator  of  scale  of  size 
of  planimetric  units  in  a  fluvially  controlled  land¬ 
scape. 

Density  of  first-order  channel  segments.  Dj,  is 
the  total  length  of  first-order  channel  segments  in 
a  basin  of  some  larger  order,  divided  by  the  area 

of  the  entire  basin,  ----- .  It  is  useful  in  determin- 
a 

ing  whether  a  drainage  density  higher  in  one  area 
than  in  another  is  due  to  a  greater  length  from 
headward  growth  of  the  first-order  channels. 

Channel  frequency,  F,  is  the  ratio  of  the  total 
number  of  channel  segments  of  all  orders  in  a 

basin  to  the  area  of  the  basin,  (Horton,  1945, 

cl 

p.  285).  T*  is  a  supplementary  measure  of  the  fine¬ 
ness  of  texture  of  the  topography,  but  may  be  in¬ 
dependent  of  drainage  density  because  it  depends 
only  on  number  of  channel  segments,  not  on  the 
lengths. 

Frequency  of  first-order  channels,  Fi,  is  the 
total  number  of  first-order  channels  divided  by  the 

area  of  the  entire  basin,  This  measure  is 

a 

useful  in  determining  whether  a  drainage  density 
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higher  in  one  area  than  In  another  has  resulted 
from  growth  of  new  channels,  or  simply  a  lengthen¬ 
ing  of  existing  channels. 

Ituggedness  number,  H,  is  the  product  of  drain¬ 
age  density  and  rcTTcFdf  basin,  lx)th  in  the  same 

un“-s:  h  =  sM)  ur  H  =  'sm r  (s,rahlcr- 1954c) 

It  is  a  dimensionless  number,  since  D  has  the  di¬ 
mension  [l  1  ],  and  r  has  the  dimension  (L|.  Areas 
of  low  absolute  relief,  but  high  drainage  density 
are  as  rugged,  when  their  small  size  is  considered, 
as  areas  with  great  relief  but  low  drainage  density. 
Areas  that  l  re  geometrically  similar  in  planime- 
trk  and  relief  features  will  have  the  same  value 
of  the  ruggedness  number. 

Relative  relief,  R,  which  is  a  measure  of  the 
general  steepness  of  the  basin  from  summit  to 
mouth,  can  be  measured  in  a  variety  of  ways. 
Schumm  (1954,  p.  22)  defined  relief  ratio  as  "the 
ratio  between  the  total  relief  of  a  basin.  .  .and 
the  longest  dimension  of  the  basin  parallel  to  the 
principal  drainage  line.  This  relief  ratio  is  a  di¬ 
mensionless  height-length  ratio  and  is  equal  to  the 
tangent  of  the  angle  formed  by  two  planes  inter¬ 
secting  at  the  mouth  of  the  basin;  one  representing 
the  horizontal,  the  other  passing  through  the  high¬ 
est  point  of  the  basin.”  In  oddly  shaped  basins,  or 
basins  with  the  highest  point  at  some  place  other 
than  at  the  position  farthest  removed  from  the 
mouth,  there  is  a  question  as  to  what  should  be 
taken  as  the  length  measurement,  and  the  interpre¬ 
tation  of  the  relief  ratio  as  the  tangent  of  the  angle 
between  two  planes  is  not  appropriate. 

For  this  study  relative  relief.  R,  is  defined  as 
the  ratio  of  the  basin  relief,  expressed  in  units  of 

miles,  to  the  length  of  perimeter:  R  =  The 

advantage  of  this  measure  over  relief  ratio  is  that 
both  relief  and  perimeter  had  been  measured  for 
other  purposes;  the  position  of  the  highest  point  of 
the  basin  with  respect  to  the  mouth  does  not  affect 
the  measure.  A  rational  relation  between  the  meas¬ 
ure  of  relief  used  here  and  Schumm’s  relief  ratio 
can  be  developed  by  assuming  various  shapes  for 
the  basin. 

Ratio  of  channel  lengths  to  basin  perimeter,  P, 
is  a  dimensionless  measure  of  fineness  of  topo¬ 
graphic  texture  relative  to  basin  size.  Two  basins 
of  different  sizes,  but  with  geometrically  similar 
planimetric  properties,  will  have  different  values 
of  drainage  density,  but  identical  values  of  the 
channel-perimeter  ratio. 

Basin  circularity,  C.  measures  the  degree  to 
which  a  basin  approximates  a  circle  in  planimetric 
outline.  Circularity  is  defined  as  the  area  of  the 
basin  divided  by  the  area  of  the  circle  with  equal 
perimeter  (Miller,  1953,  p.  8).  The  circle  can  be 


defined  uniquely  as  that  plane  figure  in  which  the 
greatest  area  is  bound  by  the  shortest  possible 
perimeter.  Therefore,  if  the  area  of  any  irregular 
figure,  such  as  a  drainage-basin  outline,  be  divided 
by  the  area  of  the  circle  having  the  same  perimeter 
the  quotient  will  be  a  rumber  between  zero  and  one 
The  measure  is  ambiguous,  though,  as  it  does  not 
discriminate  between  departures  from  perfect  cir¬ 
cularity  due  to  irregularities  In  gross  outline,  or 
to  crenulations  in  the  perimeter. 

Bifurcation  ratio,  BU:u*l  =  — — .  is  the  ratio 

- nu  ♦  i 

betw-een  the  number  of  channel  segments  of  a  given 
order  and  the  number  of  channel  segments  of  the 
next  higher  order.  According  to  Horton’s  law  o. 
stream  numbers,  values  of  B  tend  to  be  the  same 
for  all  values  of  u  within  a  given  basin  (Horton,  1945, 
p.  286).  The  value  of  B  characteristic  of  all  orders 
within  the  basin  is  the  slope  of  the  regression  line 
when  numbers  of  channel  segments  are  plotted  on 
semi-log  paper  against  order  number  (Strahler. 
1954b,  p.  345). 

Length  ratio,  L  „  - „ . i  =  — _ •  is  the  ratio  of  the 

Tu-1 

mean  lengths  of  channels  of  a  given  order,  to  the 
mean  length  of  channels  of  the  next  lower  order. 
According  to  Horton’s  law  of  stream  lengths,  the 
mean  channel  length  for  each  order  increases  as 
an  exponential  function  of  channel  order  (Horton, 

1945,  p.  236).  Whether  the  values  of  Lu:u_i  for  a 
basin  are  constant  within  a  basin,  or  change  sys¬ 
tematically  with  order  is  a  question  that  must  be 
answered  by  statistical  analysis. 

Maximum  valley-side  slope.  O.  is  the  slope  in 
degrees  of  the  valley  side  at  its  steepest  point  along 
an  orthogonal  profile.  This  definition  assumes  that 
short,  oversteepened  sections  of  the  valley  side 
near  channels,  that  are  the  obvious  result  of  recent 
rejuvenation,  are  not  included  in  the  profile.  It  is 
not  likely  that  such  sections  of  slope  would  be  in¬ 
cluded  in  map  measurements,  but  it  is  necessary 
to  exclude  them  from  the  definition,  as  valley-side 
slopes  are  also  measured  in  the  field  (see  maxi¬ 
mum  valley-side  slopes  under  Measured  Properties). 
The  measure  is  obtained  from  a  map  by  computing 
the  horizontal  distance  between  two  contour  lines 
that  are  separated  vertically  by  say,  100  feet,  and 
converting  the  tangent  to  the  angle  in  degrees.  De¬ 
pending  on  the  variability  of  slopes  within  the  basin. 
15  to  75  slope  readings  are  needed  to  give  an  ade¬ 
quate  estimation  of  the  mean  value. 

The  importance  of  the  valley-side  slopes  instud¬ 
ies  of  erosion  lies  in  the  fact  that  the  energy  avail¬ 
able  for  erosion  in  foot-pounds  per  pound  of  water 
per  linear  foot  is  eoual  to  the  sine  of  the  slope. 

Dihedral  angle  between  valley  sides.  £.  meas¬ 
ured  in  degrees,  is  the  angle  between  the  valley 
sides  measured  perpendicular  to  the  stream.  It 
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cannot  Ik-  measured  easily  either  m  the  field  or 
from  the  m;»|»,  IjuI  It  ran  be-  calculated  Figure  2 
is  an  illustration  of  tin-  dihedral  angle.  inflt-rllon 
angh  gradient.  and  valley-side  slope  Appendix  I 
gives  Hu-  geomet rlral  relation  among  the  quan'ltles 
l)  ||/  (  and  > 

Gradient  of  stream  channel,  > ,  expressed  in  de¬ 
crees.  can  be  obtained  from  Hie  map  by  measuring 
the  distance  along  the  channel  between  two  or  more 
contour  lines.  This  quantity  is  of  great  Importance 
in  the  equilibrium  of  the  Ixisin,  since  It  is  related 
to  rates  of  transport,  runoff,  and  celerity  of  flood 
<  rests . 


figure  d .  M'-lallons  among  valley-side  slope.  Inflection 
angle,  channel  gradient,  and  dihedral  i  igle  between 
vilh-y  sides  lor  an  Idealized  channel. 


Problems  ol  Field  Study 


METHODS  OF  FIELD  OBSERVATION 

Field  Investigations  were  carried  out  (1)  to  ob¬ 
serve  the  accuracy  of  the  map.  particularly  in  re¬ 
gard  to  the  extent  of  the  drainage  net.  and  (2)  to 
obtain  measurements  of  characteristics  of  the  soil 
and  vegetal  cover  that  are  believed  to  Influence 
morphometric  properties  of  the  Ixtsin 

At  representative  points  within  each  basin 
several  stations  were  established  at  which  a  ser¬ 
ies  of  measurements  and  observations  was  made. 

In  large  basins  about  10  stations  per  square  mile 
were  established,  in  small  basins  the  density  in¬ 
creased  to  200  stations  per  square  mile.  To  esti¬ 
mate  the  number  of  measurements  required  to  de¬ 
tect  a  difference  of  a  specified  amount  in  some 
property  of  interest,  among  a  number  of  areas, 
a  prior  knowledge  of  the  standard  deviation  of  the 
property  is  necessary.  Estimates  of  the  standard 
deviations  of  the  properties  studied  in  the  field, 
excepting  maximum  valley-side  slopes,  were  not 
available  at  the  beginning  of  the  study.  Without 
such  estimates  it  is  impossible  to  judge  in  ad¬ 
vance  the  necessary  number  of  readings  within 
each  basin,  or  the  number  of  basins,  needed  to  de¬ 
tect  differences  in  means  of  the  properties  studied. 
It  was  assumed,  however,  that  for  a  given  property, 
greater  variability  would  exist  among  the  means 
for  all  basins  than  among  readings  within  the 
basins;  analysis  of  the  data  proves  this  assumption 
to  be  correct  (Table  6).  Therefore,  the  emphasis 
was  placed  on  gathering  data  from  a  large  number 
of  basins,  and  fewer  measurements  were  taken 
within  each  basin. 

A  major  problem  in  the  field  study  of  each  area 
was  to  decide  which  linear  depressions  should  be 
included  as  channels,  and  which  should  be  omitted 
The  crux  of  the  problem  is  in  the  identification  of 
the  smallest  “finger-tip,”  or  first-order  channels, 
because  the  channels  of  higher  order  present  no 
difficulty  once  all  first-order  channels  are  drawn 
on  the  map.  First-order  channels  can  be  defined 
as  the  smallest  permanent  depressions  found  in  a 
connected  drainage  system,  whose  form  is  main¬ 
tained  indefinitely  against  agents  tending  to  obliter¬ 
ate  it,  by  the  action  of  water  flowing  in  a  concen¬ 
trated  stream.  However,  any  definition  depending 
on  the  knowledge  of  the  behavior  of  the  depression 
through  several  years,  or  seasons,  is  unusable  in 
studies  in  which  observations  are  made  in  each 
basin  during  a  single,  relatively  brief  period.  For 


the  purposes  of  field  identification,  true,  functional, 
erosionnl.  li  st-order  channels  are  defined  as  linear 
depressions  having  the  following  characteristics: 

(1)  Water  entering  the  depressions  is  forced  to 
flow  in  a  concentrated  stream,  and  under  no 
conditions  can  become  sheetflow. 

(2)  At  its  lower  end  the  depression  (a)  joins 
another  similar  depression  of  equal  or 
greater  size,  (b)  disperses  onto  an  alluvial 
fan  or  cone,  or  (c)  ends  at  a  sinkhole  or 
other  localized  area  of  recognizably  great 
permeability. 

(3)  Erosion  and/or  creep  have  produced  (a) 
graded  valley  sides  that  lead  down  to  the 
depression,  and  have  orientations  noticeably 
different  from  the  hillside  adjoining  the 
area  drained  by  the  depression.  The  gradient 
of  the  depression  in  this  case  is  lower  than 
the  slope  of  the  adjacent  surfaces  (Gilbert, 
1877,  p.  115).  Or  (b)  the  depression  is  steep¬ 
sided  with  rectangular  or  trapezoidal  cross- 
section.  its  depth  being  on  the  order  of  ten 
times  the  maximum  micro- relief  on  the 
older  unmodified  surfaces  near  the  depres¬ 
sion. 

This  definition  excludes  steep-gradient,  grassy 
swales  through  which  water  might,  after  a  mild 
rain,  pass  in  the  form  of  sheetflow.  However,  per¬ 
haps  such  channels  need  not  be  excluded  from  con¬ 
sideration  as  higher  order  channels,  uncommon 
though  they  may  be. 


OBSFRVED  WATERSHED  PROPERTIES 

Vegetation:  The  typical  vegetation  assemblages 
with  their  elevations  and  exposures  were  noted  in 
each  basin  (Table  3).  Identifications  of  species 
were  made  from  Southwestern  Trees  (Little,  1950). 
Certain  species  of  trees  and  bushes  provide  a 
basis  from  which  inferences  about  the  climate  of 
the  basin  can  be  drawn  (Little,  1950;  Costello,  19?). 
Estimation  of  the  precipitation-effectiveness  index 
for  each  basin  was  obtained  from  a  regression  line 
for  precipitation-effectiveness  against  elevation 
and  confirmed  by  field  data.  For  instance,  the 
saguaro  cactus,  creosote  bush,  and  paloverde  of  the 
desert  or  Lower  Sonoran  Life  Zone  indicate  a  pre¬ 
cipitation-effectiveness  index  of  less  than  15;  the 
shrub  live  oak  and  cliffrose  of  the  chaparral  or 
Upper  Sonoran  Life  Zone  indicates  a  P-E  index  of 
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35  t»  40  ami  the  Kngelmann  spruce  and  alpine  t tr  of 
tfie  llud.soniun  Life  Zone  indicate  a  P- I  index  of  90 
to  120  (see  Climatic  Data). 

Lithology  Bnsii.if  were  selected  for  study  partly 
on  ifie  Ixtsis  of  similarity  of  lithology  throughout.  In 
a  few  instances,  tiasins  were  found  to  !»»'*• 
i  i  i  n  oiite,  tug  from  that  of  the  major  part  of  the 
basin  In  addition,  in  four  lias  ins  the  rocks  lie  In 
near-horizontal,  alternating  weak  and  strong  layers, 
consisting  of  sandstone,  limestone,  and  shale,  or 
lava  and  pumice.  Descriptions  of  the  lithology  pre¬ 
vailing  in  each  basin  are  given  in  the  general  de¬ 
scription  of  basins.  If  consistent  structural  planes 
are  present  throughout  the  area,  their  strike  and 
dip  are  noted  also. 

Percentage  of  basin  lying  within  each  type  of 
cover:  Each  basin  studied  in  the  field  was  divided 
into  portions  covered  by  forest  and  bushes,  grass¬ 
land.  and  bare  soil  or  rock.  The  percent  of  total 
area  lying  within  each  type  of  cover  was  then  esti¬ 
mated.  either  from  aerial  photo  raphs,  or  from  a 
vantage  point  providing  a  general  view  of  the  entire 
basin.  It  was  found  that  on  those  maps  with  green 
overprint  representing  vegetative  cover,  there  was 
often  little  corresponded  e  with  the  actual  cover. 
Data  are  given  in  the  general  description  of  each 
basin. 

Range  evaluation:  Many  of  the  basins  studied  in 
the  field  were  being  grazed  at  the  time  the  field 
work  was  done,  or  had  recently  been  grazed.  The 
extent  of  grazing  was  estimated  from  observation 
of  the  general  length  of  grass  stubble.  A  simple 
count  of  animal  density  is  inadequate;  for  instance, 
near  Tucson,  three  animals  per  square  mile  result 
in  considerably  over-grazed  range,  whereas  in  the 
Colorado  Front  Range,  many  more  animals  can  be 
supported  without  serious  over-grazing.  Data  for 
each  basin  are  given  in  the  general  description  of 
basins. 

Presence  of  gullies:  It  is  assumed  here  that  gul- 
lies  usually  signify  readjustment  of  the  drainage 
pattern  to  a  changed  environment.  Where  numerous 
gullies  were  present,  excluding  those  few  caused 
by  some  factor  not  affecting  the  entire  basin,  note 
was  made  of  whether  they  were  primarily  on  the 
valley  slopes  or  in  the  channels.  In  a  single  case 
extensive  gullying  was  found  in  both  situations.  It 
seems  likely  that  channel  trenching  and  gu.Uying  on 
the  upland  slopes  may  be  due  to  opposite  causes. 

The  surface  characteristics  of  the  mantle  of  the 
basins  may  have  changed  greatly  with  the  advent  of 
human  activity.  Grazing,  logging,  burning,  and  min¬ 
ing  activities  leave  very  apparent  effects.  Grazing 
can  drastically  reduce  the  amount  of  cover,  and 
even  moderate  grazing  lowers  the  infiltration  capac¬ 
ity  of  the  soil.  Watershed  management  studies  in  the 
Colorado  Front  Range,  on  granitic  soils  similar  to 
the  majority  of  those  studied  here,  and  in  ponderusa 
pine  forest,  show  that  heavy  grazing  more  than 


doubles  the  amount  of  erosion  per  unit  area,  for 
each  Inch  of  runoff  (lam',  1953,  p  2H)  It  coultl 
therefore  l>e  contended  that  measurement  of  surface 
characteristics,  Including  vegetal  cover  as  they 
are  now  will  not  reflect  the  condition  of  the  l»<xo< 
v.t.c ,,  ,t  w.ts  oev eloped,  and  are  therefore  valueless 
All  evidence  available  to  the  writer  indicates  that 
Increase  of  dm  Huge  density,  brought  about  by- 
lowered  resistance  to  erosion,  concentration  of 
runoff,  decreased  Infiltration  rate,  disruption  of 
plant  cover,  etc.,  occur  very  soon  after  the  change 
of  conditions.  New  gullies  on  slo|X)8  lie  low  a 
new  road,  or  in  areas  of  recently -burned  forest 
arc  convincing  evidence  of  the  rapidity  of  the 
change.  Gullies  In  slopes  are  Included  in  the  drain¬ 
age  net  and  it  is  therefore  considered  likely  that 
measurements  of  surface  properties  do  bear  rela¬ 
tion  to  the  extent  of  channels.  Part  of  the  difficulty 
is  in  the  assumption  inherent  in  the  objection.  Ba¬ 
sins  were  not  developed  then  left  unchanging  and 
static  for  the  next  surface  geologist  to  study,  but 
must  reflect  by  continual  development  the  changes 
in  climate  and  environment  within  the  basin. 

Perhaps  there  Is  more  reason  to  believe  that 
degiee  of  slope  will  not  be  influenced  immediately 
by  changes  in  surface  properties,  as  time  is  re¬ 
quired  to  accomplish  the  regrading  of  a  slope. 
Erosion  of  a  large  amount  of  material  is  necessary 
to  change  the  average  basin  slope  even  a  degree 
or  two.  Since  rates  of  creep  are  not  well-known, 
supporting  evidence  for  this  view  is  needed. 


GENERAL  PLAN  OF  FIELD  MEASUREMENTS 

At  each  station  established  in  a  basin,  a  series 
of  measurements  was  repeated  in  as  nearly  the 
same  manner  as  possible  from  basin  to  basin.  All 
measurements  should  have  been  taken  at  randomly 
located  points,  but  unfortunately,  it  was  impossible 
to  do  this  consistently.  In  cobbly  and  bouldery  soils 
it  was  often  difficult  to  locate  places  on  which 
several  of  the  measurements  could  be  made.  In 
other  areas,  such  as  the  Chinle-shale  badlands, 
where  full  randomization  would  be  relatively  easy, 
it  was  less  important  dut  to  the  high  degree  of 
homogeneity  of  the  area. 


MEASURED  PROPERTIES 

Infiltration  capacity:  The  present  study  did  not 
include  a  general  investigation  into  the  factors  in¬ 
fluencing  the  infiltration  capacity  of  the  soil.  The 
infiltration  capacity  measurements  were  taken 
solely  with  the  intention  of  obtaining  comparative 
data,  not  necessarily  the  actual  rate  of  infiltration 
during  a  rain. 

The  infiltrometer  used  in  this  study  was  designed 
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l>V  the  wrlli'r  U>  lx-  portable,  light  in  weight,  and 
conservative  in  use  of  water.  These  requirements 
were  important  Ix-rause  In  many  areas  it  was  nei - 
essnry  to  currv  the  apparatus  and  water  r 
i  ugged  ter  rain. 

The  principle  of  measurement  Is  simple  water 
from  a  tank  Is  drawn  through  a  metering;  device, 
an  adjustable  valve,  a  sprinkler,  and  distributed  by 
Inspection  evenly  over  a  plat  of  soil  one  squat  ■?  foot 
in  area  at  the  maximum  rate  of  infiltration.  Dis¬ 
charge  Is  measured  by  the  meter  In  ounces  per 
minute,  and  Is  converted  to  Inches  (depth)  per  hour 
over  an  area  of  one  square  foot.  The  apparatus  Is 
mounted  on  a  standard  surveying  tripod,  and  weighs 
•15  pounds  when  full  of  water  (Figure  3).  The  reser¬ 
voir  tank  holds  2  1  2  gallons,  and  will  supply 
enough  water  for  two  to  five  tests,  depending  on  the 
length  of  run,  and  the  infiltration  capacity  of  the 
soil. 


Figure  3.  Inl'i Urometer  assembly  diagram. 

The  one  square-foot  area  is  bounded  by  a  stain¬ 
less  st€nl  ring  with  a  sharpened  lower  edge,  that 
can  be  driven  into  the  ground.  Doubtless  a  double¬ 


ring  arrang'-nu  nt,  with  water  filling  an  t  jtside  ring 
to  lessen  th  •  amount  of  lateral  flow,  would  lx-  de¬ 
sirable  were  accuracy  the  main  consideration.  The 
added  weight  ,.f  .ti.wiin *  i  ami  the  nut. cornu  11  y 
greater  quantity  of  water,  would  reduce  the  distance 
which  the  instrument  could  bo  carried  during  a  day  s 
work.  The  lnflltrometor  used  in  this  study  Is  defin¬ 
itely  a  reconnaissance  Instrument,  and  cannot  com¬ 
pare  In  accuracy  with  more  elaborate  equipment. 

The  tank  containing  the  water  is  made  of  stain¬ 
less  steel,  and  Is  built  on  the  principle  of  the  Marl- 
otle  bottle  (Westphal,  1952,  Z  welter  Tell,  p.  15). 

The  tank  can  lx?  closed  air  tight  by  a  rubl>er  stopper, 
the  only  air  entering  the  tank  comes  through  a 
copper  tube  extending  down  almost  to  the  bottom  of 
the  tank.  Thus,  the  only  free  surface  Is  at  the  bot¬ 
tom  of  the  tube,  and  once  the  air  has  started  bubbling 
out  of  the  tube  into  the  tank,  the  pressure-head  is 
measured  from  the  end  of  the  tube.  Therefore, 
throughout  the  course  of  the  test,  though  the  level 
of  the  water  decreases,  the  hydrostatic  head  re¬ 
mains  the  same.  This  is  necessary  because  tne 
measuring  instrument  is  accurate  under  constant 
head.  Once  a  definite  rate  of  flow  has  been  estab¬ 
lished,  it  is  not  necessary  to  make  repeated  changes 
in  the  valve  setting,  as  would  be  the  case  were  the 
head  continuously  decreasing. 

The  measuring  device  is  a  Flowrator,*  which 
is  basically  a  tapered  glass  tube  with  a  metal 
float  fitting  closely  within  the  tube.  For  protec¬ 
tion,  a  pressed  steel  cover  surrounds  the  tube 
and  supports  the  pipe  connections.  Water  flows  up¬ 
ward  through  the  tube;  as  the  discharge  increases 
the  frictional  force  of  the  water  on  the  tloat  in¬ 
creases.  and  an  upward  acceleration  is  imparted. 
Higher  in  the  tube,  the  diameter  of  the  opening 
and  the  area  of  discharge  are  greater;  the  vel¬ 
ocity  of  water  flowing  around  the  float  is  less,  and 
hence  the  frictional  force  is  less.  The  float  there¬ 
fore  seeks  a  level  at  which  its  weight  is  exactly 
balanced  by  the  upward  force  of  the  moving  water, 
and  this  position  is  a  measure  of  the  discharge. 

The  maximum  discharge  possible  is  11.2  ounces 
per  minute  or  8.4  inches  ppr  hour.  Figure  4  gives 
the  conversion  graph  from  percent  of  total  flow  to 
inches  per  hour  over  one  square  foot.  After  leav¬ 
ing  the  Flowrator,  the  water  passes  down  through 
a  sprinkler  and  is  distributed  evenly  over  the  area 
within  the  ring. 

The  surface  of  the  test  area  is  prepared  by  re¬ 
moving  loose  leaves,  twigs,  accumulated  pine  duff 
and  needles,  leaving  only  grass  on  the  largely  min¬ 
eral  soil.  At  the  beginning  of  a  test,  the  discharge 
is  set  at  75%  of  the  maximum,  8.4  ounces  per  min¬ 
ute,  or  6.3  inches  per  hour.  If  the  soil  is  dry,  it 
will  take  15  to  100  seconds  to  saturate  the  surface 

♦Manufactured  by  Fisher  and  Porter  Company,  Hatboro, 
Pennsylvania. 
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Figure  I.  Conversion  graph:  percent  of  total  flow  (100f 
11.2  ounces  per  minute)  to  inches  per  hour  over  an  area 
of  one  square  foot. 

and  to  form  excess.*  The  time  at  which  the  excess 
forms  is  recorded,  and  the  discharge  is  decreased. 
Some  practice  is  required  to  adjust  the  discharge 
to  match  the  infiltration  rate,  and  at  best  it  takes 
one  minute  to  determine  if  the  machine  is  adjusted 
correctly.  On  most  poorly-developed,  gravelly 
soils,  during  the  first  15  minutes  of  the  test,  the 
rate  of  infiltration  decreases  rather  rapidly. 
Thereafter,  the  rate  is  fairly  constant,  though  it 
may  decrease  slowly  if  considerable  clay  is 
present  in  the  soil.  It  is  usually  possible  to  adjust 
the  instrument  to  within  1%  of  flow  after  20  to  25 
minutes  have  elapsed,  giving  an  adjustment  error 
on  0.084  inches  per  hour,  and  the  rate  of  infiltra¬ 
tion  usually  changes  little  thereafter.  The  “instan¬ 
taneous”  rates  for  3,  5,  10,  15,  20,  25,  30  and 
sometimes  40  and  50  minutes  were  recorded.  The 
final,  stable  rate  (usually  the  minimum  rate)  was 
taken  as  the  infiltration  capacity  (Horton,  1942, 
p.  480).  The  values  of  infiltration  capacity  obtained 
with  this  instrument  are  probably  higher  than  the 
actual  value  prevailing  during  a  hard,  general  rain, 
since  no  allowance  is  made  for  the  lateral  flow  of 
water  when  it  passes  below  the  edge  of  the  ring. 

+  L.  1).  I, uve  (personal  communication,  1955)  ami  others 
believe  that  because  of  wetting  difficulties,  most  reliable 
results  may  be  obtained  on  pre-wetted  plats  that  have 
drained  12  to  24  hours.  The  writer  encountered  wetting 
difficulties  commonly  in  soils  with  a  deep  pine-needle 
deposit.  In  a  few  cases  a  small  amount  of  detergent  was 
added  to  the  water,  and  this  usually  shortened  the  time 
needed  to  saturate  the  surface.  Too  tew  trials  were  made 
to  determine  whether  this  altered  the  final  infiltration 
capacity. 


A  greater  source  of  error,  however,  may  bo  the 
illnturlami  of  the  soil  structure  caused  by  pound 
lug  In  the  ring,  thereby  opening  up  fissures  not 
otherwise  present.  The  values  obtained  are  believed 
lo  lx-  consistent  within  themselves,  since  \  nines  of 
infiltration  capacity  vary  greatly  between  areas  of 
different  rock  type,  and  less  so  within  areas  of 
similar  type. 

Soil  strength,  s,  Is  measured  under  instantan¬ 
eous  Impact  of  a  12  pound  Iron  shot,  dropped  a 
standardized  distance,  lo  Impart  CO  foot  pounds  of 
energy  normal  to  the  ground  surface.  This  distance 
Is  5  feet  for  horizontal  surfaces,  and  5  feet  divided 
by  the  cosine  of  the  slope  on  sloping  surfaces.  Sat¬ 
isfactory  results  can  l>e  obtained  on  slopes  up  to 
35  degrees.  The  volume  of  the  imprint  of  the  shot 
div  led  Into  60  foot-pounds  gives  a  measure  of 
strength  in  pounds  per  square  foot. 


60  ft. -lbs, 
x  ft.1 


ou  ..  : 

—  lbs.  per  ft. 
x 


The  volume  of  the  imprint  can  lx?  obtained  from  a 
graph  relating  diameter  of  the  base  of  a  spherical 
segment  to  volume.  The  diameter  of  the  imprint  is 
easily  measured.  Appendix  III  gives  the  derivation 
of  the  relation  between  the  diameter  of  imprint  and 
volume  of  imprint.  Figure  5  illustrates  the  method: 
Figure  6  is  a  grapli  relating  diameter  of  imprint 
directly  to  strength  of  soil. 


Figure  5  Principle  of  meusuring  soil  strength  by  drop¬ 
ping  a  spherical  weight  a  standardized  distance  and 
measuring  the  volume  of  imprint. 
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Figure  6.  Conversion  ol  diameter  ol  imprint  ol  12 -Ih. 
shot,  diameter  0.37b  It.,  to  soil  strength,  when  dropped 
he  distance  necessary  to  impart  fiO  ft. -lbs.  normal  to 
the  soil  surface. 


To  standardize  the  moisture  conditions  under 
which  the  strength  tests  were  maue,  two  or  three 
repetitions  of  the  test  were  made  on  the  soaked 
plat  of  ground  used  in  the  infiltration  test  to  obtain 
a  value  of  the  strength  of  the  soil  approximating 
its  condition  after  a  hard  rain.  In  addition,  three 
or  four  tests  were  made  in  the  surrounding  area 
on  the  soil  in  its  natural  condition.  In  most  areas 
the  soil  was  extremely  dry;  only  in  a  few  cases 
was  it  damp. 

It  would  have  been  desirable  to  randomize  all 
measurements  of  strength.  However,  the  measure¬ 
ment  of  wet  soil  strength  was  necessarily  limited 
to  the  area  selected  for  infiltration  measurements, 
which  were  in  turn  limited  by  the  stoniness  of  the 
ground.  Randomization  over  dry  surfaces  would 
produce  tests  on  cobbles,  logs,  etc.,  and  the  re¬ 
sults  would  not  be  indicative  of  soil  strength.  The 
measurements  were  probably  not  systematically 
biased,  as  there  was  no  way  of  estimating  the 
strength  of  the  point  directly  under  the  shot  prior 
to  making  the  test. 

Size-distribution  of  soil  was  described  at  each 
station  within  the  infiltrometer  ring  by  the  visual 
and  manual  methods  developed  by  the  Department 
ot  Soil  Mechanics  at  Columbia  University  (Bur- 
mister,  1952,  p.  48).  The  average  soil  descrip¬ 
tion  for  test  plats  within  each  basin  is  given  in 
Table  4. 

Percent  bare  area  was  obtained  by  measuring 
the  exposed,  bare,  mineral  soil  with  a  steel  meas¬ 
uring  tape  extended  for  50  feet  over  the  ground  in 
a  randomly  oriented  direction.  The  ground  at  each 
foot  mark  on  the  tape  was  inspected,  and  if  bare  of 
all  cover,  such  as  grass,  fallen  leaves,  twigs, 
limbs,  cobbles,  or  larger  rock  fragments,  it  was 


ghm  a  ciiuni  of  one.  For  rarti  fool  mark  lying 
atxnr  a  fieri  Ion  of  ground  lltat  covered  in  an> 
way,  the  count  wa.n  zero,  ojkxi  'ountmg  through  SO 
cel,  She  iroiiilv;  »,{  fi/wt  it  .f  'ha  1  log  a  Ur.  c  luir 
ground  wan  multiplied  y  .wo  to  give  an  expression 
in  percent  of  total  area. 

lhe  direction  of  traverse  am*  |*itnt  «.!  origin 
were  randomized  by  throwing  a  geologic  liamim  r  in 
the  air  with  a  spinning  motion,  and  extending  the 
tape  in  the  direction  In  which  the  pick  pointed  upon 
landing.  Two  such  readings  were  made  at  cacti  .sta¬ 
tion.  The  average  of  the  results  for  each  basin  arc 
given  m  Tabic  3. 

Roughness,  or  numln  r  and  size  of  nek  'r.ign.cnts 
lying  loose  on  the  surface,  or  partially  burled,  was 
found  by  counting  and  measuring  all  fragments  1  2 
inch  or  greater  in  mid-  ...i meter  within  a  circle  with 
a  radius  of  one  f (Kit .  The  center  of  the  circle  was 
located  by  selecting  a  numlxT  of  feet  between  9  and 
100  from  a  table  of  random  numbers,  an  orientation 
was  found  by  throwing  a  geologic  hammer  into  the 
air  in  the  same  manner  as  dcsrrit>od  above  lor  per¬ 
cent  bare  area.  Fragments  were  measured  to  the 
nearest  1  16  inch  by  small  calipers.  The  distance 
which  a  buried  fragment  extended  above  the  ground 
was  recorded.  One  circle  was  measured  at  each 
station.  The  roughness  number,  M,  is  the  average 
total  length  of  diameters  of  pebbles  in  each  circle 
for  the  entire  basin.  M  is  then 


M 


\ 


2-  f.m, 

i=  1 

N 


\ 

V 


»-  1  1 
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where  f,  is  the  number  of  rock  fragments  in  the  i’th 
circle  of  one  foot  radius,  m;  is  the  average  mid¬ 
diameter  of  rock  fragments  in  the  i’th  circle,  x  ,  t 
is  the  diameter  of  the  j’th  pebble  in  the  i’th  circle, 
and  N  is  the  number  of  circles  measured  in  each 
basin.  The  roughness  number  increases  rapidly 
with  increasing  number  of  loose  fragments,  and 
also  with  increasing  average  size  of  rock  fragments. 
Among  the  basins  studied,  those  having  a  large 
number  of  smaller  fragments  indicate  greater 
roughness  than  basins  with  a  small  number  of  lar¬ 
ger  fragments. 

Maximum  valley-side  slopes,  O,  were  measured 
in  the  field  with  an  Abney  hand  level,  calibrated  to 
read  to  1/2°.  Within  each  basin,  about  40  measure¬ 
ments  were  usually  taken,  which  would  provide  suffi¬ 
cient  probability  of  detecting  differences  of  1  1/2°  or 
2°  in  the  means  of  two  basins  when  the  standard  devi¬ 
ation  is  about  4°.  The  mean  slopes  of  all  basins  for 
which  data  were  taken  are  given  in  Table  2. 

In  measuring  slopes,  the  writer  walked  along  the 
valley  side  near  the  upper  limit  of  the  steepest 
portion  of  the  slope,  with  an  assistant  near  the  foot 
of  the  slope.  At  fairly  regular  intervals,  slope  and 
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ur tmialion  readings  ttpfr  recorded  In  Itasina  in 
»hlch  t  re  nr  hi  tic  or  rejuvenation  ol  the  channel  was 
4|ijnfpiit,  I  hi-  short,  sleep  slope  near  lit-  channel 
•a as  i^iitn t*t|,  and  the  long  slope  on  the  main  valley 
stil<-  was  mi-aiurnl. 

In  lar^e  liaslns.  (host-  o|  om-  square  milt-  or 
i;  r  »-:«*«•  r  an  •>  s  system  ul  randomized  sampling 


that  will  slill  alio-*  rapid  measurement  Is  needed 
In  such  cases  a  tendency  lo  Ignore  slopes  ol  lower 
aiijilf.  and  lake  loo  many  readmes  on  sleep 
slopes  would  Idas  l  lx  mean.  In  smaller  lias  ins, 
where  several  readmes  are  laken  on  each  valley 
Side,  Idas  is  less  likely  •’*•!  aune  a  ^r«at<  i  pc i  cent 
ol  llie  area  is  sampled. 


1. lunatic  Data 


lllllr  i  nricfill.il  aih.inic  call  l»-  m.i<!c  in  llir  in\ iiHj.;.in>n  <•!  rainfall  In  .  •  >iituuiali"n  of  I  *-a  llngn  at  arbitral 
ilv  nclccti  il  liK'.lllll«-«  »hl<  h  ha|>|>ctl  l<>  !.«•  icfilcrn  *>(  |x>|<ulati<>n  lull'll  nr  Mull  of  r  1 1 11 1 .1 1 1  lor  a  rrl.iltwh  nlc  ft 
|)crlml  at  cardulh  nrlcctcil  critical  loi'.illtli  n  mnilil  .nli.wnc  our  l.noii  Imljn  lir  iivf  than  tile  c<  •fllimi-’ltl  'll  !"I 
many  years  of  routine  nunllng*  at  n|><>r:»itear  .  locate. I  nt.iti.mn  (Turn.ij;e  att.i  Mallcrv  lhll  from  lrilro.iuction. 
In  Forest  Shreve) 


Since  1941  there  has  been  no  advance  in  system¬ 
atic  collection  of  weather  records,  especially  tem¬ 
perature  data,  in  mountainous  regions.  In  compiling 
climatic  data  for  tills  study,  considerable  reliance 
was  placed  on  records  (some  of  only  a  few  years 
duration)  obtained  from  high  mining  tamps  in  the 
Colorado  Front  Range  (1860  to  1915),  from  ranches 
in  New  Mexico  and  Arizona,  ami  from  Forest  Service 
ranger  stations  and  fire  lookouts.  Perforce,  much 
of  the  data  is  from  the  yeats  before  19.10,  and 
obtained  from  Bulletin  W,  Climatic  summary  of  the 
United  States  (U.  S.  Weather  Bureat  ,  Department 
of  Agriculture,  1933).  Although  the  Weather  Bureau 
under  the  Department  of  Commerct  has  placed  em¬ 
phasis  on  improving  techniques  of  lorecasting, 
there  is  a  great  need  for  research  into  the  effects 
of  topography,  exposure,  and  elevation  on  precipi¬ 
tation  amounts  and  intensities,  temperature,  and 
snow  accumulation,  in  connection  with  hydrologic 
studies  for  flood  control,  erosion  problems,  forest 
and  range  management,  and  other  conservation 
programs.  This  need  is  being  met  partially  by 
utilization  of  weather  records  of  Forest  Service 
Experimental  Stations,  and  by  regional  studies 
such  as  those  prepared  bv  the  Hydrologic  Services 
Division,  U.  S.  Weather  Bureau  (U.  S.  We  ither 
Bureau  Technical  Paper  No.  24,  1954,  Parts  I  and 
II).  But  these  summaries  would  be  greatly  im¬ 
proved  were  adequate  mountain  weather  data  avail¬ 
able. 

The  data  secured  consist  of  average  monthly 
precipitation,  monthly  average  of  mean  daily  tem¬ 
perature,  average  number  of  days  per  year  with 
.01  inches  or  more  precipitation,  and  elevation  of 
the  station.  In  addition,  the  depth  of  precipitation 
of  the  average  five-year,  one-hour  storm  for  the 
general  region  was  obtained  from  Technical  Paper 
No.  24,  Part  II  (U.  S.  Weather  Bureau,  1954).  Data 
are  presented  by  state,  section,  and  station  in 
Table  5. 

The  Thornthwaite  Precipitation- Effectiveness 
Index,  or  P-E  index,  is  a  measure  of  the  availabil¬ 
ity  of  moisture  to  vegetation,  and  depends  on  the 
amount  and  distribution  of  precipitation  and  evapor¬ 
ation;  the  latter  in  turn  depends  on  variations  of 


temperature  and  wind  movement  (Thornthwaite. 
19.11).  The  exact  rational  formula  is 

i  io  LJl 
i.  r  ■ 

where  P  is  the  average  precipitation  for  each 
month,  and  E  is  the  average  evaporation  for  each 
month.  The  index  I,  is  therefore  a  dimensionless 
ratio.  An  empirical  formula  that  is  an  approxima¬ 
tion  to  I.  when  no  evaporation  data  are  a  ailable, 
is 


where  T  is  the  monthly  average  daily  temperature. 

Thornthwaite  divided  the  range  of  values  of  1  as 
follows: 


Characteristic  vegetation  P-E  Index 

Rain  forest  ?  128 

Forest  64  -  127 

Grassland  32  -  63 

Steppe  16  -  31 

Desert  •  15 


The  P-E  index  for  each  station  for  which  sufficient 
data  are  available,  or  for  which  temperatures  can 
be  estimated  by  interpolation  from  nearby  stations, 
is  provided  in  Table  5. 

In  the  areas  of  the  West  for  which  data  were 
obtained,  the  precipitation  increases  with  eleva¬ 
tion,  above  a  lower  limit,  and  temperature  decreases 
(Turnage  and  Mallcry,  1941).  Within  a  limited  area, 
the  P-E  indices,  when  plotted  against  elevation, 
show  a  significant  linear  increase.  The  points  on 
the  graphs  deviate  from  linearity  in  Colorado  and 
Utah  for  elevations  greater  than  10. DUO  feet  and 
8,000  feet  respectively.  The  graphs  obtained  were 
then  used  to  estimate  the  mean  P-E  index  for  each 
basin.  Ignored  are  the  differences  of  precipitation 
effectiveness  on  north  and  south-facing  hillsides, 
the  effect  of  nearby  mountains,  arid  the  situation  of 
the  station.  The  value  of  I  is,  therefore,  only  ap¬ 
proximate,  However,  the  results  agree  well  with 
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W£«*l4Ui»n  til  Ihr  lunlltn  Figure  7  •>>><••*  r» 
t‘f4|jf’.n  f  I’-F  Index  4g4lnnt  e|cS4ltoll  * 


P  F  Index 


Figure  7  (a) 

P-  E  Index 


Figure  7  (b) 

♦Lines  or  curves  fitted  by  eye  only. 


P-E  Index 


Figure  7  (c) 


Figure  7.  Variation  of  precipitation— effectiveness  index 
with  elevation; 

(a)  in  Park  Park  Plateau  and  Front  Range,  Colo,  and 
Wasatch  Mtns.,  Utah, 

(b)  in  southeastern  Arizona  and  southwestern  New  Mex¬ 
ico, 

(c)  in  north-central  and  west-central  New  Mexico. 

Intensity-du  rat  ion- frequency  relations  are  ex- 
pressed  by  the  amount  of  rain  received  from  a 
storm  of  one  hour  duration  expected  to  occur  once 
in  five  years.  The  information  was  obtained  from 
Technical  Paper  24,  Part  II  (U.  S.  Weather  Bureau, 
1954,  p.  6-8),  by  multiplying  the  two-year  one-hour 
amounts  as  given  on  isopleth  maps,  by  1.4  which 
is  the  ratio  of  the  amounts  for  the  two  return  per¬ 
iods  in  the  western  United  States. 

Season  of  maximum  rainfall  is  secured  directly 
from  the  tabulated  amounts  of  precipitation  for 
each  month. 

Ratio  of  mean  January  precipitation  to  1/ 1 2 
mean  annual  precipitation  is  obtained  directly  from 
the  values  for  mean  annual  precipitation,  and  mean 
January  precipitation.  The  measure  is  important 
because  of  the  greater  effectiveness  of  rain  for 
erosion  when  the  vegetal  cover  is  least  extensive 
as  it  is  in  this  winter  month.  The  calculated  value 
was  found  to  be  within  the  range  of  the  value  given 
by  the  map  in  Climatic  Atlas  of  the  United  States 
(Visher,  1954,  p.  204). 
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General  Description  ol  Basins 


I  liih  M'l  tton  irtMlM  grnrr.il  a.spiTlA  of  fill  h 
Ixi.sln.  Ohi.i lin'd  largely  from  map  and  field  InHpcr- 
lion,  the  description  is  provided  (or  reference 
onlv,  no  effort  Is  made  to  ilscuss  the  r  'n'rllrutton 
of  the  Information  on  eaeh  lias in  in  rega.d  to  the 
coni' lus Ions.  Location,  physiographic  province  and 
section,  typo  of  surroundings,  grazing  history  and 
range  evaluation,  lithology  and  structure,  presence 
of  gullying,  type  and  extent  of  cover,  and  other 
features  are  liest  explained  verbally,  since  most 
of  the  information  is  non-numer leal.  Cltmatologtc 
data  arc  included  in  the  form  of  estimated  precipi¬ 
tation-effectiveness  index,  and  the  distribution  of 
rainfall  by  months.  The  basins  are  listed  alpha¬ 
betically  by  state  and  quadrangle  in  which  they  arc 
located.  Maps  of  all  basins  are  given  in  Appendix 
IV. 


BASINS  IN  ARIZONA 

Chinle  Badlands.  Cameron,  including  two  majut 
basins  in  this  locality'  are  cove  red  by  a  spec  ial  map 
made  in  August,  1953,  by  E.  D.  Koons,  S.  A. 
Scbumm,  and  A.  N.  St  -aider,  at  the  scaleone  inch 
to  50  feet  (1:600).  The  topography  is  “micro-relief," 
that  is,  an  area  of  fully-developed,  mature  topo¬ 
graphy,  relatively  rugged,  but  composed  of  ex¬ 
tremely  small  sized  units.  The  area  is  located  at 
the  edge  of  the  Navajo  section  of  the  Colorado 
Plateau,  in  the  Painted  Desert.  The  lithology  is  the 
blue-gray  Chinle  shale,  in  thickly  platy  laminae  of 
silt  and  clay,  and  fine  sand,  containing  nodules  of 
gypsum  and  thin  plates  of  iron-stained  chalcedony. 
Beds  are  essentially  flat-lying.  No  vegetation  what¬ 
soever  is  present  anywhere  in  these  basins.  The 
surface  layer  of  solt,  unstable  sand  and  clay,  one 
to  four  inches  thick,  with  numerous  drying  cracks, 
has  been  produced  from  the  original  shale  by  swell¬ 
ing  of  the  clay  during  rains,  and  partial  desiccation. 
Any  disturbance  of  the  surface  produces  a  cascade 
of  the  crust  material,  exposing  a  spot  of  bare,  un¬ 
altered  shale,  and  burying  the  foot  of  the  slope 
more  deeply.  (See  also,  Schumm,  1 . '54,  p.  80). 

The  nearest  weather  station  at  Tuba  City,  about 
18  miles  north,  is  probably  in  a  meteorologically 
similar  area.  Of  the  yearly  rainfall,  46%  occurs  in 
the  months  July  through  October.  The  region  is 
definitely  desertic,  with  a  P-E  index  of  12.8. 

Finley  and  Adams  Canyon,  and  Sycamore  Canyon 
(Harshaw  Quadrangle):  Four  fourth-order  basins 


arc  included  in  this  sect  inn.  Field  work  was  n«d  de¬ 
tailed.  only  infiltration  measurements  ’.ieri  taken 
In  these  areas.  The  kisins  In  in  the  Patagonia 
Mountains,  In  the  Basin  and  Range  province,  Mexi¬ 
can  Highland  section.  Underlying  the  liaslns  is  a 
granite  complex  that  is  deeply  weathered  in  places, 
with  many  res'dual  boulders  and  dikes  standing  in 
relief.  No  gullying  is  present.  Vegetation  varies 
from  oak  chaparral  on  south-facing  slopes,  to  oak 
woodland  and  pinon- juniper  on  the  north-facing  and 
higher  slopes.  Grass  is  plentiful  in  open  areas  and 
is  ungrazed.  About  85%  of  the  area  is  in  forest  and 
bushes.  12%  to  14%  is  grassland,  and  the  remaining 
is  hare  soil  or  rock.  The  P-E  index  is  estimated  at 
40  to  65,  depending  on  elevation  and  exposure.  The 
nearest  source  of  weather  records  is  Nogales.  14 
miles  southwest,  where  68%  of  the  yearly  precipi¬ 
tation  is  in  July  through  October. 

Tributary  Basins  to  Harshaw  Creek  on  Red 
Mounta  i  (Mt.  Hughes  Quadrangle),  located  in  Secs. 
T6,  17,  20.  2l,  Twp.  22  S..  Ra.TF  F.  are  one  fourth 
order  basin  and  two  third-order  liasins.  Red  Moun¬ 
tain  is  very  near  the  boundary  between  the  Sonoran 
Desert  and  Mexican  Highland  sections  of  the  Basin 
and  Range  province,  and  stands  prominently  2,300 
feet  above  the  Sonoita  Creek  valley.  The  rock  is 
wholly  igneous;  a  thick  series  of  mineralized 
felsite  overlies  prophyritic  basalt.  Cliffs  form 
easily  in  both  series,  and  slide  rock  slopes  are 
common.  Several  mine  drifts  are  present,  but  none 
are  active. 

The  main  cover  in  the  basins  is  oak  chaparral 
and  brush,  grading  downward  into  scmi-dese.t 
shrub  and  grassland,  and  upward  into  pinon-juniper 
woodland,  and  near  the  top  of  Red  Mountain  on  north¬ 
facing  slopes,  larger  Arizona  white  oak.  An  esti¬ 
mated  50%  of  the  area  is  in  brush  and  chaparral, 

30%  is  grassland,  and  20%  bare  rock.  The  grassland 
area  is  grazed  lightly,  and  (lie  range  appeared  to  be 
in  fair  to  good  condition. 

The  P-E  index  is  estimated  to  be  35  to  45,  de¬ 
pending  on  elevation  and  orientation.  The  nearest 
weather  station  is  at  Patagonia,  three  miles  distant. 
There,  63%  of  the  yearly  precipitation  occurs  in  the 
months  July  through  October. 

Whipple  Basin  (Prescott  Quadrangle)  is  a  fourth- 
order  basin  in  the  northwest  quarter  of  Sec.  26, 

Twp.  14  N.,  Ra.  2  W.,  at  the  north  boundary  of  the 
Whipple  Military  Reservation,  in  the  Mexican  High¬ 
land  section,  Basin  and  Range  province.  The  basin 
lies  in  an  extensive  structural  valley,  three  miles 
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froit  ifn  imrll,  «•  i wl  of  (hr  Hrarinliaw  McantaitKi  It.<- 
immediate  rrtlrf  If)  (fir  .U'r.l  In  low,  mil  mol  r  (full 
250  lrf(  Tfir  rot  k  In  lit  •(!(«•  granite.  not  greatlv 
weathered  where  c\|)(inrtl,  and  oVe-rlaUl  (iV  llllll 
well-  |  lulu  f  .1  Iril  (NHiltlrry  .t  1  lu  V  lull* .  The  reiver  In  ne¬ 
tt  malt'd  to  In-  0?  I;  g  r.i.nril.liwl,  30%  llrrft  and  Irj.shrre. 
.mil  3{  lurr  roek.  Grazing  is  ii.txlrr.itr  and  (fir 
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Prescott,  atxiul  two  mill's  to  t he  southwest, 
probably  has  very  nearly  the  same  climate.  The 
P-K  index  is  *10.  and  the  months  of  greatest  pre- 
(  ipitntlon  are  July,  August,  September,  IX'cembor. 
January,  and  February,  with  7 *1 X  of  the  annual 
precipitation  occurring  In  those  six  months. 

Hot;  Canyon  Basin  (Sonolla  Quadrangle)  is  a 
four  Hi -order  basTiTTributary  to  Mop  Canyon,  and  ex¬ 
tends  across  Sec.  20,  Twp.  20  S.,  Ha.  16  E.  The 
area  is  in  a  dissected  alluvial  fan  on  the  east  flank 
of  the  Santa  Rita  mountains,  in  the  Mexican  High¬ 
land  section,  Basin  and  Range  province.  The  im¬ 
mediate  relief  is  500  feet,  but  the  mountains  rise 
sharply  to  the  west  to  a  height  of  4,600  feet  above 
the  basin,  doubtless  causing  increased  precipita¬ 
tion  from  maritime  tropical  air  masses  from  the 
south  and  southeast.  South-facing  slopes  are  al¬ 
most  entirely  grassland,  while  the  north-facing 
slopes  have  only  small  grassland  areas.  About  55% 
of  the  basin  is  grassland,  45%  is  oak  woodland, 
and  less  than  1%  is  bare.  Grazing  is  light;  the 
range  is  fair  to  good,  with  much  long  bunch  grass 
remaining.  Gullies  are  present  in  some  small 
channels  and  in  many  places  extend  into  fresh 
slopes. 

The  nearest  weather  records  are  from  the 
Crown  C  ranch,  one  mile  east,  where  precipitation 
data  are  available.  The  P-E  index  is  estimated  to 
be  35  to  40,  depending  on  exposure.  Greatest  pre¬ 
cipitation  is  in  the  months  July  through  October, 
when  62%  of  the  annual  rainfall  occurs. 

Basins  in  Saguaro  National  Monument  (Tucson 
Quadrangle,  1904):  Two  fourth-order  basins  near 
the  west  boundary  of  the  Monument  and  one-half 
mile  north  of  the  headquarters,  were  mapped  by 
the  writer  and  assistant  in  December,  1955,  on  a 
scale  100  feet  to  the  inch  (1:1,200).  The  area  is  in 
the  Sonoran  Desert  section,  Basin  and  Range  prov¬ 
ince.  The  basins  are  in  an  area  of  low  relief,  a 
dissected  pediment,  near  the  foot  of  the  Tanque 
Verde  Mountains.  Rock  type  is  mainly  chlorite  and 
sericite  schist,  with  injected  fine-grained  granite. 
Attitude  of  (he  planes  of  schistosity  is  variable, 
but  the  strike  averages  N  10°  E,  dip  25°  to  45°  W. 
About  60%  to  70%  of  the  area  is  bare  soil  and  rock, 
the  remainder  is  covered  by  saguaro  “forest.”  The 
prickly  pear  is  over-grazed,  with  three  animals 
(horses)  per  square  mile,  and  the  range  is  very 
poor.  As  the  result  of  over-grazing  and  disturbance 


of  the  nurf.ii  r  bv  III*-  animals'  tuxives,  liullicroun 
gutlicn  tunc  U  come  established  on  the-  ftllln Ulc.n . 
amt  4  few  nm.il!  channels  deepened.  Main  channels 
arr  broad  and  filled  with  sand  and  gravel. 

The  climate  data  for  Hits  area  were  taken  froit 
Tucson  Magnetic  Observatory  and  lh«  old  Vail  and 
P.inl.wio  records.  Tin  P-1  index  is  estimated  at 
15.  Of  the  annual  rainfall,  52%  is  In  the  months 
July  lhr  -ugh  October. 


BASINS  IN  COLORADO 

Dory  Hill  Basin  and  Mesa  Gulch  (Black  Hawk 
Qua  cl  r  a  ngTeT  a  r  e  two  fourth-order  basins  In  Secs. 

3 1,  33,  Twp.  2  S..  Ra.  72  \V.,  and  Secs.  4,  5,  B,  9, 
Twp.  3  S..  Ra.  72  \V.,  In  the  Colorado  Front  Range, 
Southern  Rocky  Mountains  province.  The  basins 
lie  in  a  dissected  pediment  or  other  old  erosion 
surface  of  moderate  relief,  sloping  eastward.  The 
immediate  relief  is  1,500  feet.  Higher  mountains 
near  the  Continental  Divide,  IP  to  15  miles  to  the 
west,  reach  heights  of  4,500  feet  above  the  basin 
level. 

The  rock  underlying  the  basins  is  granite-gneiss 
and  schist,  with  a  probable  general  northward  dip. 
About  60%  of  the  area  is  grassland  and  is  grazed, 
35%  to  37%  is  in  forest,  the  remainder  is  bare.  The 
estimated  P-E  index  is  76  to  80,  according  to  the 
nearest  records  taken  from  Frances  weather  sta¬ 
tion  at  nearly  the  same  elevation.  The  months 
March  through  July  receive  54%  of  the  mean  an¬ 
nual  precipitation,  though  rainfall  is  distributed 
rather  evenly  throughout  the  year. 

Cabrini  Gulch  (Morrison  Quadrangle)  is  a  small 
third-order  basin  below  the  Mother  Cabrini  Shrine, 
off  Highway  40,  in  the  southeast  corner  of  Sec.  16, 
Twp.  4  S.,  Ra.  70  W.,  and  is  tributary  to  Mt.  Ver¬ 
non  Canyon.  It  is  about  a  mile  west  of  the  moun¬ 
tain  front,  in  the  Colorado  Front  Range,  Southern 
Rocky  Mountains  province.  The  rock  underlying 
the  basin  is  hard,  black,  biotite-hornblende  schist, 
with  infrequent veinlets  of  feldspathic  material. 
About  95%  of  the  basin  is  grassland,  the  remainder 
is  covered  by  small  trees  and  bushes.  Grazing  is 
moderate  to  heavy,  and  the  range  is  fair.  The  main 
channel  of  the  basin  is  gullied  to  a  depth  of  20  feet 
near  the  mouth  and  less  than  one  foot  near  the  head. 
A  single  gully  in  the  slopes  was  caused  by  concen¬ 
tration  of  runoff  from  the  road  which  passes 
across  the  top  of  the  basin.  The  P-E  index  is  es¬ 
timated  at  46  on  the  basis  of  elevation  alone,  al¬ 
though  this  figure  may  be  high,  considering  that  a 
large  percentage  of  the  basin  faces  south.  The 
months  of  greatest  rainfall  are  April  through 
August,  when  61%  of  the  yearly  precipitation  occurs. 

Green  Mountain  Basin  (Morrison  Quadrangle)  is 
a  small  third-order  basin  on  the  southeast  end  of 
Green  Mountain,  in  the  center  of  Sec.  19,  Twp.  4S., 
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ratio  I’irtln.oni  m-riion,  Great  I'Uln.i  province 
aUiul  lw-o  n.ilc.i  rail  of  She  mount.. in  front  Under  - 
1  %  l»»K  the  tun  In  in  a  rturnr  IxniUlcry  rungi  imrralr 
Itlu  Green  .Mountain  conglomerate  i  of  alluvial  ori¬ 
gin.  ronniiillng  of  ijjiinwn  aiul  mctan.orphtc  rtwfc 
fragment*!.  Helow  tltr  conglomerate  m  a  partially 
consolidated,  greenish-brown  t  roMstieilcIril  naiul. 
Tim  mountain  In  a  remnant  of  Tertiary  alluvium, 
protected  tiy  the  res  inlaid  conglomerate.  An  ex¬ 
posure  of  granite  rohtiles  reveals  partial  disinte¬ 
gration  to  a  depth  of  two  and  one  naif  feel.  Almost 
al*  cobbles  have  a  caliche  coating  on  the  under 
sides. 

About  95%  of  the  area  is  grassland.  3%  to  4%  has 
small  trees  and  brush,  and  the  remainder  is  tiure. 
The  basin  Is  moderately  grazed,  and  the  range  Is 
fair  to  good.  The  main  channel,  atniut  8  feet  deep, 
is  trenched  and  has  ;>  trapezoidal  cross-section. 

The  smaller  channels,  higher  In  the  Ivisln,  are 
gullied  less  deeply  and  have  V-shaped  cross-sec¬ 
tions. 

The  P-E  Index  is  estimated  at  38  to  44  The  pre¬ 
cipitation  is  assumed  to  be  the  same  as  that  of 
Morrison,  3  miles  southwest,  where  61%  of  the  an¬ 
nual  precipitation  occurs  in  the  months  of  April 
through  August. 

Delonde  Gulch  (Nederland  Quadrangle)  is  a 
fourth-order  basin  in  Secs.  3  and  4,  Twp.  1  S., 

Ra.  73  W.,  in  the  Colorado  Front  Range.  The  im¬ 
mediate  area  has  about  1,800  feet  of  relief,  while 
six  or  seven  miles  west  the  peaks  of  the  Conti¬ 
nental  Divide  rise  to  heights  of  4,000  feet  above 
the  basin.  The  rock  in  the  area  is  granite-gneiss 
and  schist,  disordered,  but  with  a  general  north¬ 
ward  dip.  Considerable  mining  was  done  immed¬ 
iately  west  of  the  basin  until  the  late  1930’s, 
though  no  mines  were  seen  in  the  basin  itself.  More 
recently  logging  has  become  the  major  activity, 
but  :ts  eventual  extent  cannot  be  predicted.  Forest 
covers  nearly  100%  of  the  basin.  The  estimated 
P-E  index  is  84,  and  the  precipitation  is  heaviest 
in  the  months  of  March  through  August,  according 
to  the  nearest  weather  records,  at  Gold  Hill. 

Hicks  Gulch  (Nederland  Quadrangle)  is  two 
miles  south  of  Delonde  Gulch,  in  Secs.  15,  16,  17, 
Twp.  1  S.,  Ra.  73  W.,  at  about  the  same  elevation, 
and  is  very  similar  to  it,  except  as  described  be¬ 
low.  Many  inactive  mine  drifts  are  present  in 
Hicks  Gulch.  To  provide  timber  for  mines,  most  of 
the  large  limber  pines  in  the  upper  part  of  the 
basin  were  logged  35  to  50  years  ago.  The  basin  is 
grazed  lightly,  and  the  range  is  good  to  excellent 
in  the  grasslands.  Aerial  photographs  show  about 
75%  of  the  basin  covered  by  forest  and  bushes,  20% 
to  23%  by  grass,  and  the  remainder  bare.  There 
are  no  gullies,  but  rejuvenation  of  the  main  channel 
is  evident  by  the  steep-walled  inner  channel  which 


t*-i  oii.m  conn  tile  r  a  hi  y  deeper  toward  the  mouth  In 
each  of  the  two  main  I  r  Itula  r  len,  nmall  Inic-IlKi' 
featuren  extern!  along  the  (funnel  for  nearly  a  quar¬ 
ter  of  a  mile,  and  are  prutubly  the  renull  of  a  late 
I’lelnloeeive  earth- flow  or  noltfluctton  episode,  The 
lunin  wan  prolubly  never  glaciated,  since  no  spur 
trimming  or  other  alteration  of  the  valley,  and  no 
cirque  or  nlvalton  cirque  in  present  The  broad, 
poorly-drained  upland  in  the  highest  part  of  the 
lunin  was  affected  greatly  by  pertglacl.il  phenomena, 
such  as  frost  heaving,  and  even  now  soil  moisture 
may  In.*  frozen  a  few  feet  lie  low  the  surface  during 
much  of  the  year.  The  P- h  index  is  estimate  I  to  tie 
79  to  85,  depending  on  elevation  and  exposure. 

Moon  Gulch  (Nederland  Quadrangle)  is  in  Secs. 

2.  37  10,  IT  Twp.  2  S..  Ra.  TT  W.,  about  four  miles 
SSE  of  Hicks  Gulch,  al  alxwt  the  same  elevation. 
Underlying  the  basin  is  the  same  granite-gneiss 
and  schist  as  in  the  two  lias  ins  farther  north  in  this 
quadrangle,  and  the  planes  of  schistosity  have  a 
general  northward  dip.  The  greater  height  of  the 
southern  part  of  the  basin  is  due  to  the  presence  of 
a  fine-grained  granite  stock,  resistant  to  eiosion. 
Around  this  stock  are  numerous  inactive  mine  ex¬ 
cavations.  Considerable  portions  of  the  lias  in  have 
been  burned  over,  some  rather  recently,  and  are 
gullied  extensively.  About  50%  of  the  basin  area  is 
in  forest  cover,  35%  burned,  and  the  remainder  is 
grassland.  The  grasslands  are  primal  ily  moderately 
grazed  meadows  near  the  main  channel,  ami  the 
range  is  in  good  to  excellent  condition;  elsewhere 
the  range  is  poor,  and  not  grazed.  The  main  channel 
has  produced  a  broad  floodplain,  bounded  by  ter¬ 
races  of  20  to  30  feet  in  height.  The  floodplain  nar¬ 
rows  rapidly  upstream,  and  is  not  apparent  more 
than  a  mile  above  the  mouth. 

The  P-E  index  is  estimated  to  be  78  to  85,  de¬ 
pending  on  elevation  and  exposure.  According  to  the 
nearest  records  of  the  Silver  Lake  weather  station, 
the  months  with  greatest  precipitation  are  March 
tiirough  July,  when  54%  of  the  annual  amount  occurs. 

Tributaries  to  Saruche  Canyon  and  Gallinas 
Creek  (Starkville  Quadrangle)  comprise  two  basins, 
one  third-order  and  one  fourth-order,  in  the  Park 
Plateau,  Raton  section,  Great  Plains  province.  The 
tributary  to  Gallinas  Creek  lies  in  Secs.  27,  34, 

Twp.  34  S.,  Ra.  64  W.;  the  tributary  to  Saruche 
Canyon  lies  partially  in  the  southwest  quarter  of 
Sec.  32,  Twp.  34  S,,  Ra.  64  W.  The  rock  is  alternat¬ 
ing  thick-bedded  sandstone  and  thinner  shale,  and 
coal  strata.  The  beds  are  flat-lying,  and  cliffs  are 
common  on  the  sandstone.  Considerable  coal  min¬ 
ing  operations  are  carried  on  in  the  area,  though 
not  in  these  basins.  The  estimated  P-E  index  is  48. 
The  nearest  source  of  weather  records  is  Madrid, 

10  miles  northwest,  where  66%  of  the  annual  pre¬ 
cipitation  occurs  during  the  period  April  through 
August. 
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BASINS  IN  NEW  MEXICO 

CoUiiiiwixkI  Canyon  (Allir  Canyon  Qaadranglr) 
lu-ti  mainly  in  Sro>.  1,  2,  3.  Twp.  1C  S. .  Ra.  12  W., 
in  Hit  Mexican  Highland  tirclion,  Basin  ami  Range 
pruumc.  linn  area  wan  not  Included  In  the  I  If  Id 
nludy.  The  P- 1  index  In  38,  entiuialed  (ruin  l hr 
elevation.  Weather  data  are  av.ulahh  Iron,  the 
Mlmhrc.n  Hanger  Station,  one  mile  away.  The  (t  real  - 
e.nt  r.Unl.Ul  occurs  in  July  through  September, 
when  52%  of  the  annual  precipitation  occurs. 

Tributaries  to  Saplllo  Creek  (All If  Canyon 
Quadrangle!  are  two,  small,  fourth-order  luslns 
in  Secs.  2(F,  21,  29,  Twp.  15  S..  Ha.  12  W.,  In  the 
Mexican  Highland  section,  Basin  and  Range  prov¬ 
ince.  These  areas  wire  not  visited  for  field  study. 
The  estimated  P-E  index  Is  *10,  and  the  rainfall 
characteristics  arc  assumed  to  be  the  same  as 
above  for  the  Mlmbres  Ranger  Station,  six  miles  to 
the  southeast. 

Sawmill  Canyon  and  Tributaries  (Big  Burro 
Mountains  Quadrangle)  a  .e  threcTTourtli-order  ba¬ 
sins  in  Secs.  11,  TJ^Fa,  7  wp.  20  S.,  Ra.  16  W.,  and 
Secs.  6,  7,  18,  Twp.  15  W.,  Ha.  20  S..  in  the  Mexi¬ 
can  Highland  section,  Basin  and  Range  province. 
They  extend  over  a  large  portion  of  the  Big  Burro 
Mountains,  where  the  area  has  a  relief  of  2,300 
feet.  The  Big  Burro  Mountains  are  an  isolated 
granite  mass,  which  is  structurally  complex,  con¬ 
taining  many  dikes  of  fine  to  coarse  texture  rang¬ 
ing  in  composition  from  granitic  to  ultramafic. 
Wherever  exposed,  the  granite  is  deeply  and  thor¬ 
oughly  weathered.  Forest  of  pines  and  oak  brush 
covers  about  -10%  of  the  area,  the  remainder  has 
sparse  bunch  grass  or  is  completely  bare,  and 
the  range  is  poor.  Much  gullying  has  occurred  on 
the  slopes,  and  some  in  the  small  channels.  The 
estimated  P-E  index  is  40  to  55,  varying  with  ex¬ 
posure  and  elevation.  From  the  daily  weather  rec¬ 
ords  kept  for  seven  years  at  the  Smithsonian  Insti¬ 
tution  Astrophysical  Observatory  sunspot  labora¬ 
tory  (Tyrone  Station)  on  top  of  Jacks  Mountain, 
near  the  head  of  these  basins,  the  P-E  index  was 
54.  Greatest  precipitation  is  from  July  through 
September,  and  a  lesser  maximum  in  December 
through  February;  74%  of  the  vearly  precipitation 
occurs  during  those  months. 

Walnut  Creek  (Big  Burro  Mountains  Quadrangle) 
is  in  the  northeast  quarter  of  Sec.  7,  Twp.  20  S., 

Ra.  15  W.,  on  the  northeast  slope  of  the  Dig  Burro 
Mountains.  These  two  third-order  basins  differ 
|>'urii  Si^nulJ  (’ll  i.yjin  in  (heir  ig  nealer  la  live  re¬ 
lief  and  abundance  of  vegetation;  forest  covers 
about  85%  of  the  area,  the  remainder  is  bare.  The 
estimated  P-E  index  is  50  to  55.  In  other  respects 
these  basins  are  essentially  like  the  others  in  the 
Big  Burro  Mountains. 

Paso  del  Norte  Basin  (Bland  Quadrangle)  occu¬ 
pies  half  of  Sec.  2,  Twp.  ^8  N.,  Ra.  4  E.,  and  is 


near  the  IfcJwiHlary  Ixlinni  the  Southern  It  oc  k  y 
Mountalnn  province  and  the  Mexican  Highland  not- 
Hun  of  the  Bunin  and  ltai.ge  province,  it  is  on  the 
edge  of  the  Jeinez  Caldera,  and  lien  in  an  area  of 
*. olcanli  ilaclle,  tunall,  an  I  pumice  t*rd»,  dipping 
nouthc.inlv.ard.  Immediately  to  the  north  and  went, 
the  Sierra  de  Ion  Vallen  drops  700  feet  into  the 
Valle  Grande,  which  in  the  caldera  proper.  Else¬ 
where  along  the  Sierra,  peaks  reach  elevations  of 
10,200  feet.  Almost  100%  of  the  Ivistn  is  forest 
covered.  The  estimated  P-E  Index  is  98.  The  near¬ 
est  weather  station  Is  Los  Alamos.  July  and  August 
are  the  months  of  maximum  rainfall,  when  35£of 
the  yearly  precipitation  occurs. 

Peters  Dome  Basin  (Bland  C^udrangle)  Is  a 
fourth-order  tributary  oTCapulIn  Canyon  immedi¬ 
ately  east  of  Sec.  20,  Twp.  18  N.,  Ra.  5  E.,  In  the 
San  Miguel  Mountains,  in  the  Mexican  Highland 
section,  Basin  and  Range  province.  The  rocks  are 
nearly  horizontal,  alternating  basalt  and  pumice. 
Though  on  all  sides  of  the  basin  the  basalt  forms 
cliffs  800  to  900  feet  in  height,  the  basin  itself  is 
free  of  larger  cliffs;  only  near  its  mouth  has  it  be¬ 
gun  to  erode  into  the  main  cliff-forming  layer. 
Forest  and  bushes,  partially  burned  50  to  75  years 
ago,  cover  about  94%  of  the  basin;  4%  is  grassland, 
and  the  remainder  is  bare.  The  basin  is  poor  range 
and  is  not  grazed;  however,  it  does  serve  as  a 
main  stock  drive- way  for  herds  from  the  west.  No 
gullies  are  present.  The  estimated  P-E  index  is 
65  to  75,  depending  on  elevation  and  exposure.  The 
nearest  source  of  weather  records  of  the  same 
elevation  is  10  miles  northeast  at  Los  Alamos,  dis¬ 
tribution  of  precipitation  is  probably  the  same  as 
for  Paso  del  Norte  basin. 

Cerro  Pavo  Basin  (Cerro  del  Grant  Quadrangle) 
is  a  large  third-order  basin  in  Secs.  7,  18,  19,  Twp. 
21  N.,  Ra.  4  E.,  in  the  southernmost  end  of  the 
Southern  Rocky  Mountains  province.  The  basin  lies 
atop  a  high  lava-capped  mesa,  which  falls  away 
with  considerable  relief  to  the  west,  north,  and 
east,  and  rises  southward  to  the  smooth,  forested 
peaks  of  the  north  rim  of  the  Jemez  Caldera.  Un¬ 
derlying  the  basin  is  gray  to  pink  dacite  or  rhyolite 
porphyry,  with  occasional  interbedded  scoria,  all 
essentially  Hat-lying.  About  8  %  of  the  basin  is 
forest  covered;  the  remainder  is  grassland.  The 
range,  grazed  by  sheep  and  cattle,  is  fair  to  good 
in  the  meadows  near  the  main  channels.  No  gullies 
ale  present,  and  small  channels  are  poorly  defined. 
The  P-E  index  is  110,  estimated  from  the  elevation. 
The  dHKrttmtun  of  rtUbfuU  Its  to  Li  simi¬ 

lar  to  that  of  Capulin  Ftanger  Station,  where  the 
months  of  greatest  rainfall  are  May  through  Septem¬ 
ber,  when  74%  of  the  annual  precipitation  occurs. 

Basins  in  Hurley  East  Quadrangle  are  located 
in  Secs.  15,  22,  23,  Twp.  18  S.,  Ra.  12  W.,  north¬ 
east  of  the  mining  town  of  Hurley.  These  t.iree 
third-order  and  one  fourth- o"der  basins  were  noi 
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Vtnlti'O  til  I  lie  ( trlil.  V  Jnibllnlucl  gcolij^;  tc  map 
(Paige,  lOlo)  nlw»ii  llul  Ternary  lunall  ami  aiulr- 
nlle,  overlying  gravel  ami  lull,  lui  11.  ihe  higher 
part*  l(  lln  Bin  Inn.  Thu  P-E  index  lit  31  to  35  n- 
U m. i ted  it  to  elevation  only  Tin-  nearest  sou  .  e  ol 
rainfall  data,  Ft.  Bayard,  records  that  S4 X>  of  the 
yearly  prerlpitat ion  Is  in  tin-  montlui  July  through 
September. 

Hairpin  Basin  M.i  ro.ta  Quadrangle)  is  a  small 
tf’lrci-orJer  Bisln  In  Sees.  1,  2.  Tup.  21  N\,  Ha.  2 
K  .  and  Is  in  Hie  southernmost  part  of  the  Southern 
Hoik)  Mountains  province.  The  main  rock  type  is 
thlckly-lxnlded  red  sandstone,  essentially  hori¬ 
zontal,  i  roppm^  out  In  low  cliffs  on  south-facing 
slopes.  About  90%  of  the  area  Is  covered  by  brush 
and  forest,  5%  is  grassland,  and  5%  Is  Bire  soil 
and  rock.  The  basin  Is  moderately  grazed,  and  the 
range  on  north-facing  slopes  is  good  to  fair.  The 
only  gullies  present  are  those  formed  as  a  result 
of  recent  road  building.  The  P- E  index  is  72,  es¬ 
timated  from  the  elevation.  Capulin  Ranger  Station 
is  the  nearest  source  of  weather  data;  there,  74% 
of  the  annual  precipitation  occurs  in  May  through 
September. 

Poleo  Dasln  (Jarosa  Quadrangle)  is  a  fourth- 
order  basin  in  Secs.  27,  28,  Twp.  22  N.,  Ra.  2  E., 
in  the  southern,  ost  part  of  the  Rocky  Mountains 
province,  near  the  northern  limit  of  the  plateau 
comprising  the  San  Pedro  Mountains.  The  topograpy 
drops  rapidly  northward  into  the  nio  Puerco  Valley. 
The  rock  types  are  reddish,  coarse-to-conglomer- 
atic  arkosic  sandstone,  gray  quartzose  sandstone, 
crinoidal  limestone,  and  shale.  In  the  greater  part 
of  the  area,  the  beds  are  flat-lying,  but  dip  east¬ 
ward  10°  to  20°  in  the  higher  western  part  of  the 
basin.  Trees  and  brush  cover  an  estimated  96%  of 
the  area,  grass  1%,  and  the  remainder  is  bare  soil 
and  rock.  Grazing  is  light  to  moderate;  the  range  is 
poor.  Gullies  are  present  in  the  main  channel  and 
one  of  the  main  tributaries.  They  are  V-shaped, 
one  to  eight  feet  deep,  cut  in  alluvium,  but  not  bed¬ 
rock.  The  estimated  P-E  index  is  77,  and  rainfall 
distribution  is  assumed  to  be  the  same  as  for  Capu¬ 
lin  Ranger  Station. 

Quartz  Hill  Basin  (Paxton  Springs  Quadrangle) 
is  a  fifth-order  basin  in  Secs.  19,  20,  29,  30,  Twp. 

10  N.,  Ra.  11  W.,  on  the  eastern  flank  of  the  Zuni 
Mountains,  in  the  Dalil  section  of  the  Colorado 
Plateau  province.  The  rock  underlying  the  entire 
basin  is  a  gneissic  granite,  with  occasional  sept¬ 
like  layers  of  muscovite-chlorite  schist,  dipping 
steeply  southwestward.  Of  several  mines  present, 
only  one  is  still  active.  The  granite  is  deeply 
weathered  in  most  places,  especially  at  the  eastern 
end  of  the  basin  where  remnants  of  arkose  origin¬ 
ally  deposited  directly  on  the  granite  are  found. 

The  greater  decomposition  here  can  be  attributed 
to  pre-arkose  weathering.  Trees  and  bushes  cover 
an  estimated  80%  of  the  basin,  grass  15%,  and  5% 


in  lure  noil  and  rock.  Tin-  area  in  grazed  moder¬ 
ately  lo  heavily,  range  In  poor  lo  fair.  The  few 
gullied  present  are  in  the  nmallent  ilia  line  In.  The 
P-E  Hide*  in  44.  estimated  from  the  elevation.  The 
ills trtbut ion  of  rainfall  Is  assumed  to  lx-  similar  to 
tliat  of  El  Morro  National  Monument.  18  miles  west, 
and  San  Rafael,  15  miles  northeast.  The  months  of 
greatest  precipitation  are  July  through  September, 
when  -18 %)  of  the  annual  rainfall  occurs. 

Camtnn  Basin  (Santa  Fe  C^i.idrangle )  is  a  small 
f ou rth-order  iuHln  located  in  the  northeast  quarter 
of  Sec.  28.  Twp.  17  N..  Ra.  9  E.,  in  the  Mexican 
Highland  section  of  the  Basin  and  Range  province, 
in  the  valley  of  the  Santa  Fe  River.  The  lithology  Is 
the  Santa  he  formation,  indurated  sand  and  gravel  of 
alluvial  origin.  About  15%  of  the  Basin  Is  covered  by 
trees  and  bushes.  30%  by  bunch  grass,  and  55%  Is 
Ixtre.  Grazing  Is  very  heavy,  mainly  by  sheep  and 
goats,  resulting  in  very  poor  range  conditions.  The 
main  channels  are  gullied,  with  vertical  sides.  The 
P-E  index  is  between  that  for  Santa  Fe  Airport  and 
Santa  Fe  City,  or  about  23.  The  distribution  of 
rainfall  is  isumed  to  be  about  the  same  as  that  of 
Santa  Fc  City,  where  44%  of  the  annual  precipitation 
is  in  July  through  September. 

Tano  Basin  (Santa  Fe  Quadrangle)  is  a  fourth- 
order  basin  essentially  like  Camino  Basin  in  lithol¬ 
ogy.  It  lies  north  of  the  city  of  Santa  Fe  proper,  but 
within  the  city  limits,  and  is  a  tributary  to  Arroyo 
de  los  Frijoles.  Trees  and  bushes  cover  80%  of  the 
basin  area;  grass  4%,  and  16%  is  bare.  The  range  is 
very  poor  to  fair,  and  has  probably  been  grazed  in 
the  recent  past,  though  not  immediately  prior  to  the 
time  of  field  study.  Most  of  the  channels  are  trenched, 
and  gullies  have  become  established  on  the  slopes. 

The  main  channel  is  deepened  to  15  feet  near  its 
head,  but  is  not  gullied  toward  the  mouth;  there,  a 
broad,  sandy  channel  is  present.  All  gully  sides  are 
nearly  vertical;  the  bottoms  are  usually  flat.  The 
P-E  index  is  estimated  at  33,  somewhat  above  that 
of  Santa  Fe.  The  distribution  of  precipitation  is  the 
same  as  for  Camino  Basin. 

Cedro  Peak  Basin  (Sedillo  Quadrangle)  is  a 
fourth-order  basin  in  Secs.  25,  36,  Twp.  10N.,  Ra. 

5  E.,  in  the  Mexican  Highland  section,  Basin  and 
Range  province,  in  the  northernmost  extension  of 
the  Man^ano  Mountains.  The  rock  is  the  Pennsyl¬ 
vanian  Magdalena  formation,  composed  of  limestone, 
feldspathic  sandstone,  breccia,  and  thin  shale  beds. 

An  anticlinal  axis  passes  north -south  through  the 
western  part  of  the  basin,  from  which  beds  dip  gently 
northeast  and  northwest.  About  97%  of  the  basin  is 
forest  covered,  2%  is  grassland,  and  the  remainder 
is  bare.  Grazing  is  light  and  the  range  is  generally 
very  poor,  with  only  small  areas  of  fair  range.  The 
main  channel  is  gullied  near  its  mouth  to  a  depth 
of  five  feet.  All  other  channels  show  signs  of  intense 
recent  scour,  down  to,  but  not  into  bedrock,  but  are 
not  actually  gullied.  Slopes  give  evidence  of  recent 


20 


rilit'ct-uMiili  nr  rill  mi* ion,  Inj t  luvr  v»-rv  few  gul- 
I it**.  A  single  gully  has  l**comr  established  down  4 
rulgr  r rent  e r<>r> ton  having  occurred  along  a  trail. 
The  1’-  K  Index  is  estimated  to  U-  30  to  -10,  depend¬ 
ing  on  elevation  and  exposure.  The  dlstrllmtton  of 
rainfall  is  assumed  to  l>e  the  same  as  lor  the 
Tljeras  Hanger  Station,  three  miles  northwest, 
where  during  the  months  of  May.  .July,  August,  and 
September  -17  (,  of  the  annual  precipitation  occurs. 

Tributaries  to  Hear  Canyon  (Twin  Sisters  Qiad- 
r.angle)  are  three  fourth-order  ‘co  ins  in  Sees. "13, 

14,  33,  2*1,  Tup.  10  S.t  Ha.  1  ,  In  the  Mexican 

Highland  section,  Basin  and  1...  ge  province.  The 
area  is  a  rugged,  dissected  section  of  lavas,  mainly 
andesite  and  tvisnlt  (Paige,  1 9 1 G ) ,  alternating  with 
gravels  and  luffs;  the  lava  beds  form  prominent 
steps  in  the  hillside.  The  area  was  not  visited  in 
the  field.  The  estimated  P-E  index  is  GO.  The  near¬ 
est  source  of  weather  records  is  Pinos  Altos, 
where  50%  of  the  yearly  precipitation  occurs  during 
July,  August,  and  September. 

Little  Cherry  Creek  (Twin  Sisters  C^iadrangle) 
is  a  fifth-order  basin  in  Secs.  15,  16.  19,  20,  21,  22, 
Twp.  16  S.,  Ra.  13,  W.,  in  the  Mexican  Highland 
section  of  the  Basin  and  Range  province.  The  esti¬ 
mated  P-E  index  is  51.  In  other  respects,  the  basin 
is  similar  to  the  others  in  Twin  Sisters  quadrangle. 
The  liasin  was  not  inspected  in  the  field. 

Mill  Creek  (Twin  Sisters  Quadrangle)  is  a  fourth- 
order  basin  similar  to  the  two  preceding  ones,  ex¬ 
cepting  the  P-E  index,  which  is  estimated  at  45. 

Wilson  Crpek  Basins  (Twin  Sisters  Quadrangle) 
are  four  fourth-order  basins  similar  to  the  three 
other  areas  in  Twin  Sisters  quadrangle,  except  that 
the  estimated  P-E  index  is  51.  Forest  covers  100% 
of  the  basin. 


BASINS  IN  UTAH 

Dry  Hollow  (Lehi  Quadrangle)  is  a  fourth-order 
basin  in  Secs.  20,  21,  Twp.  4  S.,  Ra.  1  E.,  on  the 
south  side  of  the  Traverse  Mountains,  five  miles 
from  the  foot  of  the  Wasatch  Front,  in  the  Great 
Basin  section,  Basin  and  Range  province.  The  Tra¬ 
verse  Mountains  have  not  been  affected  by  Recent 
uplift,  as  have  the  Wasatch  Mountains,  resulting  in 
less  rugged  topography  (Marsell,  1953,  p.  20).  The 
prevailing  rock  type  is  hard  buff  quartzite.  The  soil 
contains  considerable  clay  that  is  probably  wea¬ 
thered  from  limestone  and/or  lava,  though  these 
lithologies  are  not  seen  in  the  basin.  Trees  and 
brush  cover  60%  of  the  area,  grass  35%,  and  5%  is 
bare.  The  range  which  is  generally  poor,  is  heavily 
over-grazed  in  places  by  sheep  and  cattle.  While 
a  few  small  tributaries  are  gullied,  the  main  chan¬ 
nel  is  not  and  is  indistinct  in  most  places.  The 
tributaries  usually  end  in  alluvial  cones,  and  do  not 
direc  ly  join  the  main  channel.  The  entire  surface 


uf  the  l us  in  in  cum* rod  by  quartzite  fragments, 
which  caurir  mutability  uf  the  slope*  on  the  steep 
hillsides.  The  estimated  R- F  index  is  40  to  52.  de¬ 
pending  on  elevation.  The  distribution  of  precipita¬ 
tion  is  protiibly  nearly  the  same  as  that  of  Alpine, 
five  miles  east,  where  except  for  four  dr>  months, 
June  through  Septemlx*r,  ram  is  distributed  evenly 
throughout  the  year.  Only  21%  of  the  annual  precipi¬ 
tation  occurs  in  those  months. 

Hog  Hollow  (Lehi  Quadrangle)  is  a  fifth-order 
basin  similar  to  Dry  Hollow,  but  Is  nearer  the 
Wasatch  Front.  The  surface  of  the  basin  shows  con¬ 
siderable  andesitic  lava  and  much  less  quartzite. 

A  greater  percentage  of  this  area  is  covered  by 
trees  and  bushes  Hum  in  Dry  Hollow.  The  estimated 
P-E  index  is  51  to  55,  depending  on  dotation.  The 
distribution  of  rainfall  is  the  same  as  for  Alpine, 
which  Is  1  1/2  miles  southeast. 

Maple  Hollow  and  Mcrrer  Hollow  (Lehi  Quad¬ 
rangle)  are  two  (ourth-ordtTr  basins  located  in 
Secs.  15  and  16,  21,  respectively,  Twp.  4  S.,  Ra.  1 
E.,  and  are  similar  to  the  two  preceding  areas  in 
this  quadrangle,  with  a  greater  amount  of  lava 
cropping  out  than  in  Dry  Hollow.  The  estimated 
P-E  index  is  50  to  53,  depending  on  elevation.  The 
distribution  of  precipitation  is  assumed  to  be  the 
same  as  above  for  Alpine. 

Pole  Canyon  (Spanish  Fork  Peak  Quadrangle)  is 
a  large  Tifth-order  liasin  in  Secs.  13,  14,  15,  22, 

23,  24,  27,  Twp.  9  S.,  Ra.  3  E.,  in  the  Wasatch 
Mountains  of  the  Middle  Rocky  Mountains  province. 
Almost  continual  uplift  during  the  late  Pleistocene 
and  Recent  has  caused  rejuvenation,  and  accounts 
for  the  great  ruggedness  of  the  area.  The  main  lith¬ 
ologies  are  limestone  and  quartzite.  Forest  covers 
almost  the  entire  area.  The  estimated  P-E  index  is 
70  to  80,  depending  on  elevation.  The  distribution 
of  precipitation  is  assumed  to  be  similar  to  that  of 
Thistle,  where  the  rainfall  is  fairly  evenly  distrib¬ 
uted  throughout  the  year,  with  slight  maxima  in  Jan¬ 
uary  and  March,  and  a  minimum  in  June  and  July. 

Deadmans  Hollow  (Springville  Quadrangle)  is  a 
small  third-order  basin  in  Secs.  1  and  36,  Twps.  7 
and  8  S.,  Ra.  3  E.,  in  the  Wasatch  Mountains  of  the 
Middle  Rocky  Mountains  province.  Uplift  has  caused 
rejuvenation  and  rugged  topography.  The  main 
rock  type  is  buff  quartzite  and  sandstone.  Trees 
and  brush  cover  about  20%  of  the  area,  the  remain¬ 
der  has  sparse  bunch  grass  or  is  completely  bare. 
The  estimated  P-E  index  is  45  to  50,  depending  on 
elevation.  The  distribution  of  rainfall  is  probably 
the  same  as  for  Provo  Bench  station,  where  the 
occurrence  of  precipitation  is  fairly  even  through¬ 
out  the  year,  with  a  slight  maximum  in  May,  and  a 
minimum  in  June  and  July. 

Spring  Creek  (Springville  Quadrangle)  is  a  fourth- 
order  basin  in  Secs.  24,  25,  26,  Twp.  7  S.,  Ra.  3  E., 
in  the  Wasatch  Mountains  of  the  Middle  Rocky 
Mountains  province.  Uplift  of  the  mountain  range 
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Iuh  caused  rejuvenation  and  r united  Ic^mv: raphy  in 
l hr  area  The  main  rock  lyj>e  is  i;ray  limestone 
4iuJ  hull  quarlziti  T rern  and  (rushes  I'uo-r  BH  of 
flu  lusin,  (.Tarts  nearly  151,  while  only  .*  ni-^ll^t- 
hie  pari  ts  Uire.  Caltle  ^ra/.e  flu-  twain  heavily 
near  I  he  nwda  aiKl  the  main  channel,  the  rantfe  la 
poor  to  (air.  The  main  dunnel  la  gullied  to  a  depth 
of  3  feel  near  ll.a  head  and  20  feel  at  the  mouth. 

The  P-h  Index  la  estimated  at  70  to  75,  depending 
on  elevation.  Til'1  distr lliutlon  of  rainfall  is  prob¬ 
ably  the  same  as  above  for  Provo  Dench  station. 


fribwtary  to  Hobble  Creek  <Sj>r  ir^v  llle  Cjaat- 
rankle*  in  a  nmail  third-order  twain  in  Seta  3C 
ami  1,  T'*pn,  7  S.  and  ti  ttv  3  K  The  Iwnin  la 
similar  to  (Vadmann  Hollo'*.  I  he  main  lithology  in 
tniff  sandstone,  quartzite,  a nd  (•  r  a y  limestone  The 
twain  is  grazed  heavily,  and  the  ran^e  is  poor  to 
fair.  A  few  new  gullies  are  present  on  the  uistn 
slopes,  tail  old  clunmds  are  not  £uUled.  The  esti¬ 
mated  P-  K  index  is  55  to  60.  Precipitation  dtstrilAi- 
tton  la  the  same  as  for  Provo  Bench  station,  the 
nearest  source  of  weather  data. 


Methods  of  Statistical  Analysis 


ABPKOPMATKNESS  OF  STATISTICAL  METHODS 

The  important  assumptions  Inherent  in  statisti¬ 
cal  procedures  are  becoming  well-known.  Depart¬ 
ures  in  practice  iron,  r l|» itl  adherence  to  the  condi¬ 
tions  imposed  by  those  assumptions  will  produce 
more-or-less  serious  effects  on  validity  of  conclu¬ 
sions  drawn  from  the  results  of  a  test.  In  this  study, 
the  attempt  was  made  to  comply  with  sound  statis¬ 
tical  methods,  but  some  wide  departures  exist. 

Normality  of  parent  populations  is  fundamental 
in  the  parametric  methods  of  analysis  used  here. 

The  number  and  size  of  samples  taken  in  any  single 
catego.y  of  measurement  within  a  given  environ¬ 
ment  are  insufficient  to  allow  detection  of  moderate 
departures  from  normality.  One  exception  is  max¬ 
imum  valley-side  slope,  which  other  investigators 
have  found  to  be  normally  distributed  within  a 
drainage  basin  (Strahler,  1950,  p.  681-685);  this  is 
also  true  of  all  areas  studied  here.  Another  ex¬ 
ception  is  the  drainage  density  of  subordinate  ba¬ 
sins  within  a  single  larger  basin,  or  within  the 
same  rock  type;  these  distributions  do  not  appear 
to  be  significantly  different  from  normal,  although 
total  samples  are  small.  Miller  (1953,  p.  14)  ren¬ 
dered  symmetrical  the  distribution  of  areas  of 
first-order  basins  by  logarithmic  transformation. 

Sample  means  tend  to  be  normally  distributed, 
regardless  of  the  parent  distribution,  when  samples 
are  large,  and  the  parent  populations  have  finite 
variances  (Dixon  and  Massey,  1951,  p.  52).  Most 
data  used  here  are  sample  means,  and  for  tnat  rea¬ 
son  are  expected  to  be  approximately  normally  dis¬ 
tributed.  The  analysis  of  variance  and  similar  tests, 
are  remarkably  robust  to  departures  from  normal¬ 
ity  (Box,  1953,  p.  318)  for  tests  on  means.  The 
question  of  normality  is  therefore  probably  of  less 
importance,  for  the  purposes  of  this  study,  than 
other  factors. 

Objective  sampling  is  the  sine  qua  non  of  ail 
statistical  methodsjThis  requirement  of  statistics 
imposes  restrictions  on  the  methods  of  obtaining 
data  suitable  for  analysis.  Choice  of  sampling 
methods  is  the  area  of  greatest  divergence  irom 
conventional  methods  used  in  geology.  Practicing 
statisticians  have  long  realized  the  difficulty  in  ob¬ 
taining  unbUsed  samples  of  a  population  when  there 
is  freedom  of  choice  in  selecting  the  observations 
to  be  recorded.  The  training  of  a  geologist  is  at 
least  in  part  that  of  learning  to  observe  in  a  man¬ 
ner  that  assumes  a  bias  toward  representation  of 


the  true  .situation,  thus  reducing  variability  In 
measurements  and  allowing  t realm  nt  by  non- 
slattsltral  methods.  In  the  statistical  approach,  It 
must  lx.-  assumed  that  the  true  char  .cterlstics  of  the 
universe  are  unknown  and  can  lx1  estimated 
from  random  samples.  Ain  nipls  to  obtain  data  close 
to  some  value  arrived  at  by  “judgment,’  “insight,’ 
or  “experience."  will  nearly  always  result  in 
biased  samples,  and  failure  of  the  method  of  analy¬ 
sis.  In  this  study,  samples  were  taken  Ijoth  by  ran¬ 
dom  and  non- random  methods,  and  results  an  .u  - 
cordingly  subject  to  qualification. 

Study  areas  wo,  chosen  after  cursory  examina¬ 
tion  showed  then  to  be  satisfactory  in  accessibil¬ 
ity,  absence  of  cliffs,  low  In  roughness,  and  in  the 
mature  stage  of  development.  No  effort  was  made 
to  select  by  random  methods  sample  basins  in  all 
types  of  rock;  the  majority  are  in  granite  and  coarse 
clastic  sediments,  and  few  in  shale  or  limestone. 
Measurements  within  basins  were  largely  on  areas 
selected  to  be  representative  and,  as  noted  „uie, 
certain  measurements  would  mean  little  on  areas 
differing  greatly  from  the  average  of  the  entire 
basin.  Selecting  measurement  sites  in  a  manner  to 
reduce  the  variance  of  the  observations  may  not 
bias  the  observed  mean.  However,  if  analysis  of 
variance  is  used  to  test  for  significance  of  differ¬ 
ences  among  means  of  a  measurement  in  several 
basins,  too  many  significant  results  will  be  obtained 
if  inferences  are  made  about  the  entire  population 
instead  of  the  portion  that  was  measured.  This  is 
because  the  error  attributable  to  a  single  measure¬ 
ment  will  occur  in  the  denominator  of  the  F  ratio, 
and  will  cause  it  to  exceed  the  critical  value  too 
frequently.  In  regression  analysis,  selection  in  in¬ 
dependent  variables  probably  does  not  seriously 
affect  the  value  of  correlation  and  regression  co¬ 
efficient  unless  either  extreme  or  intermediate 
values  are  eliminated.  In  this  case,  as  extreme 
values  were  frequently  rejected,  the  indicated  cor¬ 
relation  will  be  too  low.  Selection  in  the  dependent 
variables  has  greater  effect.  If  the  selection  tends 
to  reduce  the  variability  of  the  dependent  variable, 
the  indicated  correlation  will  be  too  high  (Ezekiel, 
1930,  p.  266-267). 

Additivity  and  linearity  of  effects  is  an  impor¬ 
tant  assumption  in  multiple  regression  analysis, 
when  several  variables  are  examined  for  their  in¬ 
fluence  on  another.  The  plane  of  best  fit  through  a 
set  of  points  can  always  be  found,  and  in  cases  pre¬ 
sented  here  it  is  used  only  where  the  dependent 
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-.  4.'U1I<  is  1UH-4  f  '»Svr»i  plottrl  4*f4ifiat  i*4ch  iivitc- 
(irrulrlil  ,  4ft4t>l<-  *  It  is  <w>!  itelWilti'ly  klvonii  »lu  Itirf 
.jf  tvoi  rlJrfla  of  ro»  k  ly|»o,  1c£rt4l  foirr,  4f*l  i' 1 1  - 
liulii'  rlriiidila  oft  nit)l|ihulii^|('  |iriij>rfltra  4fr  4i!dl- 
l  i  *.  r.  Hours  cf,  4  s  i»  ill  tw  srrli  lit  rolllt  Idr  rlU£  rffci'lJi 
of  I l»c nt  (4rl<ifn  on  \4llrV- nUtr  tilojx-  .iikI  dr4UV4£i‘ 
<t<-ll*|l),  IU»l»*  I  tfU-4  f  Uo-l4lUl)|('  |)|olfi  oti<4  I  fled 
■* r r r  rri  ...let!  by  t4ktftg  Icagurlllumi.  indicating 
either  multtpllr.it  |\ r  rllrrln  ur  non -linear  variables. 

Mulllplr  n-grertMon  ami  correlation  ran  In-  used 
•ahere  Y  in  a  dependent  random  variable  and  X  j , 

X  ,  ...  X„are  independent  random  variables  whose 
sallies  are  known,  ami  not  subject  to  errors  In 
measurement  or  observation.  In  the  cases  to  lie 
considered  here,  va'ues  of  all  variables  were  ob¬ 
tained  by  measurement  in  the  field,  from  maps,  or 
by  meterob  gleal  observation,  and  hence  ai<>  esti¬ 
mates  subject  to  at  least  some  observational  error. 
The  effect  of  introducing  observational  error  into 
the  independent  variables  Is  to  reduce  the  indicated 
correlation  and  regression  coefficient.  The  regres¬ 
sion  Is  Y  on  estimated  Xj,  ...  X„,  rather  than  on 
true  Xj,  ...  X„  (Ezekiel,  1930,  p.  225;  Professor 
Howard  I.evene,  Personal  Communication).  The 
combineu  influence  of  selection,  tending  to  increase 
indicated  correlation  and  regression  coefficients, 
ard  introduction  of  error  into  the  X,,  tending  to 
reduce  the  coefficients,  is  complex,  and  the  effects 
on  the  conclusions  drawn  in  this  study  are  not  known. 

The  theoretical  considerations  discussed  in  the 
previous  paragraphs  are  of  great  importance  in 
confirming  hypotheses  by  statistical  methods,  when 
research  and  experimentation  have  been  properly 
controlled,  when  properties  of  experimental  mate¬ 
rial  are  well-known,  and  when  only  a  few  variables 
are  considered  at  one  time.  The  research  in  this 
study  is  primarily  investigative,  and  it  is  the  wri¬ 
ter’s  opinion  that  it  is  probably  sound  scientific 
practice  to  lean  toward  declaring  too  many  rela¬ 
tions  or  differences  to  be  significant,  rather  than 
risk  failing  to  detect  a  truly  significant  relation.  It 
should  be  understood  that  the  hypotheses  thus  pro¬ 
posed  are  not  considered  to  be  definitely  confirmed 
until  rigorous  investigation  has  produced  sustaining 
evidence. 
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The  Information  available  for  each  basin  was  re¬ 
corded  on  two  series  of  punch  cards.  Fourth  and 
fifth-order  basins  were  divided  into  third-order 
basins,  with  one  card  for  each.  Cards  were  then 
easily  sorted  into  classes  to  obtain  frequencies  in 
each  cell  in  a  two-way  table.  Usual  methods  of  an¬ 
alysis  of  grouped  data  yielded  correlation  and 
regression  coefficients  with  comparative  speed. 

Of  the  total  number  of  properties  investigated, 
those  proving  to  be  significantly  correlated  were 
considered  for  further  analysis. 

Correlation  analysis  of  grouped  data  is  limited 
by  several  considerations:  (1)  because  of  grouping, 
the  estimated  standard  errors  tend  to  be  too  large, 
and  the  correlation  coefficient  too  small.  It  is  nec¬ 
essary  to  have  a  considerable  number  of  cells  in 
the  correlation  table,  and  this  requirement  is  lim¬ 
ited  by  the  number  of  holes  in  the  cards.  (2)  The 
correlation  table  is  difficult  to  interpret  as  to 
linearity,  degree  of  correlation,  importance  of  ex¬ 
treme  values,  etc.  (Croxton,  1953,  p.  136).  Despite 
these  objections,  the  method  of  correlation  using 
grouped  data  on  punched  cards  proved  valuable  in 
facilitating  rapid  analysis.  Appendix  V  gives  tne 
correlation  tables  and  the  correlation  coefficients 
obtained,  for  a  number  of  important  relations. 


REGRESSION  AND  CORRELATION  ANALYSIS  OF 
PAIRED  VARIABLES 


FINAL  MULTIVARIATE  ANALYSIS 


Preliminary  analysis  of  grouped  data:  The  data 
of  greatest  interest  and  promise  were  grouped  into 

♦In  view  of  the  differences  in  opinion  prevalent  on  the  ef¬ 
fects,  say,  of  increased  rainfall  on  texture  of  topography, 
when  modifications  of  intensity-frequency  relations,  soil 
properties,  vegetal  cover,  and  stream  regimen  are  con¬ 
sidered,  the  question  of  non-additivity  of  effects  should 
be  the  subject  of  further  research. 


Certain  of  the  significant  relations  discovered 
from  analysi:  f  grouped  data  were  further  ana¬ 
lyzed  by  ungrouped  methods,  in  order  to  find  the 
statement  of  most  linear  function  giving  the  best 
fit  and  lowest  scatter.  The  final  results  of  these 
are  given  in  part  B  of  Appendix  V.  These  relations 

♦Keysort  cards,  manufactured  by  The  McBee  Company, 
Athens,  Ohio. 


-» r t  v  .ilbl  4 *»  lliry  4i  r,  ln.it  thru  ui«'(uliu*«r>  t  .in  Ik* 
tnrr«4r»rt|  if  combined  wttn  oilier  me.u.uri ii.enb. 
Js*-viT4l  combination:*  of  loga r tlbm lc  ami  arithmetic 
m  4let)  Were  tried  for  each  pair  ol  variables. 

Multivariate  i  nr  relation  analysis  1.4  4  method  ol 
4 nal yt inn  the  effect  uf  several  independent  variables 
on  ihe  van.ltlon  o(  4  dependent  variable.  Tilth  llietll- 
ikI  i.s  necessary  In  determining  whether  each  In¬ 
dependent  variable  ads  independently  ol  the  others, 
nr  i(  an  apparent  significant  correlation  obtained 
on  a  two-variable  analysis  is  actually  due  to  their 
txilli  Ik- 1  tut  correlated  with  a  third  variable.  Where 
all  aspects  of  an  experiment  are  not  easily  con¬ 
trolled  and  all  possible  combinations  of  factors  are 
not  found,  mult. pie  correlation  and  regression  anal¬ 
ysis  will  Ik*  of  greater  use  in  quantitative  geomor¬ 
phology,  than  the  standard  experimental  designs 
employing  the  analysts  of  variance.  (See  analysis 
of  covarlanc'*.) 

The  methods  used  m  analyzing  two  to  four  vari¬ 
ables  are  usually  presented  fully  in  elementary 
statistics  textbooks.  The  amount  of  calculation 
needed  to  prepare  a  multiple  correlation  study  in¬ 
creases  as  some  power  of  the  number  of  variables 
considered.  In  this  study,  the  final  analysis  con¬ 
sists  of  determination  of  the  effects  on  a  dependent 
variable  (drainage  density),  of  five  independent 
variables  (infiltration  capacity,  P-E  index,  etc.). 

No  explicit  solution  of  normal  equations  for  this 
many  variables  is  available.  The  normal  equations 
were  first  expressed  in  matrix  form,  the  inverse 
matrix  found  by  the  Doolittle  method,  and  the  in¬ 
verse  then  used  10  obtain  the  desired  solutions. 

(See  Anderson  and  Bancroft,  1952,  p.  168,  and 
Ezekiel,  1950,  pp.  168,  363,  367).  The  full  solution 
is  given  in  Appendix  VII. 

Correlation  analysis  of  six  variables,  to  yield 
stable  solutions,  usually  requires  more  values  on 
each  variable  than  are  used  here.  In  this  study, 
however,  each  of  the  23  values  used  represents 
more  actual  measurements  averaged  to  obtain  it. 
Drainage  density  represents  the  average  stream 
length  per  unit  of  area,  and  was  obtained  from  ten 
to  300  individual  measurements  of  stream  length; 
each  P-E  index  is  obtained  from  12  paired  values; 
infiltration  capacity  is  an  average  of  four  to  11 
readings;  percent  barp  area  is  the  average  of  two 


b>  20  measurements,  roughness  HuiiiU-r  in  meas¬ 
ured  Inm.  hundreds  of  rock  fragments  Mm  *i  l«  ll<  r 
results  would  protiahly  Ik*  obtained  by  using  the  or 
igie.al,  uiiav e  raged  figures,  as  the  gain  in  degrees  *f 
freedom  would  I*'  c oiiside  1  able  in  n.oril  i.ises  L'n- 
foriunalely,  lo  do  this  w  uld  mean  duplicating  values 
of  the  dependent  variable  Therefore  Ihe  value  of 
llu*  drainage  density  would  have  to  Ik*  repeated  4 *> 
many  times  as  there  are  infiltration  readings  in  th* 
las  In.  and  there  would  lx*  no  real  gain  in  degrees 
of  freedom.  A  solution  would  Im*  the  use  of  analysis 
of  variance  with  unequal  replication  in  each  cell 
Init  this  is  unfeasible  because  all  con  binatlons  of 
the  five  independent  variables  are  not  found. 

Analysis  of  variance  is  a  technique  of  determin¬ 
ing  the  significance  of  difference  among  several 
means  drawn  from  distinct  environments.  This  study 
treats  the  effect  of  rock  type  on  several  morpho¬ 
metric  and  soil-cover  properties  uf  the  basins  in 
each  type.  A  simple  one-way  classification  is  used, 
and  six  rock  categories  are  distinguished.  Where  sig¬ 
nificant  differences  are  found  to  exist  among  the 
means,  Scheffc ‘s  method  of  judging  contrasts  is  used 
to  test  specific  pairs  or  groups  of  means  for  signifi¬ 
cant  differences  (Scheffc,  1953,  p.  87).  The  analysis 
of  variance  tables  are  presented  in  Appendix  VI, 
and  contrasts  are  analyzed,  wherp  :  apropriate. 

The  analysis  of  covariance  is  essentially  similar 
to  the  analysis  of  variance,  with  the  addition  that 
an  associated  measurement  is  made  for  each  item, 
in  order  to  determine  how  much  of  the  differences 
in  the  means  of  the  primary  measurements  is  the 
result  of  variation  in  the  related  variable.  The 
means  are  adjusted  to  what  they  would  be  were  the 
associated  measure  the  same  for  each.  Then,  if  sig- 
nigicant  differences  are  present,  these  are  due  to 
the  differences  in  the  main  variable  of  classification 
(Dixon  and  Massey.  1951,  p.  173).  One  aim  of  the 
study  was  to  determine  if  a  consistent  difference  in 
drainage  density  could  be  shown  to  exist  among  ba¬ 
sins  with  gullied  slopes,  trenched  channels,  and  no 
gullies.  Infiltration  capacity  was  the  associated 
measurement,  and  variations  in  drainage  density 
due  to  differences  in  infiltration  capacity  were  re¬ 
moved.  The  analysis  of  covariance  table  is  given  in 
Appendix  VI  B,  and  the  results  discussed  in  t.ie  con¬ 
clusions  section. 


Inferences  Drawn  From  Results 
of  Statistical  Tests 


RESTRICTIONS  ON  IN FFRFNCES 

Certain  rest  irtions  must  lx-  placed  on  tin  con¬ 
clusions  Ixmansr  of  the  nature*  of  tin*  statistical 
im*tli(xl.s  used.  Perhaps  most  Important  Is  the  In¬ 
fluence  of  subjective  sampling  on  the  values  of  the 
correlation  coefficients  obtained  In  the  analysts. 
This  Influence  may  lx*  great,  txit  little  can  lx*  said 
of  it  before  much  more  detailed  work  is  done  on 
limited  aspects  of  this  general  study.  Rather  high 
values  for  the  critical  level  (e.g.  .05,  or  .10)  were 
used,  so  that  by  chance  alone,  one  would  expect 
some  of  the  differences  or  the  correlations  to  be 
sufficiently  great  to  be  judged  significant  erron¬ 
eously. 


RESULTS  OF  ANALYSIS  OF  VARIANCE  AND 
COVARIANCE  TESTS 

Rock  type  and  drainage  density:  The  type  of 
rock  has  an  apparent  influence  on  drainage  density, 
according  to  the  high  value  obtained  fur  F  In  the 
following  analysis  of  variance  table: 


Source  of 
Variation 

S.  S. 

d.f. 

M.  S.  F 

Among  rock 
types 

Within  rock 

57,371 

5 

11,474  56.7 

types 

12,140 

60 

202.3 

Total 

69,511 

65 

1,069 

F.  001=  4.76 


Six  categories  of  lithology  are  distinguished:  clastic 
(sandstone  and  alluvium),  shale,  limestone,  schist, 
granite  (includes  gneiss),  and  acid  volcanic  (in¬ 
cludes  pumice).  Inspection  of  the  means  of  D  for 
each  rock  type  shows  that  some  are  decidedly  dif¬ 
ferent  from  others  (Appendix  VI  A).  The  estimated 
effect  of  each  lithology  on  D  is: 

Shale:  150.5  Schist:  37.9 

Clastic:  -  7.8  Granitic:  -  3.86 

Limestone:  -12.5  Acid  volcanic:  -12.28 

These  values  are  probably  meaningless,  and  it 
should  be  noted  that  they  do  not  add  to  zero,  as  the 


sample  size  was  different  in  each  rase.  Hie  Impor 
tant  point  Is  that  shale  and  schist  have  average 
drainage  density  well  alxne  the  mean,  whereas 
limestone  and  acid  volcanic  rocks  are  well  lx*low 
tin*  mean. 

Scheffe's  method  of  judging  contrasts  was  used 
to  (inn  which  of  the  means  produce  the  significant 
results  (Scheffe ,  1953).  Significant  differences  in 
drainage  density  exist  between  the  average  for  shale 
and  schist,  the  two  lithologies  with  the  highest 
values,  and  the  average  for  granitic,  clastic,  vol¬ 
canic.  and  limestone,  and  also  between  the  average 
for  shale  and  the  average  for  schist.  No  significant 
difference  in  drainage  density  was  found  among 
granitic,  volcanic,  clastic,  or  limestone  (see  Appen¬ 
dix  VI  A). 

It  is  apparent  that  rock  type  has  perceptible  ef¬ 
fect  on  drainage  density  under  the  general  mountain- 
type  climate  in  this  study.  However,  the  sampling 
of  shale  and  limestone  basins  is  poor,  and  all  ba¬ 
sins  from  both  these  lithologies  are  taken  from 
areas  with  extremes  of  other  conditions.  For  in¬ 
stance,  all  the  shale  basins  lie  in  a  very  dry  cli¬ 
mate.  The  limestone  basins  are  in  the  Wasatch 
Range,  and  were  subjected  to  late  Pleistocene  up¬ 
lift  and  rejuvenation.  Tests  of  differing  rock  types 
within  the  same  general  climatic  type,  and  in  areas 
with  identical  histories  by  Miller  (1953,  p.  24,  and 
Table  1)  show  the  importance  of  lithology  in  deter¬ 
mining  drainage  density  in  a  humid  climate.  Further 
such  studies  may  show  that  elimination  of  variation 
due  to  climate  and  differing  geologic  histories  will 
reduce  the  error  variance  enough  to  allow  detection 
of  a  significant  difference  among  granite,  clastic 
rocks,  and  acid  volcanic  rocks,  that  in  this  study 
show  no  differences. 

Rock  type  and  infiltration  capacity:  Individual 
infiltration  capacity  readings  (not  basin  averages) 
were  tested  against  rock  type.  The  rock  types  dis¬ 
tinguished  are  the  same  as  for  drainage  density. 

The  analysis  of  variance  table  is: 


Source  of 
Variation 

S.  S. 

d.f. 

M.  S. 

F 

Among  rock  types 

67.89 

5 

13.58 

7.43 

Within  rock  types 

249.4 

135 

1.85 

Total 

317.3 

140 

2.27 

F  001(5,135)  =  4.10 
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t  hr  niur  ul  1  in  n  tgnil  ti  4C*.,  XIVil  l  indicate  a  that 
tliffe  re fu'ea  ill  inflltratlon-rapai'lly  !!.rji-.c  c\lt>l 
.Hi.oiv,  (hr  urlimn  rock  types  (!utMIUll»n  ul  tf*e 
riHV-lyjK-  means  (Appendix  VIA)  »lw»  n  llial  non  e 
differ  more  widely  from  the  grand  mean  than 
others.  The  estimated  effe  I  ul  each  litholuvty  <>*> 
infiltration  eapac  lly  lit 

Shale  - 1 .  -t  3  &li  i  rit  -1. 00 

Clastic  -0.53  Granitic.  0.81 

Limestone  0.52  Acid  volcanic  -0.31 

Much  reliance  cannot  Is*  placed  on  these  values.  It 
is  significant,  perhaps,  that  shale  and  schist  have 
infiltration  capacities  Ixrlow  the  mean,  whereas 
they  have  drainage  densities  well  alxive  the  mean. 
I.imes'one  and  granitic  rocks  have  infiltration 
capacities  alxive  the  mean,  hut  drainage  density 
l-elow  the  mean.  Add  volcanic  rocks  probably  haw 
infiltration  capacity  not  significantly  different  fron. 
the  mean,  though  drainage  density  below  the  mean, 
and  thus  act  somewhat  anomalously. 

The  high  value  of  F  is  the  result  of  the  high  in¬ 
filtration  rates  on  granitic  and  limestone  litholo¬ 
gies.  The  contrast  (granitic  and  limestone)  vs. 
(clastic  and  volcanic)  is  igniflcant  at  .05.  Also, 
the  contrast  (clastic  and  volcanic)  vs.  (shale  and 
schist),  formed  from  the  groups  with  the  lowest  in¬ 
filtration  rates,  was  judged  not  significant  at  .05 
(see  Appendix  VI  A). 

The  question  arises  as  to  whether  lithology 
affects  values  of  drainage  density  other  than 
through  the  value  of  the  infiltration  capacity.  A 
possible  method  of  approach  may  be  through  an  ex¬ 
periment  using  analysis  of  covariance,  correcting 
values  of  drainage  density  for  differences  in  infil¬ 
tration  capacity,  and  selecting  basins  with  the 
same  general  climate  and  history,  but  from  dif¬ 
ferent  lithologies.  The  data  obtained  in  this  study 
are  inadequate  to  attempt  this  analysis. 

Rock  type  and  valley-side  slopes:  The  individ¬ 
ual  slope  measurements  within  each  area  studied 
were  arranged  into  groups  corresponding  to  rock 
type,  without  regard  to  basin.  The  estimated  effect 
of  each  lithology  on  valley-slide  slope  is: 


Shale: 

7.9 

Schist: 

-4.4 

Clastic: 

1.5 

Granitic: 

-1.3 

Limestone: 

9.1 

Acid  Volcanic: 

-3.9 

Limes!"'  e  is  well  above  the  general  mean,  prob¬ 
ably  because  the  only  samples  are  fron:  basins  in 
the  Wasatch  Mountains  which  have  undergone  late 
Pleistocene  rejuvenation.  Shale  is  also  high,  per¬ 
haps  because  of  the  aridity  of  the  region  sampled. 
Schist  and  acid  volcanic  rocks  are  low.  To  test 
whether  these  differences  are  significant,  an  anal- 
ysis-of-variance  test  was  made: 


S<  ini  r  C  r  of 

Variation 

5.  8. 

d  1  M.  8  1 

Among  n 

type* 

Within  ri>c>. 

1 1.983. 48 

5  2396.69  18.99 

typed 

1 .o.29B.35 

87  4  126  19 

Total 

122.281.83 

879  . 

F  (5.874) 

4.10 

The  value  ul  F  obtained  is  significant  al  .001.  and 
the  evidence  points  to  the  conclusion  that  rock  type 
affects  valley-side  slopes. 

Be  au.se  the  average  slopes  on  shale  and  lime¬ 
stone  are  10  greater  than  on  the  other  lithologies, 
it  is  reasonably  suspected  that  the  high  val  es  are 
the  result  of  some  other  factors,  specifically  arid¬ 
ity  and  rejuvenation,  that  become  manifest  through 
the  poor  sampling  of  shale  and  limestone  basins, 
mentioned  above.  Therefore,  the  contrast  (clastic 
and  granitic)  vs.  (schist  and  volcanic)  was  analyzed, 
and  was  found  to  be  significant  at  .05.  Il  can  lie 
stated  that  lithology  probably  does  influence  valley- 
side  slopes  (see  Appendix  VI  A). 

Infiltration  capacity,  drainage  density,  and 
gullies:  Basins  for  which  infiltration-capacity  read¬ 
ings  are  available  were  grouped  according  to  whe¬ 
ther  they  had  gullied  slopes,  trenched  channels,  or 
unmodified  drainage  lines.  Infiltration  capacity  was 
used  as  a  supplementary  measmement  to  adjust  the 
drainage  densities  of  the  basins.  The  purpose  of  the 
lest  was  to  determine:  (1)  whether  basins  with  gul¬ 
lies  show  consistently  higher  drainage  densities 
than  basins  with  trenched  channels  or  unmodified 
drainage  lines,  when  variation  due  to  infiltration 
capacity  is  removed;  (2)  whether  basins  with  gullied 
slopes  show  consistently  lower  infiltration  capacity. 
The  analysis  of  covariance  table  is: 


Source  of 


Variation 

d.f. 

Sx2  £xy 

n 

sy- 

d.f. 

_sy'2 

M.S. 

Among 

means 

2 

1.4  -.74 

.42 

2 

.50 

.25 

Within 

groups 

20 

10.5  -12.03 

24.44 

19 

10.68 

.56 

Total 

22 

11.9  -12.77 

24.86 

21 

11.18 

For  X,  F 

=  1. 

33;  for  Y,  F  = 

0.17;  for  Y\  F  = 

0.45, 

where  X  is  infiltration  capacity,  Y  unadjusted 
drainage  density,  and  Y'  the  drainage  density  ad¬ 
justed  for  differences  in  infiltration  capacity.  The 
critical  value  of  F  ^(2,20)  is  3.49,  and  for  F 05 (2,19) 
is  3.52.  None  of  the  values  of  F  obtained  is  signifi¬ 
cant.  Though  there  is  an  increase  in  F  from  Y  to 


'i  ,  |(  in  not  enough  to  rwrtil  the  critical  \aluc  An 
can  Ixc  nct'n  from  the  rllrcln  (see  also  Appendix 

\  l  !M 


Edicts  of  gullying  Effects  of  gullying 

on  l>»g,  D 

on  infiltration 
capacity 

Channels 

Channels 

trenched  -0.237 

i  reached 

0.15 

Slopes  gullied  0.167 

Slopes  gullied. 

-0.28 

No  gullies  -0.035 

No  gullies: 

0.27 

The  mean  drainage  density 

for  gullied  las  Ins 

is 

higher  than  for  trenched  unmodified  lias  ins,  and  tiie 
mean  infiltration  capacity  Is  lower  than  for  the 
other  classes;  basins  with  trenched  channels  have 
lower  drainage  density  than  either  gullied  basins 
or  untrenched  and  unsullied  basins.  However  these 
differences  are  not  (treat  enough  to  lx*  judged  sig¬ 
nificant  at  .05  considering  the  amount  of  variability 
in  the  samples.  These  results  fit  theoretical  con¬ 
siderations,  and  that  may  tie  reason  to  Is'lieve  that 
with  larger  samples,  and  restricted  sanij  ling,  the 
results  would  prove  significant. 

There  is  present  a  significant  decrease  in  within- 
groups  mean  square  (variance)  when  the  effect  due 
to  differences  in  infiltration  is  removed.  The  unad¬ 
justed  mean  square  is  24.44,20  ^  1.22,  and  the 
adjusted  mean  square  is  0.56;  F  =  1.22/0.56  = 

2.18,  which  is  greater  than  the  critical  value  for 
.05.  This  indicates  that  a  significant  effect  of  in¬ 
filtration  capacity  on  drainage  density  is  present. 


RESULTS  OF  SIMPLE  CORRELATION  AND 
REGRESSION  ANALYSES 


Relations  among  morphometric  properties  of 
relief:  Preliminary  analysis  with  punched  cards 
indicated  that  significant  correlations  exist  among 
many  of  the  morphometric  properties  of  the  basins 
studied.  The  correlation  tables  thus  obtained  are 
given  in  Appendix  V  A.  Among  properties  involving 
relief,  it  was  found  that  maximum  valley-side  slopes 
are  correlated  positively  with  relative  relief  and 
ruggedness  number;  slopes  are  negatively  corre¬ 
lated  with  channel  frequency  and  drainage  density. 

/i  scatter  plot  of  valley-side  slope  (ord.nate) 
against  ruggedness  number  (abcissa)  showed  that  a 
large  number  of  points  representing  thirci-order 
basins  lay  above  the  average  line,  and  points  repre¬ 
senting  fifth-order  basins  lay  consistently  below 
the  average,  indicating  that  the  ruggedi.ess  number 
was  affected  by  the  greater  absolute  i « lief  of  fifth- 
order  basins.  For  purposes  of  estimating  valley- 
side  slopes,  a  better  fit  is  obtained  by  decreasing 


the  value  of  M  by  an  amount  depending  on  the  tan  In 
urdcj,  This  procedure  ran  I*-  explained  t>\  consult  r- 
tng  a  mode  1  an  a  ha\  ng  uniform  drainage  density 
throughout,  and  with  uniform  \  alley-side  slopes, 
for  all  orders  of  basins.  Then  the  absolute  relief 
within  a  basin  will  normally  lx1  greater  for  basins 
of  higher  orders,  simply  Ixcause  they  are  larger. 
Whereas  the  drainage  density  remains  constant 
(or  all  sizes  and  orders,  tin  ruggedness  numl>er 
will  increase  with  order.  Therefore,  to  obtain  the 
Ix'.st  relation  with  \  alley-side  slope,  which  is  con¬ 
stant.  the  value  of  H  obtained  for  each  basin  should 
!>e  divided  by  some  number  that  increases  with 
order,  so  that  a  constant  for  the  area  is  obtained. 

In  this  case  the  older  itself  was  used,  though  some 
function  of  the  order  number  might  give  Ixdler  re¬ 
sults,  perhaps  u"  or  u'  .  The  plots  of  0  against 

log  ^  and  log  R  arc  given  in  Appendix  V  A.  Tiie 

relation  of  valley-side  slopes  with  these  variables 
is  considered  below  in  connection  witli  climate  and 
soil  properties. 

The  correlations  of  O  with  log  F  and  log  D  indi¬ 
cate  that  an  inverse  relationship  exists  between 
slope  steepness  and  topographic  texture  (Appendix 
V  A)  The  absolute  value  of  the  correlation  co¬ 
efficient  of  O  with  log  D  (-.43)  is  greater  than  for 
O  with  log  F  (-.34),  indicating  that  perhaps  part 
of  the  decrease  in  texture  with  increasing  va^ey- 
side  slopes  is  the  result  of  foreshortening  of 
steeper  channels  by  projection  onto  the  map.  Fore¬ 
shortening  should  not  affect  the  number  of  channel 
segments,  and  as  channel  frequency  is  negatively 
correlated  with  O,  the  decrease  in  fineness  is 
probably  real.  A  possible  mathematical  model 
fitting  this  would  have  valley-side  slopes  inversely 
proportional  to  the  logarithm  of  both  drainage  den¬ 
sity  and  channel  frequency.  Conversely,  topographic 
texture  would  vary  as  a  negative  exponent  of  aver¬ 
age  valley-side  slope. 

In  addition,  D  and  F  are  negatively  correlated 
with  R,  as  would  be  expected,  since  R  is  highly 
correlated  with  O  (Appendix  V  A).  Because  corre¬ 
lation  is  low,  it  is  not  possible  to  determine  which 
type  of  function  gives  the  best  linear  scatter  plot. 

Which  of  the  morphometric  properties,  among 
D,  F,  R,  and  O  can  be  considered  as  controlling 
factors,  if  any  are,  is  not  definitely  known.  A 
plausible  hypothesis  is  that  R  and  either  F  or  D, 
or  both,  are  relatively  independent.  D  and  F  de¬ 
pend  on  rock  type  and  climate,  and  R  probably  de¬ 
pends  greatly  on  past  diastrophic  history,  Then 
H  u 

O  and  H  or  —  are  dependent.  Log  —  is  correlated 

positively  with  log  R,  (r  =  .666,  b  =  0.237),  as 
expected,  suggesting  the  mathematical  model  that 
ruggedness  number  corrected  for  basin  order  var¬ 
ies  about  as  the  one-fourth  power  of  relative  re¬ 
lief. 
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KFt.AtTO.N.s  AMONG  1't.AMMn  me  1*1*01* t:  II- 

TIKs  Of  tlAStNS 

Simple  rorr«'liU  *n  u|  two  \  r  t>l«- a  l»y  £ r c*j ped 
data  jihowrd  Hut  high  correlations  c\int  among 
lofjarUhn.fi  o(  13,  t)  i,  F,  and  Fj.  Three  power  func¬ 
tions  noli1  obtained  an  the  re-gresst  m  1 1 ne n  rrhliitf.; 
these  prope  rt  |e,s 

(ll  lav,  D  \  ti.  la*in  F  i  -  ,9-4,  l>  -  U.H3 , 

1.92.*  b  O.nOI)1 

(2)  Lou,  13,  vs.  Log,  13,  r  .68.  b  0.921. 

[3,  0.758  13  "ri. 

(3)  la>g  F,  vs  la>g  F  r  .97,  li  .975, 

F,  .  99  F  ■ . 

A  tost  of  significance  of  b  in  equation  (1)  alxive 
shows  that  It  is  not  significantly  different  from 
.857  (.40 's  P  .30).  and  2.303  b  Is  not  significantly 
different  fron.  2.00.  when  logarithms  are  removed. 
The  relation  F  kD  can  therefore  lx*  used  as  a 
mathematical  model  With  Increasingly  fine  tex- 
tured  topography,  F  Increases  proportionally  to  D‘ , 
showing  that  the  increase  in  total  length  of  drain¬ 
age  lines  is  due  to  an  increase  in  number  of  chan¬ 
nel  segments,  and  not  to  longer  segments  of  each 
order.  An  increase  in  drainage  density  is  then,  ac¬ 
cording  to  this  hypothesis,  a  change  of  scale  of 
topography,  resulting  in  smaller  basins  and  shorter 
channel  segments  of  each  order.  If  a  model  is  con¬ 
sidered  in  which  the  area  of  a  certain  basin  is 
constant,  that  is,  erosion  is  not  allowed  to  enlarge 
the  initial  drainage  area  by  extending  the  major 
divide,  then  the  equation  F  -  kD2can  be  trans¬ 
formed  into  I(1  =  - y  -,  with  the  restriction  that 

increases  of  3l„  are  due  solely  to  increases  in  D. 
Then  the  average  channel  lengths  of  all  orders  are 
inversely  proportional  to  the  total  length  of  chan¬ 
nels  at  a  given  instant,  during  the  time  that  the  in¬ 
crease  in  D  is  taking  place. t 

The  exponents  of  D  and  F.  in  equation  (2)  and 
(3)  respectively,  are  not  significantly  different 
from  1.00.  This  suggests  that  the  density  of  first- 
order  channels  increases  linearly  with  the  total 
drainage  density,  and  the  frequency  of  first-order 
channels  increases  linearly  with  the  total  channel 
frequency;  D1  =  kjD,  and  F1  =  k2F.  These  equa¬ 
tions  can  be  transformed  into  ly  =  k  y  3:  1  „ ,  and 
ni  =  koSnu,  where  increases  in  3  1„  and  3  n,,  are 
restricted  to  increases  in  fineness  of  topographic 
texture.  Combining  these  two  model  equations  gives 

2.308  -  log, 10. 

tThe  same  may  be  said  of  the  relations  among  several 
basins  of  the  same  planimetrie  size,  but  differing  in 
drainage  density. 


i.  -  *  AtdifdlfJo  to  Ullo  hq>o<f  8  the  JJrf 

Age  length  of  lltal'OUlrr  rhuifirls  lo  A  roftaljftl 
proportion  of  the  Average  length  of  of  All 

orders,  Uw  model  implies  then,  combining  with  l  hr 
fifnt  equation  obtained.  UvAt  JJrf  Ajjc  length  of  first 
order  channels  lo  imrrarh  prop-  rlional  lo  the 
total  length  o|  channel,  when  the  latter  IfU  TejncH 


during  an  Increase  of  drainage  density,  I 


k  « 

1 '  rt'.- 

if  area  in  constant.  A  corollary  of  thin  hypothesis 
it*  that  the  average  maximum  length  which  first- 
order  channels  attain  tx*fore  bifurcating  ts  less  tn 
areas  of  high  drainage  density  than  in  areas  of  low 
density  (Schumm.  1954.  p.  10). 

The  values  which  the  constants  of  prnfiort ton¬ 
ality  assume  In  the  imxlel  equations  stated  alxne 
are  abstract,  fundamental  properties  of  basins, 
and  probably  are  Influenced  by  many  of  the  govern¬ 
ing  agents  in  l»a.sln  rcgluien.  These  quantities 
should  lx*  of  great  Interest  for  further  investiga¬ 
tion. 


RELATION  OF  VALLEY-SIDE  SLOPES  TO 
OTHER  DASIN  CHARACTERISTICS 

The  effects  of  seven  variables  on  slope  values 
are  dealt  with  as  simple  two-variable  correlation 
and  regression  relations,  and  also  combiner’  in 
multiple  correlation  coefficients.  The  variables 
are  (logarithms  taken  where  necessary  for  linear 
plot,  Y  vs.  XJ; 

X  |  wet  soil  strength/1000 
X  2  infiltration  capacity 
X  t  Log  Ruggedness  number/order 
X  .,  Log  (P-E) 

X  s  Log  relative  relief 
X  ()  Roughness  number 
X7  5-year,  1-hour  rain  5 
Y  Average  maximum  valley-side  slope. 

The  simple,  two-variabie  correlation  coeffi¬ 
cients  obtained  from  analysis  of  both  grouped 
and  ungrouped  data  are  given  in  Appendix  V  and 
below: 

rY1  =  -.423;  G  =  -.  756  +  28.26;  sy  x  =  7.4 

r y2  =  .746;  6  =1.312  f  +  20.7;  sy  x  =  4.9 

ryi  =  .777;  6  =23.71  log  “  +  30.84;  sy  x  =4.6 

rY4  =  .289;  6  =  9.90  log  (P-E)  +  10.8;  sy  .  x  =  7.0 

rY5  =  .625;  G  =18.69  log  R  +  13.1;  sY  X  =5.7 


♦True  for  third  and  fourth-order  basins;  experimentally 
verified  since  writing. 
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Multiple  correlation 
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of  three  tarl- 
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=  .916; 

II 

CD 
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Circe 
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le  partial 
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coefficients 

were  1 

found  by 

combining  the  simple 

correlation 

coeffli 

dents: 

rn  • 

=  .620 

r  Yl.? 

=  .797 

r  ^  j  i  -.532 

-  .230 

r 

=  .804 

r  ti .  ?  ="f’14 

r  1  3  .  1 

=  .769 

r  v?.: 

=  .893 

rn  ,  =  .836 

These  coefficients  are  not  part  of  a  complete 
multiple  correlation  study.  They  are  used  only  as 
estimates  of  the  population  correlation  coefficients. 
Thus,  some  of  the  coefficients  calculated  from  lar¬ 
ger  samples  are  judged  significant  though  smaller 
than  rVo  which  is  not  significant. 

The  only  variable  not  judged  significant  is  the 
roughness  number.  However,  the  negative  value  of 
the  regression  coefficient  obtained  is  cons  stent 
with  a  theory  to  be  developed  below  to  explain  the 
inverse  relation  between  slope  angle  and  wet  soil 
strength;  both  roughness  and  strength  depend  on 
the  soil  texture  and  adequate  sampling  should  de¬ 
tect  a  positive  relation  between  them. 

The  highest  correlations  were  found  between  G 

and  log  — ,  also  O  and  f.  These  have  lowest  stand- 

o  u 

ard  error  of  estimate,  that  being  on  the  order  of 
the  usual  standard  deviation  of  measurements  with- 

pj 

in  individual  basins.  Log  —  alone  explains!  60.4% 

of  the  variation  in  O  ;  f  alone  explains  55.6%.  The 
third  highest  correlation  was  found  to  be  with  log 
R,  agreeing  with  previous  findings  (Schumm,  1954, 


♦The  computation  of  multiple  correlation  coefficients 
was  restricted  here  because  not  all  intercorrelations 
among  the  X  j  were  computed.  The  multiple  coefficients 
chosen  for  presentation  show  the  importance  of  the  in¬ 
fluence  of  infiltration  capacity,  log  ruggedness  number/ 
order,  log(P-E)  index,  and  runoff  intensity  on  valley- 
side  slopes.  These  four  variables  are  furthermore  of 
greatest  interest  in  this  study,  summarising  hydrologic, 
diastrophic,  and  climatic  factors. 

fBy  “explains,”  it  is  meant  that  knowledge  of  X  enables 
us  to  predict  Y  with  error  equal  to  (1  -  r2)  of  previous 
error.  Above,  60.4  100r2,  etc. 


|>.  23.  ami  fig.  I  Si.  laig  1<  alum  explains  39.1  {>  of 
the  variation  m  o.  Mien  expressed  in  degrees, 
'..illes-mde  slope  \nries  proportionally  to  the  loga¬ 
rithm  nl  tlie  relative  relief,  the  logarithm  of  rug 
gednens  nun. U  r  corrected  for  order,  and  infiltra¬ 
tion  capacity.  The  lowest  correlation  was  found 
with  log  li’-K),  liut  II  was  judged  significant  Iw- 
cause  It  was  leised  on  a  large  sample. 

Three  multiple  correlation  coefficients  were  cal¬ 
culated  from  the  simple  coefficients  (Croxton,  195J, 
p  1’.  2).  The  values  obtained  were  rather  high,  show¬ 
ing  that  a  large  portion  of  the  variation  in  the  valley 
side  slopes  is  explained  by  consideration  of  only 
two  independent  variables.  For  instance.  Infiltra¬ 
tion  capacity,  and  the  five-year,  one-hour  rainfall 
amount  minus  the  infiltration  capacity  plus  five  (q) 
(which  take  into  account  infiltration  and  the  relation 
of  Infiltration  to  average  intensity  of  rain  storms), 
together  explain  84.3%  of  the  variation  in  the  valley- 

side  slope.”  Infiltration  capacity  and  —  together 

explain  83.8%  of  the  variation;  infiltration  capacity 
and  log  (P-E)  together  explain  72.4%  of  the  varia¬ 
tion  in  O. 

Nine  of  the  possible  partial  correlation  coeffi¬ 
cients  were  computed,  and  excepting  ryr,  y  ,  all 
are  significantly  different  from  zero.  This  means 
that  log  R  does  not  improve  the  estimation  of  O 

LI 

when  log  is  already  used.  However,  consideration 

of  log  —  when  log  R  is  already  used  does  signifi¬ 
cantly  improve  the  estimation  of  O.  Consideration 
of  q,  when  f  is  already  taken  into  account,  explains 
about  64.4%  of  the  variation  in  O  unexplained  by  f 
alone.  In  this  connection,  a  correlation  exists  be¬ 
tween  f  and  log  (P-E),  (r21  =  .750),  a  relation  well 
known,  and  probably  best  explained  by  the  increased 
amounts  of  humus  and  bacteria  in  soils  withmoister 
environments. 

The  results  of  greatest  interest  here  are  the  re¬ 
lations  between  O  and  f,  and  O  and  q.  It  is  apparent 
that  greater  average  amounts  and  intensities  of  run¬ 
off  are  associated  with  lower  slopes.  Conversely, 
areas  with  high  infiltration  capacity  have  low  run¬ 
off,  and  ceteris  paribus,  low  rates  of  slope  erosion; 
valley-side  slopes  in  these  areas  tend  to  become 
and  remain  steep.  In  discussing  a  continuous  down- 
slope  profile,  Gilbert  (1877,  p.  117)  stated  that  the 
declivity  of  slope  is  inversely  proportional  to  the 
amount  of  water  flowing  over  it,  an  observation  that 

♦The  multiple  correlation  coefficient  and  predicted  slope 
values  obtained  from  the  use  of  infiltration  capacity  (0 
and  5-year,  1-hour  rainfall  amount  (P)  as  independent 
variables  would  be  identical  with  those  obtained  from  any 
two  independent  linear  combinations  of  f  and  P,  such  as 
f  and  q  -  (P  -  f  •  5),  used  here.  The  reason  for  using  q 
is  that  it  has  physical  meaning,  viz.,  the  intensity  of  run¬ 
off  expected  every  five  years  from  a  one-hour  rain. 
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appear n  In  hold  In  gnural  Kurlherinofr  Ku  ne-  data 
f*how  that  the  relation  Ix-Uirn  nlupr  and  f  nr  q  i» 
arithmetically  linear 

0  *  »i  -  kl  '!•  0  a:  .  k-  f 

1  he  following  working  hypotheiilji  i.s  proponed 
where  infiltration  rapacity  la  low,  greater  t*  ran  ion 
accompanies  greater  runoff,  and  channels  are 
eroded  slowly  in  relation  to  the  adjacent  slopes, 
because  much  debris  from  the  slopes  is  brought 
down  into  them.  As  a  result,  the  slopes  are  lowered 
more  rapidly  than  the  channels,  and  the  ancle  of 
slope  declines.  With  higher  Infiltration  capacity, 
little  runoff  occurs,  and  the  valley-side  slopes  are 
not  eroded  so  rapidly  as  in  the  previous  case.  Con¬ 
centration  of  the  runoff  in  channels  greatly  in¬ 
creases  its  eroding  efficiency.  An  Increment  of 
water  to  channel  flow  comes  by  interflow  (subsur¬ 
face),  and  this  addition  serves  to  increase  erosion 
in  the  channels,  but  the  water  thus  added  did  not 
have  the  opportunity  to  erode  the  valley  sides.  Thus 
when  the  channel  below  a  segment  of  slope  is 
deepened,  a  proportionate  amount  of  erosion  has 
not  occurred  on  the  slope  above.  The  tendency  is 
thus  to  remove  the  support  of  the  surface  of  the 
valley  side  afforded  by  the  toe,  thereby  allowing 
an  acceleration  of  creep,  and  consequent  steepen¬ 
ing.  The  quantity  of  interflow  is  affected  by  the 
structure  of  the  soil,  and  is  favored  by  shallow 
soils,  with  presence  of  parent  rock  or  other  im- 


Figure  8.  Hypothetical  channels  (A)  with  low  infiltration 
capacity,  high  runoff  and  erosion  rates,  low  interflow, 
debris-laden  channel  and  low  valley-side  slopes;  (B)  with 
high  infiltration  capacity,  low  runoff  and  erosion  rates, 
high  interflow,  eroded  channel,  and  high  valley -side 
slopes. 


p«- r it. cable  layer  a  nfw»rl  distance  l*-|ow  ihc  sur¬ 
face  u{  the  ground  Thin  would  toe  the  lane  in  *nui  t 
if  the  nett. lartd  and  arid  went  (f  inrilcy,  Kohler 
and  Paul  turn,  1949,  p.  3881, 

Climatic  rondltliwm  enter  fur'hcr.  Ix-caune  art 
the  average  rainfall  intensity  Increase.-;,  the  aver¬ 
age  amount  and  ir.tciiNlty  of  runoff  mere  m-.-t.  Fur¬ 
thermore.  greater  P- K  index  tends  to  increase  in¬ 
filtration  capacity,  thus  pnxtucmg  steeper  slopes. 

The  effect  of  soil  strength  on  the  average  valley 
side  slope  maintained  m  an  area  of  mature  topog¬ 
raphy  could  lx-  deduced  from  sc-.iial  reasonable 
premises,  but  the  results  would  conflict  unless  the 
quantitative  Importance  of  all  pertinent  variables 
is  known.  A  theorv  m  11  Ixr  constructed  that  Is  con¬ 
sistent  with  the  results  obtained  by  statistical 
analysis  of  the  measurements  made  In  this  study, 

Init  confirmation  of  the  theory  depends  on  further 
hydro'oglc  research.  We  know  (1)  that  O  varies 
inversely  with  S.  (wet  soil  strength),  because 
r^  |  =  -.<123  is  significantly  less  than  zero;  (2)  that 
the  correlation  between  O  and  S.  is  improved 
when  the  effect  of  inf  i  It  rat  ion  capacity  is  removed 
(r y |  o  =  -.532);  (3)  that  wet  soil  strength  is  greater 
for  coarser  soils  (confirmed  by  a  one-way  analysis- 
of-variance  test,  not  reproduced  here).  Excluding 
cliffs  from  consideration,  the  effect  of  soil  strength 
must  be  through  erosive  processes  rather  than 
stability,  because  weaker  soils  tend  to  lx-  stable  at 
lower  slope  angles,  in  opposition  to  (1)  above. 

Analysis  of  covariance  of  the  data  in  Table  4 
gives  the  following  tabic; 

d.E  vx2  vXy _ vy  2  d.f.  vy^2  M.  J7 

Among 

means  2  3.25  40.82  714.0  2  521.2  260.6 

Within 

groups  20  12.59  18.02  686.4  19  660.6  34.7 


Total  22  15.84  58.84  1400.4  21  1181.8 


where  X  =  f,  Y  =  0 .  The  above  was  obtained  by 
grouping  0  and  f  into  three  classes  depending  on 
soil  texture.  Fine-grained  soils  give  high  values  of 
0,  and  coarse-grained  soils  give  low  values  of  0. 
The  F-ratio  for  slope  means  unadjusted  for  infil¬ 
tration  capacity  is  10.4  (Significant  at  .001),  and 
the  F-ratio  for  slope  means  after  adjustment  for 
infiltration  capacity  is  7.52  (significant  at  .01). 
Therefore,  soil  texture  influences  valley-side 
slope  through  some  agency  other  than  infiltration 
capacity  alone.  The  correlation  between  0  and  Sw 
is  probably  due  to  this  agency. 

In  the  following  proposal  of  a  theory  to  explain 
the  inverse  relation  befween  O  and  Sw,  it  is  as¬ 
sumed  (1)  that  the  hypothetical  areas  under  consid¬ 
eration  have  identical  climates,  (2)  that  the  bed 
rock  near  the  surface  is  essentially  homogeneous 
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afHl  dx-n  II  I!  differ  gfratly  In  Hn-ifv4lltC.il  i  on.po- 
nliliin  f  r  t»t».  live  lop  ft  o’  l  layer,  ami  (3)  Oval  the  art-a 
tl  wn  not  lie  « ire  an  Irom  a  mmrre  of  ilelvrln 

of  conlraftUrii*  calilwr  l<>  Ifut  pnxlurrd  within  the 
area. 

In  an  area  of  low-  S„.  wilt  and  clay  predominate 
in  the  Hull,  and  j: ravel  irt  alwenl.  IX-brltt  brought  in¬ 
to  the  channelfi  by  runoff  from  the  Nlope.s  remain.*! 
in  sum  pens  ion  and  is  immediately  transported  out 
of  the  area.  Consequently,  the  channels  remain  free 
from  bed  load  and  are  subject  to  immediate  erosion 
b>  the  concentrated  .stream  of  water.  Mule  runoff 
lasts,  the  channels  are  continually  deepened  and  or 
widened,  and  are  not  refilled  as  discharge  wanes— 
silt  and  especially  clay  remaining  In  suspension 
even  In  low  velocities.  Because  of  the  greater  ef¬ 
ficiency  for  erosion  of  the  concentrated  stream 
over  slope  wash,  the  channels  are  lowered  faster 
than  the  slopes,  the  toes  of  the  slopes  are  removed, 
and  the  adjacent  slopes  are  steepened  until  a  steady 
state  is  achieved  through  a  reduction  of  the  rate  of 
erosion  in  the  channels  relative  to  the  rate  of 
erosion  of  the  slopes.  The  self-limiting  mechan¬ 
ism  is  the  result  of  increased  erosion  and  creep 
on  the  sleeper  *d  slopes,  causing  channel  gradient 
to  lx?  lowered  to  that  value  necessary  to  transport 
the  small  quantity  of  fine-caliber  bed  load.  In  such 
areas,  at  steady  state  the  slope  orthogonals  will 
intersect  the  channel  at  high  angles,  and  the  inflec¬ 
tion  angles  will  be  acute. 

In  an  area  of  high  S»  gravel  predominates  in  the 
soil.  Slope  runoff  sweeps  debris  into  the  subjacent 
channel  where  coarse  material  accumulates  as  bed 
load.  Assuming  a  case  in  which  valley-side  slopes  are 
initially  steeper  than  equilibrium  requires  (e.g. 
where  channels  are  trenched  or  rejuvenated),  the 
greater  declivity  of  the  slopes  permits  movement 
of  gravel  and  cobbles  by  rolling  and  sliding  under 
relatively  low  down-slope  discharges,  but  in  the 
channel,  higher  velocities  and  discharges  are  re¬ 
quired  for  transportation  because  the  down-stream 
gravitational  component  is  low.  Down-cutting  of  the 
bed-rock  channel  cannot  occur  until  the  bed  load 
is  removed,  and  except  where  the  channel  is  sup¬ 
plied  by  ground  water,  the  bed  load  is  replenished 
by  material  brought  down  from  the  slopes  when¬ 
ever  sufficient  discharge  occurs  to  transport  all 
of  the  bed  load  already  present.  Thus,  channel  down¬ 
cutting  is  a  relatively  rare  event,  but  th*-*  slopes  can 
be  eroded  frequently.  As  channel  discharge  wanes, 
after  a  runoff-producing  rain,  any  depressions  that 
had  formed  in  the  channel  are  refilled  by  removal 
of  debris  from  suspension.  Erosion  of  the  slope 
toes  by  lateral  corrasion  can  occur  at  the  outside 
of  sharp  channel  bends,  with  consequent  local  steep¬ 
ening  of  the  cut  bank  slope,  but  a  compensating 
lowering  of  opposite  slopes  results  in  no  quantita¬ 
tively  important  change  over  the  basin  as  a  whole 


fStrahlrr  I960,  pp.  812,  813)  Tin-  n  loped  ^n-  there 
lore  lowered  at  a  greater  rate  than  the  (hanneln, 
and  reach  a  nlcady  state  at  mime  relatively  low 
angle,  when  the  down- slope  component  of  gravity  Is 
so  adjusted  aw  to  allow  only  that  amount  of  roarne  de  - 
tirls  to  enter  the  channel  which  can  lx*  transported 
by  stiean.  flow  The  channel  gradient  necessary 
to  transport  the  large  amount  of  coarse  lx-d  load  is 
relatively  stee|>,  thus  the  angle  txdween  slope  oi  - 
thogonals  and  channel  will  lx*  very  acute,  and  the 
inflection  angles  obtuse. 

In  extreme  cases,  if  the  amount  of  debris  sup¬ 
plied  to  the  i  hannels  is  quite  great.  It  seems  prob¬ 
able  that  the  channel  will  aggrade  rapidly  and  as¬ 
sume  a  braiding  habit.  The  slope  toes  along  the 
entire  course  of  the  channel  may  then  be  cut  by 
lateral  corrasion  and  the  slopes  lx-  maintained  at  a 
steep  angle  until  the  maximum  breadth  possible  for 
the  drainage  area  has  lx?en  attained.  Basins  with 
this  environment  were  not  included  in  the  areas 
studied. 

The  mechanism  proposed  above  for  determining 
the  relative  rates  of  slope  and  channel  erosion,  i.e. 
the  abundance  of  bed  load,  is  the  sail,  as  that  pro¬ 
posed  to  explain  the  relation  of  slope  angles  to  in¬ 
filtration  rates.  However,  areas  most  favorable  for 
the  accumulation  of  bed  load  because  of  low  infiltra¬ 
tion  capacity-high  runoff  intensity  are  usually  least 
favorable  for  the  production  of  coarse  debris  be¬ 
cause  of  the  prevalence  of  silt  and  clay.  Therefore, 
the  effect  of  infiltration  capacity  is  generally  op¬ 
posite  the  effect  of  debris  availability,  though  both 
depend  on  soil  texture  and  act  through  the  same 
mechanism.  This  explains  the  indicated  increased 
correlation  between  O  and  Sw  when  the  effect  of  in¬ 
filtration  capacity  is  removed.  Infiltration  capacity 
depends  also  on  climate,  so  the  net  effect  of  soil 
texture  cannot  be  judged  qualitatively  except  per¬ 
haps  in  extreme  cases.  The  equilibrium  valley- 
side  slopes  for  given  conditions  of  climate,  infiltra¬ 
tion  capacity,  soil  texture,  etc.  are  best  estimated 
quantitatively  by  use  of  regression  formulas. 

Much  further  investigation  is  necessary  before 
definite  hypotheses  can  be  formulated  to  include 
more  than  a  few  of  the  variables  at  one  time  for  t’.e 
entire  environmental  conditions  of  soil,  climate, 
history,  structure,  and  slope. 


TEXTURE  OF  TOPOGRAPHY 
AS  AFFECTED  BY  CLIMATIC 
AND  MANTLE  CHARACTERISTICS 

Information  from  23  basins  subjected  to  full  field 
investigation  was  arranged  for  multiple  regression 
and  correlation  analysis.  Drainage  density  (log  D) 
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The  multiple  correlation  coefficient  is  Rv  12i.is  = 
0.965,  and  thus  93.2%  of  the  variation  in  L)  is  ex¬ 
plained  by  including  only  climatic  and  surface 
p*operties.  If  the  effects  of  certain  morphologic 
features,  such  as  R  or  O  were  considered  in  addi¬ 
tion,  possibly  a  greater  amount  of  variation  would 
be  explained.  An  analysis  of  variance  test  for  the 
significance  of  Rv  (Anderson  and  Bancroft, 

1952,  p.  163,  172)  is: 


Source  of 


Variation 

d.f. 

S.  S. 

M.  S. 

F 

Regression 

5 

4.489 

0.898 

46.29 

Error 

17 

.329 

0.0194 

Total 

22 

4.818 

0.219 

Ry.123.tf,  is  significant  at  .001. 


*  "  cse  particular  independent  variables  are  used  be- 
c;  ,se  they  summarize  the  climatic,  hydrologic,  and  sur¬ 
face  properties  that  determine  the  topographic  texture 
of  a  region.  From  two-variable  plots,  it  was  found  that 
log  D  is  a  linear  function  of  the  logarithms  of  f,  M,  and 
q,  so  the  log  transformation  is  used  for  these.  Log  D  is 
linear  with  b/  100,  anil  therefore  the  logarithmic  trans¬ 
formation  was  not  used  lor  it. 


Multiple  cor  relation  ffuu-nt.n  m-. 
log  (I’-Kj,  log  l.  and  h-g  q  .re 

It,  .9*6 

It,  i:  ,960 

St  levied  partial  correlation  coefficients  are  * 

a)  Three  variables 


ril  : 

-.842 

rn  : 

-  698 

r>i  i 

-.969 

rl'  2 

-  ,064(  N.S. ) 

r  >  i  •  * 

-  102 

rtv: 

.3  53  (N.S.) 

ri?  I 

-  .  l58(N.S.)f 

r3t  l 

191  (N.S.) 

r  i  t  1  “ 

.34  1  (N’.S.) 

r  >  V  2 

.704 

rl  t.l  - 

-.440 

r  v  r> .  l  = 

-.214(N.S.) 

rW2  = 

-  .234(N.S.) 

Four  va 

rlables: 

r  11.21  ' 

-.689 

r  1 5 . 2 1  = 

-.463 

r  VI  .25 

-.883 

r 12.15  = 

-.518 

r  VI.  12  = 

.306(N.S.) 

r  Y  t.  12  = 

- .  424  (.10) 

r  V  5 . 1 2  = 

-.533 

r  VS. 21  = 

.0223(N.S.) 

Five  variables: 

r  Y 1 . 12  5 

-  -.766 

r  Y1. 125  : 

=  ,0139(N.S.) 

The  following  inferences  may  be  made  from  the 
quantities  tabulated  and  presented  above: 

1.  Log  D  =  Y  shows  significant  correlations 
with  (1)  log  (P-E)  index,  (2)  log  percent  bare  area/ 
100,  (3)  log  infiltration  capacity,  and  (4)  log  runoff 
intensity,  in  decreasing  order  of  importance  as  de¬ 
termined  from  the  simple  correlation  coefficients. 
Lol,  roughness  number  alone  appears  to  have  no 
correlation  with  Y.  These  variables,  together,  ex¬ 
plain  all  but  6.8%  of  the  vt  i : ° t  i on  in  Y,  and  we  may 
state  that  the  climatic,  hydro  ;gic,  and  surface 
properties  of  an  area  are  of  overwhelming  impor- 

*The  partial  correlation  coefficients  presented  were 
chosen  mainly  to  clarify  the  relative  importance  of  the 
P-E  index,  infiltration  capacity,  and  runoff  intensity  in 
determining  drainage  density. 

t(N.S.)  denotes  not  significant  at  .10;  (.10)  denotes  sig¬ 
nificant  at  10  but  not  .05;  all  others  significant  at  .05  or 
less . 
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Umr  in  m-ti-r  mining  topographic  U-xlurr.  Further - 
itiurc  their  influcnccri  ran  (*•  detcrlcd  by  relatively 
simple  method*  il  field  mcan'iremeutn. 

2  The  '.aim  i(  II  ^  jn<  indicate.*  that  uf  the  93.21 
of  i In*  uruilun  1  ’>  ^  cxplati  ■  .1  by  r<in*lderat  ton  of 
all  fi\c  Independent  variables,  92  L  In  explained  by 
X  .  X  ,  and  X'.  percent  lure  area  and  surface 
roughness  are  very  minor  influences,  when  the 
others  aie  already  considered  (note  also,  r ^ . 

.01391 

3.  1  lie  regression  coefficients,  which  are  par¬ 
tial  derivatives  of  Y  with  respect  to  each  Independ¬ 
ent  variable,  with  the  others  held  constant,  (fives 
the  following  decreasing  order  of  \  irtance  to 
the  independent  variables:  (1)  log  runoff  Intensity, 

(2)  log  (P-E)  index,  (3)  log  Infiltration  capacity. 

(4)  log  percent  lure  area  100.  and  (5)  log  rough¬ 
ness  number.  The  discrepancy  in  order  A  impor¬ 
tance  results  m  part  from  Intercorrelations  among 
the  independent  variables,  and  in  part  from  the 
scale  of  the  measures. 

4.  The  partial  correlation  coefficients  indicate 
that  when  X(  is  already  considered,  the  only  vari¬ 
able  that  is  significant  is  X  t  (ruughness).  If  this 
is  not  a  spurious  result,  then  when  the  P-E  index 
is  held  constant,  drainage  density  is  inversely  re¬ 
lated  to  roughness.  (Or,  areas  of  greater  rough¬ 
ness  are  associated  with  lower  drainage  density.) 

It  should  be  noted  that  X  t  snows  the  lowest  simple 
correlation  with  Xj  oi  all  the  independent  variables. 
Conversely,  when  the  effects  of  X2,  X-,,  Xr.  X*  and 
X^,  Xn  and  Xs,  and  Xt,  X3,  and  Xr,  are  removed, 

Xj  still  has  a  significant  correlation  with  drainage 
density.  The  conclusion  d  awn  is  that  the  P-E  index 
influences  drainage  density  through  the  agencies 
of  infiltration  capacity,  amount  of  cover,  and  run¬ 
off  intensity,  but  probably  not  through  surface 
roughness.  In  addition,  the  P-E  index  influences 
drainage  density  by  other  mechanisms,  not  ex¬ 
plained  here. 

5.  The  infiltration  capacity  of  an  area  is  deter¬ 
mined  by  the  P-E  index  of  the  region,  the  surface 
cover,  and  by  the  type  of  rock  from  which  the  soil 
was  developed.  The  amount  of  surface  cover  is  de¬ 
termined  by  the  P-E  index  also.  The  runoff  inten¬ 
sity  depends  on  average  rainfall  intensity  (a  func¬ 
tion  of  the  same  factors  that  determine  P-E  index), 
infiltration  capacity,  and  hence  on  surface  cover 
and  the  P-E  index.  In  equilibrium  conditions,  the 
drainage  density  will  attain  and  stay  at  a  value  de¬ 
pendent  on  all  these  variables,  and  others  not  con¬ 
sidered  here.  However,  if  say,  the  infiltration 
capacity  and  percent  bare  area  were  changed  by 
fire,  grazing,  cultivation,  etc.,  to  values  not  in 
agreement  with  the  P-E  index,  a  change  in  drainage 
density  would  result,  and  the  new  value  would  no 
longer  agree  with  that  expected  from  consideration 
of  the  P-E  index  alone. 


l»  F  N F HA  1.1  / A T ION'S  UN  IN\  E>"I  IUAT1V  f 
PHOCEDfUEi) 


Thv  following  nr>'  generalt/.nllcms  which  seen 
appropriate  concerning  the  philosophy  of  the  invcs- 
t igal ion  as  a  whole 

1.  It  has  ln'rii  sho*'n  that  when  various  proper 
ties  of  drainage  lias  ins  are  analyzed,  occurrences 
are  found  that  would  not  reasonably  I*-  expected  to 
result  from  operation  of  chance,  or  sampling  error 
alone. 

2.  It  was  proposed  that  lliest  non-chance  occur¬ 
rences  represent  In  some  manner  true  laws  of  na¬ 
ture,  hence  mathematical  mcxlels  were  suggested 
consistent  with  what  is  known  about  these  laws 
from  sampling  and  statistical  analysis. 

3.  It  is  the  nature  of  statistical  laws,  such  as 
obtained,  that  a  single  observation  contrary  to  the 
supposed  law  will  not  refute  it.  as  the  element  of 
random  t  rror  is  always  pri  sent  in  any  observation 
or  measurement,  and  is  here  presumed  to  exist  in 
nature  when  a  large  system  Is  considered. 

4.  Because  it  is  impossible  to  consider  all  as¬ 
pects  of  the  environment  of  a  drainage  system,  or 
any  physical  system  whatsoever,  we  will  never 
know  that  all  factors  influencing  observations  are 
included  in  the  analysis.  It  seems  most  unlikely 
that  enough  factors  could  ever  be  considered  at 
once  to  predict,  completely  without  error,  the 
lengths,  gradients,  and  positions  of  all  channels  in 
a  basin;  adherence  to  the  notion  that  perfect  pre¬ 
diction  could  be  obtained  by  finite  extension  of 
purely  analytical  differential  equations  is  then 
rather  useless.  A  better  point  of  view  is  to  acknowl¬ 
edge  that  some  natural  agencies  operate  in  a  man¬ 
ner  that  cannot  be  predicted  except  by  the  use  of 
probability  theory,  and  then  build  a  theory  based  on 
statistical  mechanics  rather  than  on  Cartesian- 
Newtonian  mechanics. 

5.  Some  of  the  non- random  occurrences,  the 
correlations  obtained  above,  may  reflect  the  agency 
of  cause  and  effect.  However,  because  of  the  com¬ 
plexity  of  the  interconnections  of  causative  factors 
and  the  presence  of  random  error,  simple  cause 
and  effect— meaning  relations  of  forces,  motions, 
and  bodies— seems  inappropriate,  and  instead  a 
complex  feedback  among  cause-like  agencies  and 
effect-like  agencies  seems  more  pertinent  (in  this 
case,  properties  of  climate,  lithology,  history, 
vegetation,  mantle,  and  form),  when  the  totality  of 
environment  is  included. 

6.  This  does  not  mean  that  physical  laws  of 
cause  and  effect  do  not  operate,  or  that  methods  of 
analysis  by  differential  equations  cannot  be  applied. 
Simple  mechanical  models  of  a  natural  system  are 
satisfactory  only  in  those  very  small  parts  limited 
both  in  time  and  physical  boundaries. 


Summary  o!  Conclusions 


Hu-  (opographli  irxiurt-  ami  average  vallcv- 
side  slopes  of  fluvinlly-iont  rolled  landscapes  arc 
i  Hated  to  Hie  causative  and  determinative  factors 
<>f  r 'llmate  and  lithology.  However,  it  is  difficult, 
or  even  meaningless  to  consider  Isolated  elements 
of  the  environment  of  a  drainage  Ixisin  in  connec¬ 
tion  ’a  ith  'hese  geometric  properties  of  the  has  In 
Ik-c -au.se  the  geologic  and  climatic  elements  act  In 
complex  ways  through  the  agencies  of  vegetal 
growth,  soil  formation,  runoff  and  erosion.  Infil¬ 
tration,  and  soil  creep.  Thus  In  this  study,  as  many 
elements  as  practicable  were  analyzed  concur¬ 
rently. 

Drainage  density  is  a  negative  power  function 
of  the  P- L  index,  a  measure  of  the  availability  of 
moisture  to  vegetation,  but  also  depends  on  the 
percentage  of  bare  area,  infiltration  capacity,  soil 
strength,  and  Intensity  of  runoff.  These  in  turn  de¬ 
pend  on  the  P-E  index  and  are  the  agents  by  which 
climate  influences  drainage  density.  Lithology  in¬ 
fluences  drainage  density  by  determining  'he 
amount  of  clay  in  the  soil  (though  this  depends 
also  on  the  climate),  and  so  affects  the  infiltration 
capacity  and  soil  strength.  Perhaps,  also,  lithol¬ 
ogy  indirectly  affects  climate,  because  in  the  areas 


studied  moist  climates  are  confined  to  mountain¬ 
ous  regions,  and  thus  are  found  primarily  on  re¬ 
sistant  rinks.  Lastly,  recency  of  uplift  influences 
drainage  density  and  valley-side  slopes  lvcau.se 
areas  of  great  immediate  relief  have  large  relief 
ratios  and  generally  coarse  texture.  To  develop 
fine  texture  In  such  an  area  would  require  Impos¬ 
sibly  seep  slopes  (or  equilibrium  conditions. 

The  valley-side  slopes  of  drainage  basins  vary 
with  lithology,  l>elng  steepest  on  shale  and  perhaps 
limestone,  and  gentlest  on  schist  and  acid  volcanic 
lithologies.  Steep  slopes  are  associated  with  large 
relief  ratio,  high  infiltration  rate,  low  wet  soil 
strength,  low  runoff  intensity,  and  high  P-E  index. 
This  suggests  that  steep  slopes  are  produced  and 
maintained  in  areas  favoring  channel  deepening 
over  slope  erosion.  In  such  areas,  water  passing 
into  the  soil  by  infiltration  has  no  opportunity  to 
erode  the  slope  and  contribute  detritus  to  the  chan¬ 
nel  load,  but  emerges  directly  into  the  channel  and 
is  immediately  available  for  channel  erosion.  Thus, 
basins  with  steep  slopes  should  be  expected,  on  the 
average,  to  have  low  drainage  density.  This  rela¬ 
tion  was  found  in  the  correlation  analysis  of  the 
two  variables. 
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Appendix  I 


DKim  AT  ION  OF  RELATIONS  AMONG  GLUME- 

TIUCAL  Ft  FME'*  la  OF  VALLEYS 

Valley-side  slopes,  inflection  ancle,  and  channel 
gradient  are  related  by  laws  governing  lines  and 
planes,  assuming  that  the  \ alley  sides  are  planar, 
ard  therefore  meet  in  a  straight  line  that  Is  the 
channel.  Consider  the  contour  map  (a)  showing  a 
portion  of  a  basin  with  valley  sides  ,-j  and  r2,  with 
unequal  slopes  Oj  and  O?  ,  intersecting  in  a  line 
(channel)  A-). 

The  projection 
of  the  channel, 

A  on  the  map 
plane  is  taken 
as  the  Y  axis, 
the  Z  axis  is 
vertical,  and 
the  X  axis  is 
perpendicular 
to  the  YZ  plane, 
as  shown  in  (a). 

If  rj  and  r2  are 
drawn  in  the 
XY  plane  from 

a  point  ‘  on  the  Y  axis  perpendicular  to  a  contour 
line,  say  200,  on  each  side  of  the  valley,  two  right 
triangles  are  formed,  ABO  and  ACO,  and  the  fol¬ 
lowing  relations  are  true: 

tj  =  90  -  ! 


,  the  slope  of  -j  ,  is  the  angle  between  p\  and 
the  Z  axis,  rj  is  the  projection  of  pj  on  the  XY 
plane,  and  *  [  Is  the  angle  between  r(  and  the  Y 
axis.  If  a.  and  bj  are  the  angles  Ik* tween  p  j  and 
the  X  and  Y  axes  respectively,  then: 

«•<>«  a  j  =  Aj  =  sin  O  j  sin  ‘  , 
cos  bj  =  Pj  -  sin  Oj  cos  fj, 
cos  O  j  c  j  -  cos  O  j. 

o»  k\,  m,  and  cj  ,  are 
the  direction  co¬ 
sines  of  -7  !  . 

If  similarly  th 
perpendicular,  p2  , 
to  is  drawn,  as 
shown  in  figure  (c), 
the  projection  of 
p2  on  the  XY  plane 
is  r2,  angle  a2  is 
between  p,  and  the 
negative  end  of  the 
X  axis,  angle  b2  is 
between  p2  and 
the  Y  ax.s,  and  angle  02  is  between  p2  and  the  Z 
axis.  By  geometry,  the  direction  cosines  of  72 
are: 

A,  =  -  sin  Go  sin  £ 2  , 

(3)  p  2  =  sin  02  cos  £  2, 


(2) 


(1)  e2  =  90  -  -2 


c2  -  cos  02. 


&  =  ci  +  e2  =  180  -  (^i  +  £2), 

where  Cj  ,  e2  are  the  angles  the  contour  line  makes 
with  the  projected  channel  (Y  axis)  and  ^  are 
the  acute  angles  wi  ich  rj  and  r2  make  with  the  Y 
axis  (axil  angles),  and  is  the  irflection  angle  of 
the  contour  line. 

If  a  perpendicular  t  -r,  ,  n  b:  drawn,  as  in  fig¬ 
ure  (b), 


By  (1)  it  is  possible  to  eliminate  and  |2  from 
the  direction  cosines,  and  obtain  expressions  in 
f  j  and  e2  ,  Gj  and  02  .  Foi  n i : 


(4) 


*1  = 

sin 

°1 

cos 

el 

M  ,  = 

sin 

Ol 

sin 

6 1 

V  L  = 

cos 

Gj: 
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ami  (of  ? » 


(S) 


,  -  *  nln  U  •.  roa  <  »  , 
i,  -  -  ain  l);  Dili  (  <  , 


(IIS 


ion  O 

com  • U ,  • nln  ()  Mill  ■ 


i  -  -  i  tin  O’  . 

A,,  the  intersection  of  :  ,  Anti  r-  must  lie  in  the 
V Z  plane,  because  by  definition  the  Y  axis  is  the 
projection  of  A  The  direction  cosines  of  the  Y 7. 
plane  are: 


(G) 


11  YZ  0 
*'  17  °  °- 


sm  ■(  *  cot  o 


sin 


ros  •  > 


( 1  -11  sin  < 


tan) 
tan  0 1 


The  direction  numbers  of  A ,  can  also  lx1  ex¬ 
pressed  In  terms  of  the  direction  cosines  of 
and  the  YZ  plane,  and  by  analogy: 


Since  Ails  perpendicular  to  both  p\  and  the  X  axis, 
then  the  direction  numbers  of  A-,  can  t>o  found  by 
evaluating  the  determinants: 


Pi  I'l 

1  PI 

•'l 

Mvz  "yz 

0 

0 

pi 

X1 

cos  0,  sin  0,  cos  e 

PY7. 

*Y7. 

0  1 

(7) 


(8)  m-,  = 


(9)  n3  = 

=  -sin  0  jSin  e , 
(Olmsted,  1947,  p,  23). 

The  direction  cosines  of  A3are: 


=  0. 


cos  O,  , 


*1 

Pi 

sin  Ojcos  e. 

sin  OjSine, 

*  YZ 

P  YZ 

1 

0 

115) 


sin  < 


tan  > 
tan  0-. 


Assuming  that  e,  and  to  will  always  be  acute, 
which  is  invariably  true  of  natural  channels,  only 
the  positive  signs  can  be  used  in  (14)  and  (15).  From 
(14): 


(16) 


tan  e 


tan  y 


tan  Oj  cos  tj  ’ 


and  from  (15) 

(17)  tane0  = 


tan  y 


tan  02  cose  0  ’ 

Adding  (16)  and  (17),  and  solving  for  tan  y , 

tany  tany 


tan  e  j  +  tan  e , 


tan  Oicosei  tan  Oncos  ej 


(10)  =  0, 

cos  O  i 

(11)  p3  =  - — —  > 

+  (cos2  03 +sin2Oj  sin2e  i  )  ^ 

-  sin  O.sintj 

(12)  i>3= - - - 

+  (cos  2  O  i  +  sin  2  G  i  sin2  e  l  )  /J 


tan  y 


(tant  ] +tane  o )  tan  Oitan  Oncost ,  coseo 


tan  Oo  cos £2  +  tan  Oi  cose  i 


(18) 


tanei  +  taneo 
tany  =  - = - 

cot  O,  cot  Oo 

- ±-  +  - - 

cose  j  cose0 


Since  p3  is  the  cosine  of  the  angle  between  A3  and 
the  Y  axis,  cos'1  p3  is  the  gradient  of  the  channel, 

'/• 

cos  Oj 

(13)  cos  y  =  p 3  - — - - — - — - 

±  (cos20  j+sin  203  sin  2e  |  )  ! 


Equation  (18)  is  the  general  expression  of  the 
relations  among  channel  gradient,  valley-side 
slopes,  and  the  components  of  the  inflection  angle 
e3  and  e  2  ,  where  e  \  +  e  2  =  &  ■ 

In  the  special  case  where  Oj  =  0?,  that  is,  the 
two  valley  sides  have  the  same  slope,  c  3  t  ->  = 
&  ,  and  (18)  can  be  simplified  to: 


Squaring  and  simplifying, 


(19) 


tany  =  sin 


~T 


tan  O. 


/ 


II  i  in  the  o  tu-dral  an^lt-  twt*cen  th<-  plam-n 
• ,  ami  .  U.  relation  anions  (  .  <  ■  .  •  -  ,  Oj  ,  ami 
IK  if*  easily  obtained.  Hy  a  theorem  In  solid  analy- 
Ur  ^ eon. t  ( ry, 

(20  I  tut  ^  ^  K ,  •  i.  |  •  l  |  1'  -  , 

where  ,  ainl  ;  art*  I K«*  direction  cosines  of  the 

planes.  Front  (•!)  ami  (5), 

ros  (  =  -  sin  0[Sln  O-cos  t  jcosej  -sin 

OjSln  O-sliKjSlix-  .  cos  0|('os  0-. 

liy  collecting  terms  and  converting  to  functions  of 


<  i  •  t  -  .  and  O;  •  0-,  the  following  expression 
obta  ined. 

cos  (  -  costO|  •  O-l-cos’  i  ,(<  j  .  *  -  I icoslO  | -  O 

sin '  ‘  ,(>  j  •  <  , ), 
or 

(21)  cos  (  cos  (Oj  •  O-  )-cos'  ';d/«cos(0|  -  O-)- 

sin'  ‘;ip. 

In  the  case  that  Oj  =  0->,  (21)  can  lx>  simplified  t 

(22)  cos  C,  =  cos‘0  -  cos«/fsln:0. 


Appendix  II 


INDEX  OF  STUDY  AH  FAS 

The  quadrangles  listed  below  correspond  to  numbered  squares  within  each  state,  as 
given  on  the  Index  maps. 

Arizona 

1.  Cl.lnle  Badlands  (Special  Map) 

2.  Prescott 

3.  Saguaro  National  Monument  (Special  Map) 

4.  Sonolta 

5.  Mt.  Hughes 

6.  Harshaw 

Colorado 

1.  Nederland 

2.  Black  Hawk 

3.  Morrison 

4.  Starkville 

New  Mexico 

1.  Jar os a 

2.  Cerro  del  Grant 

3.  Bland 

4.  Santa  Fe 

5.  Paxton  Springs 

6.  Sedillo 

7.  Twin  Sisters 

8.  Allie  Canyon 

9.  Hurley  East 

10.  Big  Earro  Mountains 

Utah 

1.  Lehi 

2.  Sprir.gville 

3.  Spanish  Fork  Peak 
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Appendix  III 


DF1HIVATION  Of  ItKl.ATION  BETWEEN  01  AM K- 
TKH  OF  IMPRINT  AND  VOLUME  OF  IMPRINT, 
AS  USED  IN  SOIL  .STRENGTH  TEST 

The  Imprint  of  the  shot  on  ^  planar  surface  is  a 
spherical  segment  with  depth  h,  diameter  d.  and 
volume  V.  The  radius  of  the  shot,  r,  Is  0.188  feet. 
An  expression  for  V  In  r  and  d  must  lx-  found. 


o 


Is  the  volume  of  a  segment  txdow  a  lira  at  X;  ,  with 
a  radius  of  base  ♦  y,  ,  The  Integration  of  (-1)  Is 
easily  done:  let  u  r*  -  y*.  then  y-  r’  -  u, 
and  du  =  -2  y  dy. 


X 


X  =  .  (r2-y2)S 


1  y  I  >'  ?(r?  ->  <-2  V  dy)  =  y  |  (r  ■*  -u)u’  ‘du. 

By  usual  methods: 

(5)  I  =  .7(r  2u'’  -  1/3  u‘‘ )  +  C, 

and  convert  ini;  back  to  y’s 

(G)  I  =  -  [  r  2{r  2  -y  2  )'J  -  1/3 (r  2-y  2)‘‘]  +  C. 

Then  from  (4), 

V  =  I"  r 


[_r2(r2-y2,S  -  l/3(r2-y2f‘ 


yj 


(a) 

Figure  (a)  is  a  semi-circle  with  a  center  at  0,  and 
radius  r.  Since  only  the  portion  in  the  positive  side 
of  thr  Y  axis  is  considered,  its  equation  is 

(1)  X  =  +  (r2-y?)‘\ 

For  a  circle, 


=  r,  {  2  /3  r^  - 


r2(r=-y^  -l/3(r  2-y  j)1’ 
(7)  V  =  2/3  r,  r^-3r2(r?-y])'T  |  (r2-y2)\ 


or 


Since  the  relation  between  diameter  and  volume  of 
the  imprint  is  the  main  concern,  let 


(2)  dX  =  -  . 

X  (r~-y2)'’ 

and  the  element  of  volume,  dV,  is 

(3)  dV  -  uy-  dX  = - ^ 

(r2_y2)'l 

since  a  sphere  is  being  considered,  which  is  ob¬ 
tained  by  rotating  the  circle  about  the  X  axis.  By 
integration, 


d  =  2  y, 
d 


then 

(8)  V  =  2/3  nr  3 -  7rr2(r  2-  +  |(r2- 

For  the  particular  shot  used,  r  =  0.188  feet,  and 
(8)  becomes 
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*8 

'  V  -  0.0138  -  0.1 10?  .0352  -  ^  . 

*  1  0-J7  .03  52  -  j  ) 

In  »if  .soil  :»i renjjili  IcM.  i),«.  Hlrcnxth.  X 


(10) 


CO 

V 


ttluatiuii  (9)  ran  !*•  unt.l  to  n.tntrui  i 
in  icrnu  at  d.  (Set-  K.^urc  C). 


K'rapli  ctf  ,s 


Appendix  IV 


MAPS  AND  DRAINAGE  NETS  OF  STUDY  AREAS 


:*car  i  4ttirfnn  Arir«c-» 


Chtnlr  lUUUmln  North  t  fourth  ofilcr) 
Chtr.lr  lUiJiA/nb  South  iThtnl  -ftlrr) 


**  jj  .3*  ^ 

Ccfr.®'cf>  A'\icr*C 

!,,*«  '  i  ^  ;  f 


rr\ 


I  mile 


Mt.  Hughes  Quadrangle,  Arizona 


Tributaries  to  Harshaw  Creek 


S  1/2  17,  T  22  S,  R  16  E  (Fourth  order) 
Nt  1/4  2U,  T  22  S,  R  16  E  (Third  order) 
W  1/2  16,  T  22  S,  R  16  E  (Fourth  order) 


Harshaw  Quadrangle,  Arizona 

Finley  and  Adams  Canyon  1,  II  (Fourth  01  del) 
Sycamore  Canyon  I,  II  (Fourth  order) 


MCprUlwSwaH 


Cerrn  ill*  I  Gr.ml  Quad  rankle  New  Mexico 


Ci*rr<i  Pavo  Basin  (Third  order) 


•  y/ 


Bland  Quadranttle,  New  Mexico 

Above:  Paso  del  Norte  Basin  (Fourth  order) 
Below:  Peters  Dome  Basin  (Fourth  order) 


Hurley  East  Quauranjtle.  New  Mexico 

V.  1  2  22.  T  18  S.  H  12  W.  I.  II  (Tlnrd  order) 

Martin  Canyon 

SI  2  15.  T  18  S,  R  12  W  (Fourth  order) 

E  1-2  22.  T  18  N.  H  12  V,  (Third  order) 


W4« 


I  mile 


f 


1  Mf- 


ZC'CLi  1\ 


Jar  >i«j  Mcitfo 

Al»nr  Hairpin  |ljmn  Third  order l 
IHir*  IViImi  Banin  i  Forth  order 


I  m  v 

Pant  >n  Spring"  Quadrangli  New  Wr\u  < 
Quart/  Hill  !U»m  (Hllh  order) 


Santa  Fe  Quadrangle,  New  Mexico 


Above:  Cairdno  Basin  (Fourth  order) 
Below:  Tano  Basin  (Fourth  oruer) 


Seel i  1 1 n  Quadrangle,  New  Mexico 


Cedro  Peak  Basin  (Fourth  order) 


i  mflgi  '  W 


Twin  Sisters  Quadrangle.  New  Mexico 


Twin  Sisters  Quadrangle.  New  Mi'xim 


Bear  Canyon  1,  n,  III,  (Fourth  order) 


Wilson  Creek  I.  II.  III.  IV  (Fourth  order) 


Twin  Sisters  Quadrangle,  New  Mexico 


Little  Cherry  Creek  (Fifth  order) 
Mill  Creek  (Fourth  order.) 


2  000  ft 


■  ■rifc-r 

r.’n  nl  > r  ilt  r 
Crii'V  T  i r  1  ird'  T 


Spanish  Fork  Peak  Quadrangle.  Utah 


Polo  Canyon  (Fifth  order) 


A|>j>cii'li\  \ 


I  \U)  \AKIAIU  i  ('l)HHI. I  A  I  ION’  AM)  MKl.lt f  ,S:.K)N  AMU  YSI  s 
A  (.'■  >1  rrlaiimi  lallt  s  »l  illations  analwi'd  In  g  t  ouprd- data  incthisls. 

1  I  "(I,  Dtainagi  Ivnsto,  \s.  P-F  Index. 

2.  Log,  Dr. i mace  IVnsttv  vs.  Infiltration  Capai  lly. 

3.  L<ft;r  Drainage  Density  vs.  peri fnt  Hire  Area. 

•1 .  l.og^  Drainage  Density  vs.  Valley-snJe  .Slojw. 

5.  Log  Channel  I'riqurm  y  vs.  Infiltration  Capacity. 

G.  Ian;  Channel  Frequency  \s.  P- F  Index. 

7.  L< *g  Channel  Frequency  is.  Percent  Hare  Area. 

H.  Log  Channel  Fr .  ;  vs.  Ian;  Wet  Soil  Strength. 

9.  Lin;  Relative  Rcnef  vs.  Ian;  Channel  Frequency. 

10.  Lot;  Frequency  of  First-order  Channels  vs.  Lin;  Channel  Frequency. 

11.  Log  Ruggedness  Numlier  Corrected  for  Order  vs.  Log  Heh.tive  Relief. 

12.  Log(.  Density  of  First-order  Channels  vs.  Log,,  Drainage  Den  itv. 

13.  Lot;,.  Drainage  Density  vs.  Channel  Frequency. 

14.  Valley-side  Slope  vs.  Infiltration  Capacity,  Fourth-order  Basins. 

15.  Valiey-sicle  Slone  vs.  Lot;  P-E  Index. 

1G.  Valley-side  Slope  vs.  Wet  Soil  Strength. 

17.  Valley-side  Slope  vs.  Log  Channel  Frequency. 

18.  Valley-side  Slope  vs.  Log  Ruggedness  Number  Corrected  for  Order. 

19.  Valley-side  Slope  vs.  Log  Relative  Relief. 
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3.  7502 


1.5078  2.0654  2.6190  5  1746  3  7302 
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H  Natter  I *1  < >t r»  .11  w I  H«'r; r«- m sl<>"  i  <d  1  1  Relations .  I'ng reopen  Data. 

1  Drainage  iVnstti  is,  p- h  |mle\. 

Drainage  Density  \s.  Runoff  Intensity. 

3  Channel  Pie<pien<\  '.s.  percent  Hare  Aren.  Kach  point  represent.-!  one  major 
i  «a  sin. 

•1.  Regiession  of  Cltannel  frequency  on  Drainage  IX-nstty.  Kadi  point  repre¬ 
sents  one  thud-order  Imstn. 

rj.  Reg res s ion  of  \  alley-side  Slope  on  Infiltration  Capacity. 

G.  Regression  of  Vallev-side  Slope  on  Runoff  Intensity. 

7.  Regression  of  \  alley-side  Slope  on  Relative  Relief. 

H.  Regression  of  \all<  v-sule  SR  pe  on  Ruggedness  Numt>er  Corrected  for  Order. 
It.  Regression  of  Vallev-side  Slope  on  Wet  Soil  Strength. 

10.  Regression  of  \  alloy-side  Slope  on  Ratio  of  Mean  January  Precipitation  to 

1  12  Mean  Annual  Precipitation. 

11.  Regression  ol  Drainage  Density  on  Percent  Bare  Area. 
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Appendix  \  1  A 


A’.AtYs*so»  \AHtANCl 

}..'  «•(  i  I  lith  s  dr  Aitv-n'r  di-nnit-.  All  fourth-  >r  1«-  r  tasms,  and  thtnl-  >rd<  i  I  as  lies  tf  Uu\  arc  not 

in<  juded  in  a  J  u r t h •  <> j  <!«•  r  oaMjii  a r <*  Koyfi  cs  of  data  for  this  stud\. 


()l  'u  r\rd  Draina^r  1  H-ris dies  (■rouped  I > l.ith>>l>>i;  • 


>hni< 

Clastic 

Ltniosti  in 

179.2 

24. 1 

14.4 

171.2 

20.6 

10.0 

3  50.4 

1  1.2 

24.4 

13.1 

V  ) 

15.1 

N  2 

x  c. 

16.2 

1  7  r. 

12.11 

X  12.2 

11.4 

11.1 

33.5 

17.1 

17.1 

14.7 

I4..H 

219.7 
N  13 
X  16.9 


Schist 

(ir.Tiit  ic 

Ac  id  \  oh  amt 

80.8 

17.4 

1  5.8 

86.6 

19.5 

14.1 

1  12.7 

15.2 

13.8 

22.9 

17.5 

14.9 

10.4 

16.4 

14.4 

313.4 

11.0 

9.9 

5.7 

4.4 

«’  5 

7.5 

13.1 

7.5 

15.0 

-  62.7 

5.0 

11.9 

14.8 

9.2 

13.7 

13.8 

14.9 

12.9 

17.1 

13.2 

15.9 

12.1 

36.9 

11.2 

58.1 

13.2 

39.3 

12.6 

32.2 

13.2 

39.9 

13.1 

32.9 

12.0 

438.4 

12.3 

10.4 

N  =  21 

286.5 

X  =  20.9 

N  =  23 

X  =  12.5 
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19  t  43  9  1 

1  329  74 

(  1  a  it  ti 
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1*  I 

lo  o 

•  7 .  ft 

4  173  47 

3.7  12  33 

3  0  3  0 

t  ,r  an  it  ii 

2  1 

438  4 

20  9 

•  J  i 

12.  <47  18 

1  52  .  2 

Irt  »  7  '.> 

I  d  (».«  n|.  air 

> 

.’4  4 

12.2 

-12.5 

3  07  3  0 

2 97.0ft 

9  0  8 

\<  id  \  ' >ii  .inn 

2 1 

.*86.  5 

12.5 

•  12.3 

3 .  o‘.<3 . 7  J 

3,508.79 

5.08 

I •  tain 

00 

I.UJJ.H 

.  . 

1  51.9 

1 09.90b.Hh 

97  705.51 

... 

I  \  i  s  40 .  J  .  -8  H 

Analysis- of-  \  artniu  «•  I  able 


Source 

of  Variation 

S.  S. 

d.f. 

M.  .5. 

1 

Among  i 

olumn  means 

57.371.03 

5 

1  1.47  4.2  1 

56.7 

Within  t 

olutnns 

12  140.37 

00 

202.34 

Total 

69,311.40 

G5 

1,069.41 

V  „  ,  (5,  GO)  3.339 

Thorefire,  a  significant  difference  among  means  of  drainage  density  for  each  lithology  exists  at  the 
.01  level  of  significance,  and  also  at  the  .001  level. 


Contrasts*  judged  significant  at  the  .05  level. 


(Granitic,  clastic,  volcanic,  and  limestone)  vs. 

(shale  and  schist)  (-97.97  -  O'  -58.85) 

Shale  vs.  Schist  (71.58  -  O  •  153.46) 

Contrasts  judged  not  significant  at  the  .05  level: 

Clastic  vs.  granitic.  (-21.24  O'  29.20) 

Limestone  vs.  granitic.  (-44.92  O'  27.52) 

Limestone  vs.  volcanic.  (-35.33  •  O'  36.81) 

Clastic  and  granitic  vs.  limestone  and  volcanic 

(-5.98  •  O'  19.82) 


Effect  of  lithology  on  infiltration  capacity:  Individual  measurements  of  infiltration  capacity  grouped  ac¬ 
cording  to  lithology,  regardless  of  basin  or  climate,  are  sources  of  data.  (Individual  infiltration  meas¬ 
urements,  grouped  by  basin,  are  given  in  Table  7.) 


♦A  contrast  is  a  <ii I lerence  in  means  or  averages  of  means  of  several  groups.  In  Schelde's  notation,  the  true  (popula¬ 
tion)  contrast  is  denoted  by  W,  and  the  estimated  contrast  calculated  I  rom  sample  data  is  denoted  by  t3.  His  notation 
is  used  in  this  section,  and  must  not  be  contused  with  the  identical  notation,  W,  used  lor  valley-side  slopes  in  the 
previous  sections. 
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S.  S. 

d.f. 

M.  S. 

1 

Among  t 

olumn  means 

67.89 

5 

13.58 

7.43 

Within  i 

oluinns 

249.4 

135 

1.85 

Total 

317.29 

140 

F  1)f,  j(f>.  135)  *1.42. 

A  significant  difference  exists  among  the  infiltration  means  for  the  various  types  of  lithologies  at  the  .001 
lev  el  of  significance. 

The  contrast  (granitic  and  limestone)  vs.  (clastic  and  volcanic >,  formed  from  the  groups  with  highest 
and  intermediate  infiltration  capacities,  was  judged  significant  at  the  .05  level  (0.450  •  O  •  2.042.  in 
Scheff/'s  notation).  The  contrast  (clastic  and  volcanic)  vs.  (shale  and  schist),  formed  from  the  groups  with 
intermediate  and  lowest  infiltration  capacities,  was  judged  not  significant  at  the  .05  level  (-0.780  •  O  • 
2.100). 

Effect  of  lithology  on  valley-side  slope:  Individual  field  measurements  of  slopes,  separated  according 


to  lithologic  type  irrespective  of  basin 

boundaries, 

are  source: 

;  of  data. 
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Totals  880 
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(XLX, ,)  2  N  396,555.17 

Analysis-of 

-Variance  Table 

Source  of  Variation 
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d.f. 

M.  S. 
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11,983.48 
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Within  groups 

110,298.35 
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Total 

122,281.83 
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nel  condition  arc  availalilc. 
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j  > 

Total  sum  of  products  =  -12.77. 

Within-groups  sum  of  products  -  -12.03. 

Among- means  sum  of  products  =  -0.74. 

Total  sum  of  squares:  for  X’s  11.92,  for  Y’s  24.86. 
Among-means  sum  of  squares:  for  X’s  1.4,  for  Y’s  0.42. 
Within-groups  sum  of  squares:  for  X’s  10.52,  for  Y’s  24.44. 


el 


Vats.  !  5<|uatrs  4li»al  (f'.f  rej?r«  n »  on  lsn<- 

(.  I  2  77  I 

I  ml  s.  h.  JCbC  •  11  l rt 

I  l  >2 

(-  1  2  03  I 

V.  ttlun-v;  fmi|)h  > .  .V  24  ^  4  -  JO  |0.0« 


I  In-  iltspi  t  iion  u(  I  hr  V  \alu<:>  altout  tin'  rt*i;rc!t!ti<)ti  line  willi  the  coefficient  t>,  (total)  in  therefore  11.18 
I  hi-  dispersion  of  the  'i  \ .dues  ;iUml  the  rt'itri'vitoii  line  with  coefficient  h.  (within)  Is  therefore  10  68. 
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that  tlie  reduction  is  not  (treat  enough  to  !*•  judged  significant. 
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Solution  of  Normal  Equations  Relating  Drainage 
Density  to  Soil  Properties  and  Climate  Elements 
by  the  Doolittle  Method. 
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Morrison  Quadrangle 
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100 


South -facing •  shrub  live  oak  chaparral, 
manzanita  . 

North-facing  emory  oak,  sllverleaf  oak,  netleaf 
oak,  Arizona  white  oak.  alligator  juniper,  manzanita, 
madrono,  Mexican  pi  non,  Mexican  blue  oak, 
abundant  grass. 

Below  1800:  catclaw  acacia. 

1800-5800:  shrub  live  oak  chaparral,  manzanita, 
ocotilla,  yucca. 

Above  5800:  Mexican  blue  oak,  Arizona  white 
oak,  alligator  juniper,  Mexican  pinon. 


Arizona  white  oak,  shrub  live  oak,  one-seed  56 

juniper,  alligator  juniper,  hairy  mountain  ma¬ 
hogany,  singleleaf  pinon,  ponderosa  pine  (small), 
cliffrose,  cholla. 

South-facing:  mainly  grass.  38 

North-facing:  emory  oak,  Arizona  white  oak,  3 

Mexican  blue  oak,  century  plant,  jumping  cholla, 
manzanita,  alligator  juniper. 

Saguaro,  mesquite,  paloverde,  ocotilla,  jumping  83 

cholla,  staghorn  cholla,  creosote  bush,  prickly 
pear,  fishhook  barrel  cactus,  brittle  bush,  Christ¬ 
mas  cactus,  sparse  bunch  grass. 


Lodgepole  pine,  blue  spruce,  Douglas-fir. 


Grass,  sage,  yucca,  poison  ivy,  wild  rose,  juni-  40 

per,  prickly  pear,  chokecherry,  hairy  mountain 
mahogany,  gooseberry,  bigtoo^h  maple,  barrel 
cactus,  Russian  thistle. 

Birchleaf  mountain  mahogany,  prickly  pear,  50.2 

grass,  sage,  yucca,  common  chokecherry  (on 
north  slopes  and  gullies),  barrel  cactus,  sun¬ 
flowers,  gooseberry. 


Aspen,  lodgepole  pine,  ponderosa  pine,  limber 
pine,  white  fir,  blue  spruce,  mountain  common 
juniper,  grass,  Douglas-fir,  willows. 
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Lodgopolr  J>tne,  limber  pine,  while  fir.  blue  27.8 

spruce,  aspen,  mountain  common  Juniper,  sage, 
bunch  grass,  klnnlkinnlk 


PI  non -juniper  woodland. 


Shrub  live  oak,  Gambel  oak,  alligator  juniper, 
ponderosa  pine,  yucca,  barrel  cactus,  bunch 
grass,  prickly  pear,  eholla,  Douglas-fir. 
South-facing:  shrub  live  oak  chaparral. 
North-facing:  ponderosa -pine  forest. 

South-facing:  shrub  live  oak  chaparral,  Gambel 
oak,  alligator  juniper,  yucca,  barrel  cactus, 
bunch  grass,  prickly  pear,  eholla. 

North-facing:  ponderosa  pine,  Douglas-fir,  shrub 
live  oak,  Gambel  oak,  alligator  juniper. 

Ponderosa  pine,  shrub  live  oak,  Gambel  oak, 
Douglas-fir,  walnut,  alligator  juniper,  prickly 
pear,  bunch  grass,  eholla. 


Ponderosa  pine,  shrub  Gambel  oak,  white  fir, 

Douglas-fir,  aspen,  Rocky  Mountain  juniper, 
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Ponderosa  pine,  shrub  Gambel  oak,  aspen,  white  25.8 

fir,  Douglas-fir,  Rocky  Mountain  juniper,  birch¬ 
leaf  mountain  mahogany,  mountain  common  juni¬ 
per,  pinon. 
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A  I'lflJ  mtaigrtmcft!!1 

1)  Within  b»olno 

S  Hf*0.  it  22,  *2  -  lf»d*T.  n  l  •:  ° 

Within  lltho|o£ir« 

N  HhO.  k  G.  »•  126.19,  »  II  .2" 

II  SU|i  mnnurrmrnbi 
Within  Ixmintt 

N  113.  k  2  1,  »•’  3111.  s»  5.6  0 

2  Portent  Haro  Area 

Within  basins  (Individual  Hold  measurements) 

N  *  222,  k  -  21.  s-  206.65.  s  11.1% 

3.  Infiltration  Capacity 

A.  Within  basins  (individual  field  measurements) 

N  122,  k  22.  s2  1.16,  s  1.08  in  hr. 

D.  Within  lithologies  (Individual  field  measurements) 

N  =  122,  k  =  6,  s2  1.31,  s  1.11  ln/hr. 

•I.  Floughness  Number 

Within  lithologies 

N  22,  k  =  6,  s2  =  628.8,  s  25.08  inches 

5.  Drainage  Density 

A.  Within  homogeneous,  contiguous  areas 

N  100,  k  =  24,  s2  =  68.09,  s  =  8.25  per  mile 
R.  Within  lithologies 

N  =  106,  k  =  6,  s  2  =  407.33,  s  =  20.18  per  mile 

6.  Relative  Relief 

Similar  basins  within  homogeneous,  contiguous  areas 
N  =  106,  k  --  25,  s2  =  3.24,  s  =  1.80 

7.  Ruggedness  Number 

Similar  basins  within  homogeneous,  contiguous  areas 
N  =  108,  k  =  26,  s2  =  0.596,  s  =  0.772 

*N  =  Total  number  of  observations 
k  =  Number  of  basins,  lithologic  classes,  etc. 
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Pine  duff,  R 

u  .J 

i  n 

10 

Grass;  leaves  R 

1.8 

6 

6 

Grass;  U 

Pine  duff,  R 

1.7 

2.0 

3 

10 

Pine  duff,  R 

"  n  n 

3.6 

12 

*  *  » 

1.7 

11 

*  „  „ 

1.5 

9 

”  ft  r> 

1.2 

1.4 

18 

9 

Grass,  U 

1.7 

8 

2.1 

1.9 

11 

9 

Pine  duff,  R 

ft  V 

2.0 

1.8 

1.8 

3.3 

3.5 
1.9 

6.5 


11 

6 

5 

8 

12 

7 

13 


Bare,  U 

ft  rt 

"  n 
»  rt 

Pine  duff,  R 
Bare,  U 

ft  it 


PAXTON  SPRINGS  QUADRANGLE 
Quartz  Hill  Basin  (granite) 
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J  MH  1.  '  |>» 

UJUlJi 

r. nt  r 

m  t  •  >p  i;  i  >» 

type  <>»  si  it> 

CAP 

pi  at  r* 

,V\NTA  |  |.  Atm  VMU.J. 

1  attains)  =  Iii  unutti) 

:  o 

<, 

liar i  ,  l 

i  t 

10 

(iriro,  l 

:  j 

i 

Hare  t' 

Taft..  lUolft  (allutluti  > 

i  <• 

:  * 

Hare.  1 

l  i 

i  »> 

l  . 

i 

1 

Sparse  grass.  1 

MIDII  l.i  *  Wl  ADIt  \ Ml  1.1. 

Coir  Pei)  Il.mln  (sandstone,  shale) 

2  l 

10 

Bar.  .  1 

i 

1 1 

Hare.  U 

:  o 

1.7 

Pine  lull ,  It 

i  i 

• 

»'> 

1 

i  i 

_•  \i 

i 

1  I 

. 

i  , 

f> 

* 

2..1 

JS 

"  "  " 

l  t.ih 

I. Kill  QUADRANGLE 

Dry  Hollow  (<|uart/.lte,  volcanic) 

<>.o 

26 

Grass,  1 

■i  . 

•i  ■ 

20 

Bare,  L 

_ .  a 

0.0 

i) 

. 

if. 

12 

- 

2.1 

21 

Grass,  U;  oak  lit. 

SPRINGVILLK  QUADRA NG l. E 

Hobble  Cr.Trib.  (sandstone) 

2.3 

•>  ) 

Grass,  U 

l  G 

11 

Bare,  U 

Spring  Creek  (limestone) 

2.9 

1  1 

Grass,  U 

2.6 

30? 

"  " 

•1.3 

28 

0  8 

20 

Leaves,  grass,  U 

1.8 

2G 

Grass,  U 

2.7 

21 

Leases,  grass,  U 

1.0 

19 

Leaves,  weeds,  U 

1.1 

*» 

Grass;  leaves,  R 

Additional  measurements  on  residual  graiite  soils: 

Arizona 

HARSHAW  (PATAGONIA  MTNS.)  QUAD,  (granite) 

2.3 

? 

1.9 

Bare 

i)  .  / 

5.3 

Grass,  R 

IRON  SPGS.  (SIERRA  PRIETA)  QUAD,  (granite) 

2.5 

Pine  duff,  R 

2.0 

WWW 

1.8? 

Rare.  U 

2.0 

Pine  duff,  R 

2.6 

? 

6.2 

Pine  duff,  R 

•1 .5 

Bare,  U 

1.0 

Pine  duff,  R 

TABLE  7  (1) 


11  \SIN 

INfcILT 

SLOPE  cfc 

TYPfc  nfc  SI  It  fc  ACfc 

CAP 

PL\T  (”) 

\i  1  t.im.i 

KIHK  LAND  (SIERRA  I’lUKTA)  QUAD.  (granite  ) 

.*..1 

2.1 

Pin.  .lull ,  :i 

Ni-m  Mi'xicu 

HIC.  lit  UHL)  MTNS.gUAUHANGI.fc;  (uraniU-) 

•1 ... 

Li’avfji,  It 

7.M 

Bar.-,  L 

(.0 

”  ** 

3.7 

"  ** 

8  .3  • 

8.3- 

Pino  .lull.  It 

fc'ORT  DA  YARD  (PINOS  ALTOS)  QUAD,  (granite) 

3 .3 

Bari-,  U 

1.8 

1.3 

TWIN  SISTERS  (PINOS  ALTOS  MTNS.)  QUAD. 

2.1 

Bari-,  U 

(granite*) 

.  "i 

